Ignore:
Timestamp:
Apr 24, 2017, 11:39:08 AM (8 years ago)
Author:
slebonnois
Message:

SL: petit correctif sur makelmdz

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F

    r1572 r1702  
    518518
    519519! Compute geopotential (physics might need it)
    520          CALL geopot  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )
     520         call tpot2t(ijp1llm,teta,temp,pk)
     521         tsurpk = cpp*temp/pk
     522         CALL geopot  ( ip1jmp1, tsurpk, pk, pks, phis, phi )
    521523
    522524           jD_cur = jD_ref + day_ini - day_ref +                        &
    … …  
    803805               
    804806!              ! Ehouarn: re-compute geopotential for outputs
    805                CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
     807               call tpot2t(ijp1llm,teta,temp,pk)
     808               tsurpk = cpp*temp/pk
     809               CALL geopot(ip1jmp1,tsurpk,pk,pks,phis,phi)
    806810
    807811               IF (ok_dynzon) THEN
    … …  
    947951
    948952!              ! Ehouarn: re-compute geopotential for outputs
    949                CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
     953! ADAPTATION GCM POUR CP(T)
     954               call tpot2t(ijp1llm,teta,temp,pk)
     955               tsurpk = cpp*temp/pk
     956               CALL geopot(ip1jmp1,tsurpk,pk,pks,phis,phi)
    950957
    951958               IF (ok_dynzon) THEN
Note: See TracChangeset for help on using the changeset viewer.