Changeset 1577 for LMDZ5/trunk/libf/dyn3d/leapfrog.F
- Timestamp:
- Oct 20, 2011, 5:06:47 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3d/leapfrog.F
r1529 r1577 224 224 1 CONTINUE 225 225 226 jD_cur = jD_ref + day_ini - day_ref + int (itau * dtvr / daysec) 227 jH_cur = jH_ref + & 226 jD_cur = jD_ref + day_ini - day_ref + & 227 & int (itau * dtvr / daysec) 228 jH_cur = jH_ref + start_time + & 228 229 & (itau * dtvr / daysec - int(itau * dtvr / daysec)) 230 jD_cur = jD_cur + int(jH_cur) 231 jH_cur = jH_cur - int(jH_cur) 229 232 230 233 … … 372 375 ! rdaym_ini = itau * dtvr / daysec 373 376 ! rdayvrai = rdaym_ini + day_ini 374 jD_cur = jD_ref + day_ini - day_ref 375 $ + int (itau * dtvr / daysec) 376 jH_cur = jH_ref + & 377 & (itau * dtvr / daysec - int(itau * dtvr / daysec)) 377 ! jD_cur = jD_ref + day_ini - day_ref 378 ! $ + int (itau * dtvr / daysec) 379 ! jH_cur = jH_ref + & 380 ! & (itau * dtvr / daysec - int(itau * dtvr / daysec)) 381 jD_cur = jD_ref + day_ini - day_ref + & 382 & int (itau * dtvr / daysec) 383 jH_cur = jH_ref + start_time + & 384 & (itau * dtvr / daysec - int(itau * dtvr / daysec)) 385 jD_cur = jD_cur + int(jH_cur) 386 jH_cur = jH_cur - int(jH_cur) 378 387 ! write(lunout,*)'itau, jD_cur = ', itau, jD_cur, jH_cur 379 388 ! call ju2ymds(jD_cur+jH_cur, an, mois, jour, secondes) … … 652 661 ! if (planet_type.eq."earth") then 653 662 ! Write an Earth-format restart file 654 CALL dynredem1("restart.nc", 0.0,663 CALL dynredem1("restart.nc",start_time, 655 664 & vcov,ucov,teta,q,masse,ps) 656 665 ! endif ! of if (planet_type.eq."earth") … … 763 772 IF(itau.EQ.itaufin) THEN 764 773 ! if (planet_type.eq."earth") then 765 CALL dynredem1("restart.nc", 0.0,774 CALL dynredem1("restart.nc",start_time, 766 775 & vcov,ucov,teta,q,masse,ps) 767 776 ! endif ! of if (planet_type.eq."earth")
Note: See TracChangeset
for help on using the changeset viewer.