Changeset 4042 for LMDZ6/trunk/libf


Ignore:
Timestamp:
Dec 1, 2021, 9:15:05 PM (3 years ago)
Author:
lguez
Message:

Bug fix: define pk also for guide_plevs == 1

We will need pk later in guide_interp if guide_T and not guide_teta.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90

    r4041 r4042  
    10311031    endif ! of if (first)
    10321032   
     1033    CALL pression_loc( ijnb_u, ap, bp, psi, p )
     1034    if (disvert_type==1) then
     1035       CALL exner_hyb_loc(ijnb_u,psi,p,pks,pk)
     1036    else ! we assume that we are in the disvert_type==2 case
     1037       CALL exner_milieu_loc(ijnb_u,psi,p,pks,pk)
     1038    endif
     1039
    10331040! -----------------------------------------------------------------
    10341041! Calcul niveaux pression modele
     
    10461053        ENDDO
    10471054    ELSE
    1048         CALL pression_loc( ijnb_u, ap, bp, psi, p )
    1049         if (disvert_type==1) then
    1050           CALL exner_hyb_loc(ijnb_u,psi,p,pks,pk)
    1051         else ! we assume that we are in the disvert_type==2 case
    1052           CALL exner_milieu_loc(ijnb_u,psi,p,pks,pk)
    1053         endif
    10541055        unskap=1./kappa
    10551056!$OMP BARRIER
Note: See TracChangeset for help on using the changeset viewer.