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

SL: rectification du precedent commit pour les routines qui ne devaient pas etre concernées..... Désolé

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r1702 r1703  
    518518
    519519! Compute geopotential (physics might need it)
    520          call tpot2t(ijp1llm,teta,temp,pk)
    521          tsurpk = cpp*temp/pk
    522          CALL geopot  ( ip1jmp1, tsurpk, pk, pks, phis, phi )
     520         CALL geopot  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )
    523521
    524522           jD_cur = jD_ref + day_ini - day_ref +                        &
    … …  
    805803               
    806804!              ! Ehouarn: re-compute geopotential for outputs
    807                call tpot2t(ijp1llm,teta,temp,pk)
    808                tsurpk = cpp*temp/pk
    809                CALL geopot(ip1jmp1,tsurpk,pk,pks,phis,phi)
     805               CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
    810806
    811807               IF (ok_dynzon) THEN
    … …  
    951947
    952948!              ! Ehouarn: re-compute geopotential for outputs
    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)
     949               CALL geopot(ip1jmp1,teta,pk,pks,phis,phi)
    957950
    958951               IF (ok_dynzon) THEN
Note: See TracChangeset for help on using the changeset viewer.