Ignore:
Timestamp:
Oct 22, 2010, 11:27:25 AM (14 years ago)
Author:
Ehouarn Millour
Message:

Implemented modifications to enable running with only one tracer for planet types different from "earth". Rem: If flag 'planet_type' is set to "earth" (default behaviour) then there must be at least 2 tracers for the dynamics to function properly.

These updates do not induce any changes in model outputs with respect to previous revisions.

EM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/branches/LMDZ5V1.0-dev/libf/dyn3dpar/leapfrog_p.F

    r1438 r1446  
    234234
    235235c$OMP MASTER
    236       dq=0.
     236      dq(:,:,:)=0.
    237237      CALL pression ( ip1jmp1, ap, bp, ps, p       )
    238238      CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
     
    14881488c$OMP MASTER
    14891489
    1490               if (planet_type.eq."earth") then
     1490!              if (planet_type.eq."earth") then
    14911491! Write an Earth-format restart file
    14921492                CALL dynredem1_p("restart.nc",0.0,
    14931493     &                           vcov,ucov,teta,q,masse,ps)
    1494               endif ! of if (planet_type.eq."earth")
     1494!              endif ! of if (planet_type.eq."earth")
    14951495
    14961496!              CLOSE(99)
     
    16811681
    16821682              IF(itau.EQ.itaufin) THEN
    1683                 if (planet_type.eq."earth") then
     1683!                if (planet_type.eq."earth") then
    16841684c$OMP MASTER
    16851685                   CALL dynredem1_p("restart.nc",0.0,
    16861686     .                               vcov,ucov,teta,q,masse,ps)
    16871687c$OMP END MASTER
    1688                 endif ! of if (planet_type.eq."earth")
     1688!                endif ! of if (planet_type.eq."earth")
    16891689              ENDIF ! of IF(itau.EQ.itaufin)
    16901690
Note: See TracChangeset for help on using the changeset viewer.