Changeset 1382 for LMDZ4/branches
- Timestamp:
- May 6, 2010, 3:44:09 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/gcm.F
r1380 r1382 274 274 CALL iniacademic(vcov,ucov,teta,q,masse,ps,phis,time_0) 275 275 endif 276 276 !#ifdef CPP_IOIPSL 277 277 if (planet_type.eq."earth") then 278 278 #ifdef CPP_EARTH 279 279 ! Load an Earth-format start file 280 280 CALL dynetat0("start.nc",vcov,ucov, 281 & teta,q,masse,ps,phis, time_0) 282 #endif 283 ! SW model also has Earth-format start files 284 ! (but can be used without the CPP_EARTH directive) 285 if (iflag_phys.eq.0) then 286 CALL dynetat0("start.nc",vcov,ucov, 287 & teta,q,masse,ps,phis, time_0) 288 endif 289 281 . teta,q,masse,ps,phis, time_0) 282 #endif 290 283 endif ! of if (planet_type.eq."earth") 291 292 284 c write(73,*) 'ucov',ucov 293 285 c write(74,*) 'vcov',vcov … … 498 490 #ifdef CPP_IOIPSL 499 491 time_step = zdtvr 500 if (ok_dyn_ins) then 501 ! initialize output file for instantaneous outputs 502 ! t_ops = iecri * daysec ! do operations every t_ops 503 t_ops =((1.0*iecri)/day_step) * daysec 504 t_wrt = daysec ! iecri * daysec ! write output every t_wrt 505 t_wrt = daysec ! iecri * daysec ! write output every t_wrt 506 CALL inithist(day_ref,annee_ref,time_step, 507 & t_ops,t_wrt) 508 ! CALL inithist_p(dynhist_file,day_ref,annee_ref,time_step, 509 ! . t_ops, t_wrt, histid, histvid) 510 endif 511 512 IF (ok_dyn_ave) THEN 513 ! initialize output file for averaged outputs 514 t_ops = iperiod * time_step ! do operations every t_ops 515 t_wrt = periodav * daysec ! write output every t_wrt 516 CALL initdynav(day_ref,annee_ref,time_step, 517 & t_ops,t_wrt) 492 IF (mpi_rank==0) then 493 if (ok_dyn_ins) then 494 ! initialize output file for instantaneous outputs 495 ! t_ops = iecri * daysec ! do operations every t_ops 496 t_ops =((1.0*iecri)/day_step) * daysec 497 t_wrt = daysec ! iecri * daysec ! write output every t_wrt 498 t_wrt = daysec ! iecri * daysec ! write output every t_wrt 499 CALL inithist(day_ref,annee_ref,time_step, 500 & t_ops,t_wrt) 501 endif 502 503 IF (ok_dyn_ave) THEN 504 ! initialize output file for averaged outputs 505 t_ops = iperiod * time_step ! do operations every t_ops 506 t_wrt = periodav * daysec ! write output every t_wrt 507 CALL initdynav(day_ref,annee_ref,time_step, 508 & t_ops,t_wrt) 518 509 ! CALL initdynav_p(dynhistave_file,day_ref,annee_ref,time_step, 519 510 ! . t_ops, t_wrt, histaveid) 520 END IF 511 END IF 512 ENDIF 521 513 dtav = iperiod*dtvr/daysec 522 514 #endif
Note: See TracChangeset
for help on using the changeset viewer.