Changeset 4228
- Timestamp:
- Jul 26, 2022, 3:36:45 PM (2 years ago)
- Location:
- LMDZ6/trunk
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/DefLists/context_lmdz.xml
r4145 r4228 63 63 <axis id="time_year" unit="day" /> 64 64 <axis id="presnivs" standard_name="Vertical levels" unit="Pa"/> 65 <axis id="presinter" standard_name="Pressure at interface of vertical levels" unit="Pa"/> 65 66 <axis id="Ahyb" standard_name="Ahyb comp of Hyb Cord" unit="Pa"/> 66 67 <axis id="Bhyb" standard_name="Bhyb comp of Hyb Cord" unit=""/> … … 119 120 </grid> 120 121 122 <grid id="grid_glo_presinter"> 123 <domain domain_ref="dom_glo" /> 124 <axis axis_ref="presinter" /> 125 </grid> 126 <grid id="grid_out_presinter"> 127 <domain domain_ref="dom_out" /> 128 <axis axis_ref="presinter" /> 129 </grid> 130 121 131 <grid id="grid_out_plev17"> 122 132 <domain id="dom_out" /> -
LMDZ6/trunk/libf/dyn3d_common/comvert_mod.F90
r2602 r4228 10 10 11 11 PUBLIC :: ap,bp,presnivs,dpres,sig,ds,pa,preff,nivsigs,nivsig, & 12 aps,bps,scaleheight,pseudoalt,disvert_type, pressure_exner 12 aps,bps,scaleheight,pseudoalt,disvert_type, pressure_exner, & 13 presinter 13 14 14 15 REAL ap(llm+1) ! hybrid pressure contribution at interlayers 15 16 REAL bp (llm+1) ! hybrid sigma contribution at interlayer 16 17 REAL presnivs(llm) ! (reference) pressure at mid-layers 18 REAL presinter(llm+1) ! (reference) pressure at interlayers 17 19 REAL dpres(llm) 18 20 REAL sig(llm+1) -
LMDZ6/trunk/libf/dyn3d_common/disvert.F90
r2786 r4228 11 11 use assert_m, only: assert 12 12 USE comvert_mod, ONLY: ap, bp, aps, bps, nivsigs, nivsig, dpres, presnivs, & 13 pseudoalt, pa, preff, scaleheight 13 pseudoalt, pa, preff, scaleheight, presinter 14 14 USE logic_mod, ONLY: ok_strato 15 15 … … 35 35 ! dpres(llm) !--- PRESSURE DIFFERENCE FOR EACH LAYER 36 36 ! presnivs(llm) !--- PRESSURE AT EACH MID-LAYER 37 ! presinter(llm+1) !--- PRESSURE AT EACH INTERFACE 37 38 ! scaleheight !--- VERTICAL SCALE HEIGHT (Earth: 8kms) 38 39 ! nivsig(llm+1) !--- SIGMA INDEX OF EACH LAYER INTERFACE … … 355 356 max(ap(l+1)+bp(l+1)*preff, 1.e-10)) 356 357 ENDDO 358 DO l=1, llmp1 359 presinter(l)= ( ap(l)+bp(l)*preff) 360 write(lunout, *)'PRESINTER(', l, ')=', presinter(l) 361 ENDDO 357 362 358 363 write(lunout, *) trim(modname),': PRESNIVS ' -
LMDZ6/trunk/libf/dynphy_lonlat/inigeomphy_mod.F90
r3895 r4228 25 25 USE nrtype, ONLY: pi 26 26 USE comvert_mod, ONLY: preff, ap, bp, aps, bps, presnivs, & 27 scaleheight, pseudoalt 27 scaleheight, pseudoalt, presinter 28 28 USE vertical_layers_mod, ONLY: init_vertical_layers 29 29 IMPLICIT NONE … … 233 233 ! copy over preff , ap(), bp(), etc 234 234 CALL init_vertical_layers(nlayer,preff,scaleheight, & 235 ap,bp,aps,bps,presnivs,p seudoalt)235 ap,bp,aps,bps,presnivs,presinter,pseudoalt) 236 236 237 237 !$OMP END PARALLEL -
LMDZ6/trunk/libf/phy_common/vertical_layers_mod.F90
r2786 r4228 15 15 REAL,SAVE,ALLOCATABLE :: presnivs(:) ! reference pressure at mid-layer (Pa), 16 16 ! based on preff, ap and bp 17 REAL,SAVE,ALLOCATABLE :: presinter(:) ! reference pressure at interface (Pa), 18 ! based on preff, ap and bp 17 19 REAL,SAVE,ALLOCATABLE :: pseudoalt(:) ! pseudo-altitude of model layers (km), 18 20 ! based on preff and scaleheight 19 21 20 !$OMP THREADPRIVATE(preff,scaleheight,ap,bp,aps,bps,presnivs,p seudoalt)22 !$OMP THREADPRIVATE(preff,scaleheight,ap,bp,aps,bps,presnivs,presinter,pseudoalt) 21 23 22 24 … … 24 26 25 27 SUBROUTINE init_vertical_layers(nlayer,preff_,scaleheight_,ap_,bp_,& 26 aps_,bps_,presnivs_, p seudoalt_)28 aps_,bps_,presnivs_, presinter_, pseudoalt_) 27 29 IMPLICIT NONE 28 30 INTEGER,INTENT(IN) :: nlayer ! number of atmospheric layers … … 34 36 REAL,INTENT(IN) :: bps_(nlayer) ! hybrid coordinate at mid-layer 35 37 REAL,INTENT(IN) :: presnivs_(nlayer) ! Appproximative pressure of atm. layers (Pa) 38 REAL,INTENT(IN) :: presinter_(nlayer+1) ! Appproximative pressure of atm. layers (Pa) 36 39 REAL,INTENT(IN) :: pseudoalt_(nlayer) ! pseudo-altitude of atm. layers (km) 37 40 … … 41 44 ALLOCATE(bps(nlayer)) 42 45 ALLOCATE(presnivs(nlayer)) 46 ALLOCATE(presinter(nlayer+1)) 43 47 ALLOCATE(pseudoalt(nlayer)) 44 48 … … 50 54 bps(:) = bps_(:) 51 55 presnivs(:) = presnivs_(:) 56 presinter(:) = presinter_(:) 52 57 pseudoalt(:) = pseudoalt_(:) 53 58 -
LMDZ6/trunk/libf/phylmd/iophy.F90
-
Property
svn:keywords
changed from
Author Date Id Revision
toId
r4103 r4228 1249 1249 1250 1250 nlev=SIZE(field,2) 1251 IF (nlev.EQ.klev+1) THEN 1252 nlevx=klev 1253 ELSE 1254 nlevx=nlev 1255 ENDIF 1251 nlevx=nlev 1252 ! IF (nlev.EQ.klev+1) THEN 1253 ! nlevx=klev 1254 ! ELSE 1255 ! nlevx=nlev 1256 ! ENDIF 1256 1257 1257 1258 IF (SIZE(field,1) == klon) then -
Property
svn:keywords
changed from
-
LMDZ6/trunk/libf/phylmd/phys_output_mod.F90
-
Property
svn:keywords
changed from
Author Date Id Revision
toId
r4170 r4228 45 45 USE mod_grid_phy_lmdz, only: klon_glo,nbp_lon,nbp_lat 46 46 USE print_control_mod, ONLY: prt_level,lunout 47 USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps, pseudoalt 47 USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps, pseudoalt, presinter 48 48 USE time_phylmdz_mod, ONLY: day_ini, itau_phy, start_time, annee_ref, day_ref 49 49 #ifdef CPP_XIOS … … 347 347 CALL wxios_add_vaxis("presnivs", & 348 348 levmax(iff) - levmin(iff) + 1, presnivs(levmin(iff):levmax(iff))) 349 CALL wxios_add_vaxis("presinter", & 350 klev + 1, presinter(1:klev+1)) 349 351 CALL wxios_add_vaxis("Ahyb", & 350 352 levmax(iff) - levmin(iff) + 1, aps(levmin(iff):levmax(iff)), positif='down', & -
Property
svn:keywords
changed from
Note: See TracChangeset
for help on using the changeset viewer.