Changeset 1382 for LMDZ4/branches


Ignore:
Timestamp:
May 6, 2010, 3:44:09 PM (15 years ago)
Author:
yann meurdesoif
Message:

Correction sortie dynamque en parallèle

EM+YM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/gcm.F

    r1380 r1382  
    274274          CALL iniacademic(vcov,ucov,teta,q,masse,ps,phis,time_0)
    275275        endif
    276 
     276!#ifdef CPP_IOIPSL
    277277        if (planet_type.eq."earth") then
    278278#ifdef CPP_EARTH
    279279! Load an Earth-format start file
    280280         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
    290283        endif ! of if (planet_type.eq."earth")
    291 
    292284c       write(73,*) 'ucov',ucov
    293285c       write(74,*) 'vcov',vcov
     
    498490#ifdef CPP_IOIPSL
    499491      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)
    518509!         CALL initdynav_p(dynhistave_file,day_ref,annee_ref,time_step,
    519510!     .        t_ops, t_wrt, histaveid)
    520       END IF
     511        END IF
     512      ENDIF
    521513      dtav = iperiod*dtvr/daysec
    522514#endif
Note: See TracChangeset for help on using the changeset viewer.