Changeset 97 for codes/icosagcm/trunk/src/time.f90

Timestamp:

08/04/12 00:13:44 (12 years ago)

Author:

ymipsl

Message:

Put conforming IO for dcmip testcase

YM

File:

• codes/icosagcm/trunk/src/time.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/time.f90

@@ -1 +1 @@
 MODULE time_mod
-
+  use prec
   PRIVATE
 
@@ -7 +7 @@
   INTEGER,SAVE :: it
 
-  PUBLIC create_time_counter_header, update_time_counter, close_time_counter
+  REAL(rstd),SAVE :: dt
+  REAL(rstd),SAVE :: write_period
+  INTEGER,SAVE :: itau_out
+  INTEGER,SAVE :: itaumax
+  
+  PUBLIC create_time_counter_header, update_time_counter, close_time_counter, init_time
 
 
 
 CONTAINS
+  
+  SUBROUTINE init_time
+  USE earth_const
+  USE ioipsl
+  IMPLICIT NONE
+  REAL(rstd) :: run_length
+  
+    dt=90.
+    CALL getin('dt',dt)
+
+    itaumax=100
+    CALL getin('itaumax',itaumax)
+
+    run_length=dt*itaumax
+    CALL getin('run_length',run_length)
+    itaumax=run_length/dt
+    PRINT *,'itaumax=',itaumax
+    dt=dt/scale_factor
+
+    write_period=0
+    CALL getin('write_period',write_period)
+    write_period=write_period/scale_factor
+    itau_out=FLOOR(.5+write_period/dt)
+    PRINT *, 'Output frequency (scaled) set to ',write_period, ' : itau_out = ',itau_out
+    
+    CALL create_time_counter_header
+    
+  END SUBROUTINE init_time
 
   SUBROUTINE create_time_counter_header
@@ -21 +54 @@
   INTEGER :: timeid, dtid
   REAL(rstd) :: dt
+  CHARACTER(LEN=255) :: time_frequency
+  
     status = NF90_CREATE('time_counter.nc', NF90_CLOBBER, ncid)
     status = NF90_DEF_DIM(ncid,'time_counter',NF90_UNLIMITED,timeid)
@@ -28 +63 @@
     status = NF90_PUT_ATT(ncid,time_counter_id,"calendar","noleap")
     status = NF90_DEF_VAR(ncid,'mdt',NF90_DOUBLE,varid=dtid)
+
+    WRITE(time_frequency,*) write_period
+    PRINT*,TRIM(time_frequency)
+    status = NF90_PUT_ATT(ncid,NF90_GLOBAL,"time_frequency",TRIM(time_frequency)//"s")
+
     status = NF90_ENDDEF(ncid)