Ignore:
Timestamp:
Dec 13, 2019, 11:41:41 AM (5 years ago)
Author:
slebonnois
Message:

SL: correction probleme sur le geopotentiel dans le cas Venus, ie Cp(T)... surtout problematique pour la haute atmosphere

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r1824 r2195  
    537537
    538538! Compute geopotential (physics might need it)
    539          CALL geopot  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )
     539! GEOP CORRECTION
     540! ADAPTATION GCM POUR CP(T)
     541!        CALL geopot  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )
     542         call tpot2t(ijp1llm,teta,temp,pk)
     543         tsurpk = cpp*temp/pk
     544         CALL geopot  ( ip1jmp1, tsurpk, pk, pks, phis, phi )
    540545
    541546           jD_cur = jD_ref + day_ini - day_ref +                        &
    … …  
    825830               
    826831!              ! Ehouarn: re-compute geopotential for outputs
    827                CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
     832! GEOP CORRECTION
     833! ADAPTATION GCM POUR CP(T)
     834!              CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
     835               call tpot2t(ijp1llm,teta,temp,pk)
     836               tsurpk = cpp*temp/pk
     837               CALL geopot(ip1jmp1,tsurpk,pk,pks,phis,phi)
    828838
    829839               IF (ok_dynzon) THEN
    … …  
    979989
    980990!              ! Ehouarn: re-compute geopotential for outputs
    981                CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
     991! GEOP CORRECTION
     992! ADAPTATION GCM POUR CP(T)
     993!              CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
     994               call tpot2t(ijp1llm,teta,temp,pk)
     995               tsurpk = cpp*temp/pk
     996               CALL geopot(ip1jmp1,tsurpk,pk,pks,phis,phi)
    982997
    983998               IF (ok_dynzon) THEN
Note: See TracChangeset for help on using the changeset viewer.