Changeset 4379


Ignore:
Timestamp:
Jan 10, 2023, 8:14:07 PM (16 months ago)
Author:
fhourdin
Message:

preff : un pas en avant, un pas en arriere

Location:
LMDZ6/trunk/libf
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/dyn3d/iniacademic.F90

    r4376 r4379  
    112112     ! specific initializations for the shallow water case
    113113     kappa=1
    114   endif
    115 
    116 
    117   if ( ( .not. read_start ) .and. ( iflag_phys > 2 ) ) then
    118         ! 101080 : specific value for CMIP5 aqua/terra planets
    119         ! "Specify the initial dry mass to be equivalent to
    120         !  a global mean surface pressure (101325 minus 245) Pa."
    121         preff=101080.
    122         call getin('psurf',preff)
    123         pa=preff/2.
    124114  endif
    125115

  • LMDZ6/trunk/libf/dyn3d_common/conf_planete.F90

    r2600 r4379  
    3030
    3131!Reference surface pressure (Pa)
    32 preff=101325.
     32! 101080 : specific value for CMIP5 aqua/terra planets
     33! "Specify the initial dry mass to be equivalent to
     34!  a global mean surface pressure (101325 minus 245) Pa."
     35preff=101080.
    3336CALL getin('preff', preff)
     37
    3438! Reference pressure at which hybrid coord. become purely pressure
    3539! pa=50000.
    3640pa=preff/2.
    3741CALL getin('pa', pa)
     42
    3843! Gravity
    3944g=9.80665
     45
    4046CALL getin('g',g)
    4147! Molar mass of the atmosphere
     48
    4249molmass = 28.9644
    4350CALL getin('molmass',molmass)
    4451! kappa=R/Cp et Cp     
     52
    4553kappa = 2./7.
    4654CALL getin('kappa',kappa)
     55
    4756cpp=8.3145/molmass/kappa*1000.
    4857CALL getin('cpp',cpp)
    4958! Radius of the planet
     59
    5060rad = 6371229.
    5161CALL getin('radius',rad)
     62
    5263! Length of a standard day (s)
    5364daysec=86400.
    5465CALL getin('daysec',daysec)
     66
    5567! Rotation rate of the planet:
    5668! Length of a solar day, in standard days
    5769daylen = 1.
     70
    5871CALL getin('daylen',daylen)
    5972! Number of days (standard) per year:
     73
    6074year_day = 365.25
    6175CALL getin('year_day',year_day)
    6276! Omega
    6377! omeg=2.*pi/86400.
     78
    6479omeg=2.*pi/daysec*(1./daylen+1./year_day)
    6580CALL getin('omeg',omeg)

  • LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90

    r4376 r4379  
    116116     ! specific initializations for the shallow water case
    117117     kappa=1
    118   endif
    119 
    120   if ( ( .not. read_start ) .and. (iflag_phys>2) ) then
    121      ! 101080 : specific value for CMIP5 aqua/terra planets
    122      ! "Specify the initial dry mass to be equivalent to
    123      !  a global mean surface pressure (101325 minus 245) Pa."
    124      preff=101080.
    125      call getin('psurf',preff)
    126      pa=preff/2.
    127118  endif
    128119
Note: See TracChangeset for help on using the changeset viewer.