Changeset 3991 for trunk/LMDZ.COMMON/libf/evolution/pemredem.F90
- Timestamp:
- Dec 16, 2025, 4:39:24 PM (3 months ago)
- File:
-
- 1 edited
-
trunk/LMDZ.COMMON/libf/evolution/pemredem.F90 (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/evolution/pemredem.F90
r3989 r3991 1 1 MODULE pemredem 2 !----------------------------------------------------------------------- 3 ! NAME 4 ! pemredem 5 ! 6 ! DESCRIPTION 7 ! Write PEM-specific NetCDF restart files. 8 ! 9 ! AUTHORS & DATE 10 ! L. Lange 11 ! JB Clement, 2023-2025 12 ! 13 ! NOTES 14 ! Inspired by phyredem from the PCM. Handles time-independent and 15 ! time-dependent variables for restart functionality. 16 !----------------------------------------------------------------------- 2 17 3 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 4 !!! 5 !!! Purpose: Write specific netcdf restart for the PEM 6 !!! 7 !!! 8 !!! Author: LL, inspired by phyredem from the PCM 9 !!! 10 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 11 18 ! DECLARATION 19 ! ----------- 12 20 implicit none 13 21 22 contains 23 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 24 14 25 !======================================================================= 15 contains 16 !======================================================================= 26 SUBROUTINE pemdem0(filename,lonfi,latfi,cell_area,ngrid,nslope,def_slope,subslope_dist) 27 !----------------------------------------------------------------------- 28 ! NAME 29 ! pemdem0 30 ! 31 ! DESCRIPTION 32 ! Create physics restart file and write time-independent variables. 33 ! 34 ! AUTHORS & DATE 35 ! L. Lange 36 ! JB Clement, 2023-2025 37 ! 38 ! NOTES 39 ! 40 !----------------------------------------------------------------------- 17 41 18 SUBROUTINE pemdem0(filename,lonfi,latfi,cell_area,ngrid,nslope,def_slope,subslope_dist) 19 20 ! create physics restart file and write time-independent variables 21 use soil, only: mlayer_PEM 42 ! DEPENDENCIES 43 ! ------------ 44 use soil, only: mlayer_PEM 22 45 use iostart_pem, only: open_restartphy, close_restartphy, put_var, put_field, length 23 46 47 ! DECLARATION 48 ! ----------- 24 49 implicit none 25 50 51 ! ARGUMENTS 52 ! --------- 26 53 character(*), intent(in) :: filename 27 54 integer, intent(in) :: ngrid, nslope … … 31 58 real, dimension(ngrid,nslope), intent(in) :: subslope_dist ! undermesh statistics 32 59 33 ! Create physics start file 60 ! CODE 61 ! ---- 34 62 call open_restartphy(filename) 35 63 … … 54 82 55 83 END SUBROUTINE pemdem0 84 !======================================================================= 56 85 57 86 !======================================================================= 58 59 87 SUBROUTINE pemdem1(filename,i_myear,nsoil_PEM,ngrid,nslope,tsoil_slope_PEM,inertiesoil_slope_PEM, & 60 88 icetable_depth,icetable_thickness,ice_porefilling,m_co2_regolith,m_h2o_regolith,h2o_ice,co2_ice,layerings_map) 89 !----------------------------------------------------------------------- 90 ! NAME 91 ! pemdem1 92 ! 93 ! DESCRIPTION 94 ! Write time-dependent variables to restart file. 95 ! 96 ! AUTHORS & DATE 97 ! L. Lange 98 ! JB Clement, 2023-2025 99 ! 100 ! NOTES 101 ! 102 !----------------------------------------------------------------------- 61 103 62 ! write time-dependent variable to restart file 104 ! DEPENDENCIES 105 ! ------------ 63 106 use iostart_pem, only: open_restartphy, close_restartphy, put_var, put_field 64 107 use soil, only: inertiedat_PEM, do_soil … … 66 109 use layered_deposits, only: layering, nb_str_max, map2array, print_layering, layering_algo 67 110 111 ! DECLARATION 112 ! ----------- 68 113 implicit none 69 114 115 ! ARGUMENTS 116 ! --------- 70 117 character(*), intent(in) :: filename 71 118 integer, intent(in) :: nsoil_PEM, ngrid, nslope … … 80 127 type(layering), dimension(ngrid,nslope), intent(in) :: layerings_map ! Layerings 81 128 129 ! LOCAL VARIABLES 130 ! --------------- 82 131 integer :: islope 83 132 character(2) :: num … … 85 134 real, dimension(:,:,:,:), allocatable :: layerings_array ! Array for stratification (layerings) 86 135 136 ! CODE 137 ! ---- 87 138 ! Open file 88 139 call open_restartphy(filename)
Note: See TracChangeset
for help on using the changeset viewer.
