Changeset 413 for codes/icosagcm/trunk/src/timeloop_gcm.f90
- Timestamp:
- 06/10/16 17:49:07 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/timeloop_gcm.f90
r407 r413 31 31 32 32 CHARACTER(len=255) :: def 33 34 ! IF (xios_output) itau_out=1 33 34 CALL init_caldyn 35 36 IF (xios_output) itau_out=1 35 37 IF (.NOT. enable_io) itau_out=HUGE(itau_out) 36 38 … … 41 43 STOP 42 44 END IF 43 44 nqdyn = 1 ! one dynamical tracer = theta for the moment45 45 46 46 def='ARK2.3' … … 141 141 CALL init_sponge 142 142 CALL init_observable 143 CALL init_caldyn144 143 CALL init_guided 145 144 CALL init_advect_tracer … … 220 219 CALL trace_on 221 220 221 IF (xios_output) THEN ! we must call update_calendar before any XIOS output 222 CALL xios_update_calendar(1) 223 END IF 224 CALL write_output_fields_basic(.TRUE., f_phis, f_ps, f_mass, f_geopot, f_theta_rhodz, f_u, f_W, f_q) 225 222 226 DO it=itau0+1,itau0+itaumax 223 227 224 228 IF (is_master) CALL print_iteration(it, itau0, itaumax, start_clock, rate_clock) 229 225 230 IF (xios_output) THEN 226 CALL xios_update_calendar(it)231 IF(it>itau0+1) CALL xios_update_calendar(it-itau0) 227 232 ELSE 228 233 CALL update_time_counter(dt*it) … … 318 323 f_ps,f_dps,f_u,f_theta_rhodz,f_phis) 319 324 CALL check_conserve(f_ps,f_dps,f_u,f_theta_rhodz,f_phis,it) 320 ENDIF 325 ENDIF 321 326 322 327 IF (mod(it,itau_out)==0 ) THEN 323 328 CALL transfert_request(f_u,req_e1_vect) 324 CALL write_output_fields_basic( f_phis, f_ps, f_mass, f_geopot, f_theta_rhodz, f_u, f_W, f_q)329 CALL write_output_fields_basic(.FALSE.,f_phis, f_ps, f_mass, f_geopot, f_theta_rhodz, f_u, f_W, f_q) 325 330 ENDIF 326 331
Note: See TracChangeset
for help on using the changeset viewer.