Changeset 1357 for LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/gcm.F
- Timestamp:
- Apr 14, 2010, 4:03:19 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/gcm.F
r1333 r1357 73 73 #include "description.h" 74 74 #include "serre.h" 75 #include "com_io_dyn.h"75 !#include "com_io_dyn.h" 76 76 #include "iniprint.h" 77 77 #include "tracstoke.h" 78 #ifdef INCA 79 ! Only INCA needs these informations (from the Earth's physics) 78 80 #include "indicesol.h" 79 81 #endif 80 82 INTEGER longcles 81 83 PARAMETER ( longcles = 20 ) … … 319 321 . ' restart ne correspondent pas a celles lues dans ' 320 322 write(lunout,*)' gcm.def' 321 322 323 323 write(lunout,*)' annee_ref=',annee_ref," anneeref=",anneeref 324 write(lunout,*)' day_ref=',day_ref," dayref=",dayref 325 write(lunout,*)' Pas de remise a zero' 324 326 ENDIF 325 327 326 c $$$if (annee_ref .ne. anneeref .or. day_ref .ne. dayref) then327 c $$$write(lunout,*)328 c $$$. 'GCM: Attention les dates initiales lues dans le fichier'329 c $$$write(lunout,*)330 c $$$. ' restart ne correspondent pas a celles lues dans '331 c $$$write(lunout,*)' gcm.def'332 c $$$write(lunout,*)' annee_ref=',annee_ref," anneeref=",anneeref333 c $$$write(lunout,*)' day_ref=',day_ref," dayref=",dayref334 c $$$if (raz_date .ne. 1) then335 c $$$write(lunout,*)336 c $$$. 'GCM: On garde les dates du fichier restart'337 c $$$else338 c $$$annee_ref = anneeref339 c $$$day_ref = dayref340 c $$$day_ini = dayref341 c $$$itau_dyn = 0342 c $$$itau_phy = 0343 c $$$time_0 = 0.344 c $$$write(lunout,*)345 c $$$. 'GCM: On reinitialise a la date lue dans gcm.def'346 c $$$endif347 c $$$ELSE348 c $$$raz_date = 0349 c $$$endif328 c if (annee_ref .ne. anneeref .or. day_ref .ne. dayref) then 329 c write(lunout,*) 330 c . 'GCM: Attention les dates initiales lues dans le fichier' 331 c write(lunout,*) 332 c . ' restart ne correspondent pas a celles lues dans ' 333 c write(lunout,*)' gcm.def' 334 c write(lunout,*)' annee_ref=',annee_ref," anneeref=",anneeref 335 c write(lunout,*)' day_ref=',day_ref," dayref=",dayref 336 c if (raz_date .ne. 1) then 337 c write(lunout,*) 338 c . 'GCM: On garde les dates du fichier restart' 339 c else 340 c annee_ref = anneeref 341 c day_ref = dayref 342 c day_ini = dayref 343 c itau_dyn = 0 344 c itau_phy = 0 345 c time_0 = 0. 346 c write(lunout,*) 347 c . 'GCM: On reinitialise a la date lue dans gcm.def' 348 c endif 349 c ELSE 350 c raz_date = 0 351 c endif 350 352 351 353 #ifdef CPP_IOIPSL … … 461 463 462 464 #ifdef CPP_IOIPSL 463 if ( 1.eq.1) then464 465 time_step = zdtvr 465 t_ops = iecri * daysec 466 t_wrt = iecri * daysec 467 ! CALL inithist(dynhist_file,day_ref,annee_ref,time_step, 468 ! . t_ops, t_wrt, histid, histvid) 469 470 ! IF (ok_dynzon) THEN 471 ! t_ops = iperiod * time_step 472 ! t_wrt = periodav * daysec 473 ! CALL initdynav(dynhistave_file,day_ref,annee_ref,time_step, 474 ! . t_ops, t_wrt, histaveid) 475 ! END IF 466 if (ok_dyn_ins) then 467 ! initialize output file for instantaneous outputs 468 ! t_ops = iecri * daysec ! do operations every t_ops 469 t_ops =((1.0*iecri)/day_step) * daysec 470 t_wrt = daysec ! iecri * daysec ! write output every t_wrt 471 CALL inithist(day_ref,annee_ref,time_step, 472 & t_ops,t_wrt) 473 endif 474 475 IF (ok_dyn_ave) THEN 476 ! initialize output file for averaged outputs 477 t_ops = iperiod * time_step ! do operations every t_ops 478 t_wrt = periodav * daysec ! write output every t_wrt 479 CALL initdynav(day_ref,annee_ref,time_step, 480 & t_ops,t_wrt) 481 END IF 476 482 dtav = iperiod*dtvr/daysec 477 endif478 479 480 483 #endif 481 484 ! #endif of #ifdef CPP_IOIPSL
Note: See TracChangeset
for help on using the changeset viewer.