Ignore:
Timestamp:
Oct 19, 2017, 12:18:04 PM (7 years ago)
Author:
Laurent Fairhead
Message:

Modifications for convergence of LMDZ6 and DYNAMICO physics
Tested with 1 year simulation of LMDZOR_v6 configuration

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/DYNAMICO-conv/libf/dyn3d/gcm.F90

    r2622 r3030  
    241241       'GCM: AVANT iniacademic AVANT AVANT AVANT AVANT'
    242242  if (.not.read_start) then
     243     annee_ref=anneeref
    243244     CALL iniacademic(vcov,ucov,teta,q,masse,ps,phis,time_0)
    244245  endif
     
    377378       tetagdiv, tetagrot , tetatemp, vert_prof_dissip)
    378379
     380  !  numero de stockage pour les fichiers de redemarrage:
     381
     382  !-----------------------------------------------------------------------
     383  !   Initialisation des I/O :
     384  !   ------------------------
     385
     386
     387  if (nday>=0) then
     388     day_end = day_ini + nday
     389  else
     390     day_end = day_ini - nday/day_step
     391  endif
     392  WRITE(lunout,300)day_ini,day_end
     393300 FORMAT('1'/,15x,'run du jour',i7,2x,'au jour',i7//)
     394
     395#ifdef CPP_IOIPSL
     396  call ju2ymds(jD_ref + day_ini - day_ref, an, mois, jour, heure)
     397  write (lunout,301)jour, mois, an
     398  call ju2ymds(jD_ref + day_end - day_ref, an, mois, jour, heure)
     399  write (lunout,302)jour, mois, an
     400301 FORMAT('1'/,15x,'run du ', i2,'/',i2,'/',i4)
     401302 FORMAT('1'/,15x,'    au ', i2,'/',i2,'/',i4)
     402#endif
     403
    379404  !-----------------------------------------------------------------------
    380405  !   Initialisation de la physique :
     
    391416#endif
    392417  ENDIF ! of IF ((iflag_phys==1).or.(iflag_phys>=100))
    393 
    394   !  numero de stockage pour les fichiers de redemarrage:
    395 
    396   !-----------------------------------------------------------------------
    397   !   Initialisation des I/O :
    398   !   ------------------------
    399 
    400 
    401   if (nday>=0) then
    402      day_end = day_ini + nday
    403   else
    404      day_end = day_ini - nday/day_step
    405   endif
    406   WRITE(lunout,300)day_ini,day_end
    407 300 FORMAT('1'/,15x,'run du jour',i7,2x,'au jour',i7//)
    408 
    409 #ifdef CPP_IOIPSL
    410   call ju2ymds(jD_ref + day_ini - day_ref, an, mois, jour, heure)
    411   write (lunout,301)jour, mois, an
    412   call ju2ymds(jD_ref + day_end - day_ref, an, mois, jour, heure)
    413   write (lunout,302)jour, mois, an
    414 301 FORMAT('1'/,15x,'run du ', i2,'/',i2,'/',i4)
    415 302 FORMAT('1'/,15x,'    au ', i2,'/',i2,'/',i4)
    416 #endif
    417418
    418419  !      if (planet_type.eq."earth") then
Note: See TracChangeset for help on using the changeset viewer.