Changeset 1703 for trunk


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é

Location:
trunk/LMDZ.COMMON/libf
Files:
3 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

  • trunk/LMDZ.COMMON/libf/dyn3d_common/geopot.F

    r1702 r1703  
    2121c     .. ngrid,teta,pk,pks,phis sont des argum. d'entree pour le s-pg ..
    2222c              phi               est un  argum. de sortie pour le s-pg .
    23 c
    24 c  This computation (with teta = cp T / pk !) is identical to
    25 c     delta phi = R/RMD T/p delta p         (r=R/RMD=cpp*kappa)
    2623c
    2724c=======================================================================

  • trunk/LMDZ.COMMON/libf/dyn3dpar/leapfrog_p.F

    r1702 r1703  
    862862c$OMP BARRIER
    863863! Compute geopotential (physics might need it)
    864 
    865          call tpot2t_glo_p(teta,temp,pk)
    866          ijb=ij_begin
    867          ije=ij_end
    868 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    869          do l=1,llm
    870            tsurpk(ijb:ije,l)=cpp*temp(ijb:ije,l)/pk(ijb:ije,l)
    871          enddo
    872 !$OMP END DO
    873 c$OMP MASTER
    874 
    875 !        CALL geopot_p(ip1jmp1,teta,pk,pks,phis,phi)
    876          CALL geopot_p( ip1jmp1, tsurpk, pk, pks, phis, phi )
    877 c$OMP END MASTER
     864         CALL geopot_p  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )
    878865
    879866           jD_cur = jD_ref + day_ini - day_ref
Note: See TracChangeset for help on using the changeset viewer.