Changeset 4042

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:

• LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90

@@ -1031 +1031 @@
     endif ! of if (first)
     
+    CALL pression_loc( ijnb_u, ap, bp, psi, p )
+    if (disvert_type==1) then
+       CALL exner_hyb_loc(ijnb_u,psi,p,pks,pk)
+    else ! we assume that we are in the disvert_type==2 case
+       CALL exner_milieu_loc(ijnb_u,psi,p,pks,pk)
+    endif
+
 ! -----------------------------------------------------------------
 ! Calcul niveaux pression modele 
@@ -1046 +1053 @@
         ENDDO
     ELSE
-        CALL pression_loc( ijnb_u, ap, bp, psi, p )
-        if (disvert_type==1) then
-          CALL exner_hyb_loc(ijnb_u,psi,p,pks,pk)
-        else ! we assume that we are in the disvert_type==2 case
-          CALL exner_milieu_loc(ijnb_u,psi,p,pks,pk)
-        endif
         unskap=1./kappa
 !$OMP BARRIER