Changeset 63 for codes/icosagcm/trunk/src

Timestamp:

07/31/12 21:43:20 (12 years ago)

Author:

ymipsl

Message:

add output for tracer

YM

File:

• codes/icosagcm/trunk/src/timeloop_gcm.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/timeloop_gcm.f90

@@ -38 +38 @@
   INTEGER :: matsuno_period
   INTEGER :: itaumax
-
+  REAL(rstd) ::write_period
+  INTEGER    :: itau_out
+  
   dt=90.
   CALL getin('dt',dt)
@@ -58 +60 @@
   IF (TRIM(scheme)=='leapfrog') matsuno_period=itaumax+1
 
+  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 allocate_field(f_phis,field_t,type_real)
   
@@ -95 +103 @@
   DO it=0,itaumax
     PRINT *,"It No :",It,"   t :",dt*It
-
+    IF (mod(it,itau_out)==0 ) THEN
+      CALL writefield("q",f_q)
+    ENDIF
+    
     CALL guided(it*dt,f_ps,f_theta_rhodz,f_u,f_q)
     CALL caldyn(it,f_phis,f_ps,f_theta_rhodz,f_u, f_dps, f_dtheta_rhodz, f_du)
     CALL advect_tracer(f_ps,f_u,f_q)
+
     
     SELECT CASE (TRIM(scheme))