Changeset 4376

Timestamp:

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

Author:

fhourdin

Message:

Imposing surface pressure from .def with iniacademic

Location:

LMDZ6/trunk/libf

Files:

• dyn3d/iniacademic.F90 (modified) (3 diffs)
• dyn3dmem/iniacademic_loc.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/dyn3d/iniacademic.F90

@@ -18 +18 @@
   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
@@ -114 +114 @@
   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
@@ -120 +130 @@
   ! Initialize pressure and mass field if read_start=.false.
   IF (.NOT. read_start) THEN
-     ! 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(:)=101080.  
-     else
-        ! use reference surface pressure
-        ps(:)=preff
-     endif
+     ps(:)=preff
 
      ! ground geopotential

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