Changeset 4054 for LMDZ6/trunk

Timestamp:

Jan 10, 2022, 1:52:20 PM (2 years ago)

Author:

lguez

Message:

Only compute pk if we are going to need it

In revsion revision 4042, wa added the computation of pk for
guide_plevs == 1. We can restrict a little better the cases where we
need pk.

File:

• LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90

@@ -1030 +1030 @@
 !$OMP END MASTER
     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
+
+    if (guide_plevs /= 1 .or. guide_t .and. .not. guide_teta &
+         .or. guide_q .and. guide_hr) then
+       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
+    end if
 
 ! -----------------------------------------------------------------