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/dyn3dpar/leapfrog_p.F

    r1564 r1702  
    862862c$OMP BARRIER
    863863! Compute geopotential (physics might need it)
    864          CALL geopot_p  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )
     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
     873c$OMP MASTER
     874
     875!        CALL geopot_p(ip1jmp1,teta,pk,pks,phis,phi)
     876         CALL geopot_p( ip1jmp1, tsurpk, pk, pks, phis, phi )
     877c$OMP END MASTER
    865878
    866879           jD_cur = jD_ref + day_ini - day_ref
Note: See TracChangeset for help on using the changeset viewer.