- Timestamp:
- Jan 8, 2023, 2:30:23 AM (21 months ago)
- Location:
- LMDZ6/trunk/libf
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d/iniacademic.F90
r4325 r4376 18 18 USE comconst_mod, ONLY: cpp, kappa, g, daysec, dtvr, pi, im, jm 19 19 USE logic_mod, ONLY: iflag_phys, read_start 20 USE comvert_mod, ONLY: ap, bp, preff, p resnivs, pressure_exner20 USE comvert_mod, ONLY: ap, bp, preff, pa, presnivs, pressure_exner 21 21 USE temps_mod, ONLY: annee_ref, day_ini, day_ref 22 22 USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0 … … 114 114 endif 115 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. 124 endif 125 116 126 CALL iniconst 117 127 CALL inigeom … … 120 130 ! Initialize pressure and mass field if read_start=.false. 121 131 IF (.NOT. read_start) THEN 122 ! surface pressure 123 if (iflag_phys>2) then 124 ! specific value for CMIP5 aqua/terra planets 125 ! "Specify the initial dry mass to be equivalent to 126 ! a global mean surface pressure (101325 minus 245) Pa." 127 ps(:)=101080. 128 else 129 ! use reference surface pressure 130 ps(:)=preff 131 endif 132 ps(:)=preff 132 133 133 134 ! ground geopotential -
LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90
r4325 r4376 19 19 USE comconst_mod, ONLY: cpp, kappa, g, daysec, dtvr, pi, im, jm 20 20 USE logic_mod, ONLY: iflag_phys, read_start 21 USE comvert_mod, ONLY: ap, bp, preff, p resnivs, pressure_exner21 USE comvert_mod, ONLY: ap, bp, preff, pa, presnivs, pressure_exner 22 22 USE temps_mod, ONLY: annee_ref, day_ini, day_ref 23 23 USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0 … … 118 118 endif 119 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. 127 endif 128 120 129 CALL iniconst 121 130 CALL inigeom … … 133 142 134 143 ! surface pressure 135 if (iflag_phys>2) then 136 ! specific value for CMIP5 aqua/terra planets 137 ! "Specify the initial dry mass to be equivalent to 138 ! a global mean surface pressure (101325 minus 245) Pa." 139 ps_glo(:)=101080. 140 else 141 ! use reference surface pressure 142 ps_glo(:)=preff 143 endif 144 ps_glo(:)=preff 144 145 145 146 ! ground geopotential
Note: See TracChangeset
for help on using the changeset viewer.