Changeset 63 for codes/icosagcm/trunk
- Timestamp:
- 07/31/12 21:43:20 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/timeloop_gcm.f90
r50 r63 38 38 INTEGER :: matsuno_period 39 39 INTEGER :: itaumax 40 40 REAL(rstd) ::write_period 41 INTEGER :: itau_out 42 41 43 dt=90. 42 44 CALL getin('dt',dt) … … 58 60 IF (TRIM(scheme)=='leapfrog') matsuno_period=itaumax+1 59 61 62 write_period=0 63 CALL getin('write_period',write_period) 64 write_period=write_period/scale_factor 65 itau_out=FLOOR(.5+write_period/dt) 66 PRINT *, 'Output frequency (scaled) set to ',write_period, ' : itau_out = ',itau_out 67 60 68 CALL allocate_field(f_phis,field_t,type_real) 61 69 … … 95 103 DO it=0,itaumax 96 104 PRINT *,"It No :",It," t :",dt*It 97 105 IF (mod(it,itau_out)==0 ) THEN 106 CALL writefield("q",f_q) 107 ENDIF 108 98 109 CALL guided(it*dt,f_ps,f_theta_rhodz,f_u,f_q) 99 110 CALL caldyn(it,f_phis,f_ps,f_theta_rhodz,f_u, f_dps, f_dtheta_rhodz, f_du) 100 111 CALL advect_tracer(f_ps,f_u,f_q) 112 101 113 102 114 SELECT CASE (TRIM(scheme))
Note: See TracChangeset
for help on using the changeset viewer.