Changeset 45

Timestamp:

07/28/12 02:56:28 (12 years ago)

Author:

dubos

Message:

Introduced run_length
Increased time step in run.def for test case 3 (gravity wave)

Location:

codes/icosagcm/trunk

Files:

• param_sets/dcmip2012/test3/run.def (modified) (1 diff)
• src/timeloop_gcm.f90 (modified) (3 diffs)

codes/icosagcm/trunk/param_sets/dcmip2012/test3/run.def

@@ -33 +33 @@
 
 # timestep : real (default=90s)
-dt = 45.
+dt = 90.
 
 # number of timestep (default 100)
-itaumax = 10000
+# itaumax = 10000
+
+# run length : real (default=dt*itaumax)
+run_length=450000
 
 # output field period : integer (default none)

codes/icosagcm/trunk/src/timeloop_gcm.f90

@@ -33 +33 @@
   REAL(rstd),POINTER :: du(:,:), dum1(:,:), dum2(:,:)
   REAL(rstd),POINTER :: dtheta_rhodz(:,:),dtheta_rhodzm1(:,:),dtheta_rhodzm2(:,:)
-  REAL(rstd) :: dt
+  REAL(rstd) :: dt, run_length
   INTEGER :: ind
   INTEGER :: it,i,j,n
@@ -44 +44 @@
   dt=90.
   CALL getin('dt',dt)
-  dt=dt/scale_factor
   
   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
@@ -101 +106 @@
   
   DO it=0,itaumax
-    PRINT *,"It No :",It
+    PRINT *,"It No :",It,"   t :",dt*It
 
     CALL guided(it*dt,f_ps,f_theta_rhodz,f_u,f_q)