Changeset 4376 for LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90

Timestamp:

Jan 8, 2023, 2:30:23 AM (17 months ago)

Author:

fhourdin

Message:

Imposing surface pressure from .def with iniacademic

File:

• LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90

@@ -19 +19 @@
   USE comconst_mod, ONLY: cpp, kappa, g, daysec, dtvr, pi, im, jm
   USE logic_mod, ONLY: iflag_phys, read_start
-  USE comvert_mod, ONLY: ap, bp, preff, presnivs, pressure_exner
+  USE comvert_mod, ONLY: ap, bp, preff, pa, presnivs, pressure_exner
   USE temps_mod, ONLY: annee_ref, day_ini, day_ref
   USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0
@@ -118 +118 @@
   endif
 
+  if ( ( .not. read_start ) .and. (iflag_phys>2) ) then
+     ! 101080 : specific value for CMIP5 aqua/terra planets
+     ! "Specify the initial dry mass to be equivalent to
+     !  a global mean surface pressure (101325 minus 245) Pa."
+     preff=101080.
+     call getin('psurf',preff)
+     pa=preff/2.
+  endif
+
   CALL iniconst
   CALL inigeom
@@ -133 +142 @@
 
      ! surface pressure
-     if (iflag_phys>2) then
-        ! specific value for CMIP5 aqua/terra planets
-        ! "Specify the initial dry mass to be equivalent to
-        !  a global mean surface pressure (101325 minus 245) Pa."
-        ps_glo(:)=101080.  
-     else
-        ! use reference surface pressure
-        ps_glo(:)=preff
-     endif
+     ps_glo(:)=preff
 
      ! ground geopotential