Changeset 6189
- Timestamp:
- Apr 23, 2026, 5:22:06 PM (11 days ago)
- Location:
- LMDZ6/trunk/libf
- Files:
-
- 3 edited
-
phylmd/open_climoz_m.f90 (modified) (4 diffs)
-
phylmd/physiq_mod.F90 (modified) (2 diffs)
-
phylmdiso/physiq_mod.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/open_climoz_m.f90
r6122 r6189 9 9 !------------------------------------------------------------------------------- 10 10 ! 11 SUBROUTINE open_climoz(ncid, press_in_cen, press_in_edg,time_in, daily, adjust)11 SUBROUTINE open_climoz(ncid, press_in_cen, time_in, daily, adjust) 12 12 ! 13 13 !------------------------------------------------------------------------------- … … 28 28 INTEGER, INTENT(OUT):: ncid !--- "climoz_LMDZ.nc" identifier 29 29 REAL, allocatable, intent(out):: press_in_cen(:) !--- at cells centers 30 REAL, allocatable, INTENT(OUT):: press_in_edg(:) !--- at the interfaces (pressure intervals)31 30 REAL, allocatable, intent(out):: time_in(:) !--- records times, in days since Jan. 1st 32 31 LOGICAL, INTENT(IN) :: daily !--- daily files (calendar dependent days nb) … … 56 55 CALL NF95_GW_VAR(ncid, varid, time_in) 57 56 ntim = SIZE(time_in) 58 59 !--- BUILD EDGES OF PRESSURE INTERVALS: HALFWAY IN LOGARITHMS60 ALLOCATE(press_in_edg(nlev+1))61 press_in_edg=[0.,SQRT(press_in_cen(1:nlev-1)*press_in_cen(2:nlev)),HUGE(0.)]62 57 63 58 !--- CHECK RECORDS NUMBER AND DISPLAY CORRESPONDING INFORMATION … … 94 89 IF (.NOT. is_mpi_root .or. .not. is_omp_root) then 95 90 ALLOCATE(press_in_cen(nlev )) 96 ALLOCATE(press_in_edg(nlev+1))97 91 ALLOCATE(time_in(ntim)) 98 92 end IF 99 93 100 94 CALL bcast(press_in_cen) 101 CALL bcast(press_in_edg)102 95 CALL bcast(time_in) 103 96 -
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r6188 r6189 1060 1060 INTEGER,SAVE :: ncid_climoz ! NetCDF file identifier 1061 1061 REAL, ALLOCATABLE, SAVE :: press_cen_climoz(:) ! Pressure levels 1062 REAL, ALLOCATABLE, SAVE :: press_edg_climoz(:) ! Edges of pressure intervals1063 1062 REAL, ALLOCATABLE, SAVE :: time_climoz(:) ! Time vector 1064 1063 !$OMP THREADPRIVATE(ncid_climoz, press_cen_climoz) 1065 !$OMP THREADPRIVATE( press_edg_climoz,time_climoz)1064 !$OMP THREADPRIVATE(time_climoz) 1066 1065 1067 1066 CHARACTER(LEN=13), PARAMETER :: vars_climoz(2) & … … 2102 2101 2103 2102 IF (read_climoz >= 1) CALL open_climoz(ncid_climoz, press_cen_climoz, & 2104 press_edg_climoz,time_climoz, ok_daily_climoz, adjust_tropopause)2103 time_climoz, ok_daily_climoz, adjust_tropopause) 2105 2104 ! 2106 2105 !IM betaCRF -
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r6188 r6189 1135 1135 INTEGER,SAVE :: ncid_climoz ! NetCDF file identifier 1136 1136 REAL, ALLOCATABLE, SAVE :: press_cen_climoz(:) ! Pressure levels 1137 REAL, ALLOCATABLE, SAVE :: press_edg_climoz(:) ! Edges of pressure intervals1138 1137 REAL, ALLOCATABLE, SAVE :: time_climoz(:) ! Time vector 1139 1138 !$OMP THREADPRIVATE(ncid_climoz, press_cen_climoz) 1140 !$OMP THREADPRIVATE( press_edg_climoz,time_climoz)1139 !$OMP THREADPRIVATE(time_climoz) 1141 1140 CHARACTER(LEN=13), PARAMETER :: vars_climoz(2) & 1142 1141 = ["tro3 ","tro3_daylight"] … … 2241 2240 2242 2241 IF (read_climoz >= 1) CALL open_climoz(ncid_climoz, press_cen_climoz, & 2243 press_edg_climoz,time_climoz, ok_daily_climoz, adjust_tropopause)2242 time_climoz, ok_daily_climoz, adjust_tropopause) 2244 2243 ! 2245 2244 !IM betaCRF
Note: See TracChangeset
for help on using the changeset viewer.
