Changeset 4379 for LMDZ6

Timestamp:

Jan 10, 2023, 8:14:07 PM (23 months ago)

Author:

fhourdin

Message:

preff : un pas en avant, un pas en arriere

Location:

LMDZ6/trunk/libf

Files:

• dyn3d/iniacademic.F90 (modified) (1 diff)
• dyn3d_common/conf_planete.F90 (modified) (1 diff)
• dyn3dmem/iniacademic_loc.F90 (modified) (1 diff)

LMDZ6/trunk/libf/dyn3d/iniacademic.F90

@@ -112 +112 @@
      ! specific initializations for the shallow water case
      kappa=1
-  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
 

LMDZ6/trunk/libf/dyn3d_common/conf_planete.F90

@@ -30 +30 @@
 
 !Reference surface pressure (Pa)
-preff=101325.
+! 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('preff', preff)
+
 ! Reference pressure at which hybrid coord. become purely pressure
 ! pa=50000.
 pa=preff/2.
 CALL getin('pa', pa)
+
 ! Gravity
 g=9.80665
+
 CALL getin('g',g)
 ! Molar mass of the atmosphere
+
 molmass = 28.9644
 CALL getin('molmass',molmass)
 ! kappa=R/Cp et Cp      
+
 kappa = 2./7.
 CALL getin('kappa',kappa)
+
 cpp=8.3145/molmass/kappa*1000.
 CALL getin('cpp',cpp)
 ! Radius of the planet
+
 rad = 6371229. 
 CALL getin('radius',rad)
+
 ! Length of a standard day (s)
 daysec=86400.
 CALL getin('daysec',daysec)
+
 ! Rotation rate of the planet:
 ! Length of a solar day, in standard days
 daylen = 1.
+
 CALL getin('daylen',daylen)
 ! Number of days (standard) per year:
+
 year_day = 365.25
 CALL getin('year_day',year_day)
 ! Omega
 ! omeg=2.*pi/86400.
+
 omeg=2.*pi/daysec*(1./daylen+1./year_day)
 CALL getin('omeg',omeg)

LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90

@@ -116 +116 @@
      ! specific initializations for the shallow water case
      kappa=1
-  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