Changeset 45 for codes/icosagcm/trunk/src/timeloop_gcm.f90
- Timestamp:
- 07/28/12 02:56:28 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/timeloop_gcm.f90
r32 r45 33 33 REAL(rstd),POINTER :: du(:,:), dum1(:,:), dum2(:,:) 34 34 REAL(rstd),POINTER :: dtheta_rhodz(:,:),dtheta_rhodzm1(:,:),dtheta_rhodzm2(:,:) 35 REAL(rstd) :: dt 35 REAL(rstd) :: dt, run_length 36 36 INTEGER :: ind 37 37 INTEGER :: it,i,j,n … … 44 44 dt=90. 45 45 CALL getin('dt',dt) 46 dt=dt/scale_factor47 46 48 47 itaumax=100 49 48 CALL getin('itaumax',itaumax) 49 50 run_length=dt*itaumax 51 CALL getin('run_length',run_length) 52 itaumax=run_length/dt 53 PRINT *,'itaumax=',itaumax 54 dt=dt/scale_factor 50 55 51 56 write_period=0 … … 101 106 102 107 DO it=0,itaumax 103 PRINT *,"It No :",It 108 PRINT *,"It No :",It," t :",dt*It 104 109 105 110 CALL guided(it*dt,f_ps,f_theta_rhodz,f_u,f_q)
Note: See TracChangeset
for help on using the changeset viewer.