| 1 | module phyredem |
|---|
| 2 | |
|---|
| 3 | implicit none |
|---|
| 4 | |
|---|
| 5 | contains |
|---|
| 6 | |
|---|
| 7 | subroutine physdem0(filename,lonfi,latfi,nsoil,ngrid,nlay,nq, & |
|---|
| 8 | phystep,day_ini,time,airefi, & |
|---|
| 9 | alb,ith,pzmea,pzstd,pzsig,pzgam,pzthe) |
|---|
| 10 | ! create physics restart file and write time-independent variables |
|---|
| 11 | use comsoil_h, only: volcapa, mlayer |
|---|
| 12 | use geometry_mod, only: cell_area |
|---|
| 13 | use surfdat_h, only: zmea, zstd, zsig, zgam, zthe, & |
|---|
| 14 | emisice, emissiv, & |
|---|
| 15 | iceradius, dtemisice, phisfi |
|---|
| 16 | use iostart, only : create_restartphy, close_restartphy, & |
|---|
| 17 | put_var, put_field, length |
|---|
| 18 | use mod_grid_phy_lmdz, only : klon_glo |
|---|
| 19 | use planete_mod, only: year_day, periastr, apoastr, peri_day, & |
|---|
| 20 | obliquit, z0, lmixmin, emin_turb |
|---|
| 21 | use comcstfi_mod, only: rad, omeg, g, mugaz, rcp |
|---|
| 22 | use time_phylmdz_mod, only: daysec |
|---|
| 23 | |
|---|
| 24 | implicit none |
|---|
| 25 | |
|---|
| 26 | character(len=*), intent(in) :: filename |
|---|
| 27 | real,intent(in) :: lonfi(ngrid) |
|---|
| 28 | real,intent(in) :: latfi(ngrid) |
|---|
| 29 | integer,intent(in) :: nsoil |
|---|
| 30 | integer,intent(in) :: ngrid |
|---|
| 31 | integer,intent(in) :: nlay |
|---|
| 32 | integer,intent(in) :: nq |
|---|
| 33 | real,intent(in) :: phystep |
|---|
| 34 | real,intent(in) :: day_ini |
|---|
| 35 | real,intent(in) :: time |
|---|
| 36 | real,intent(in) :: airefi(ngrid) |
|---|
| 37 | real,intent(in) :: alb(ngrid) |
|---|
| 38 | real,intent(in) :: ith(ngrid,nsoil) |
|---|
| 39 | real,intent(in) :: pzmea(ngrid) |
|---|
| 40 | real,intent(in) :: pzstd(ngrid) |
|---|
| 41 | real,intent(in) :: pzsig(ngrid) |
|---|
| 42 | real,intent(in) :: pzgam(ngrid) |
|---|
| 43 | real,intent(in) :: pzthe(ngrid) |
|---|
| 44 | |
|---|
| 45 | real :: tab_cntrl(length) ! nb "length=100" defined in iostart module |
|---|
| 46 | |
|---|
| 47 | ! Create physics start file |
|---|
| 48 | call create_restartphy(filename) |
|---|
| 49 | |
|---|
| 50 | ! tab_cntrl() contains run parameters |
|---|
| 51 | tab_cntrl(:)=0 ! initialization |
|---|
| 52 | !cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
|---|
| 53 | ! Fill control array tab_cntrl(:) with paramleters for this run |
|---|
| 54 | !cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
|---|
| 55 | ! Informations on the physics grid |
|---|
| 56 | tab_cntrl(1) = float(klon_glo) ! number of nodes on physics grid |
|---|
| 57 | tab_cntrl(2) = float(nlay) ! number of atmospheric layers |
|---|
| 58 | tab_cntrl(3) = day_ini + int(time) ! final day |
|---|
| 59 | tab_cntrl(4) = time -int(time) ! final time of day |
|---|
| 60 | |
|---|
| 61 | ! Informations about Mars, used by dynamics and physics |
|---|
| 62 | tab_cntrl(5) = rad ! radius of Mars (m) ~3397200 |
|---|
| 63 | tab_cntrl(6) = omeg ! rotation rate (rad.s-1) |
|---|
| 64 | tab_cntrl(7) = g ! gravity (m.s-2) ~3.72 |
|---|
| 65 | tab_cntrl(8) = mugaz ! Molar mass of the atmosphere (g.mol-1) ~43.49 |
|---|
| 66 | tab_cntrl(9) = rcp ! = r/cp ~0.256793 (=kappa dans dynamique) |
|---|
| 67 | tab_cntrl(10) = daysec ! length of a sol (s) ~88775 |
|---|
| 68 | |
|---|
| 69 | tab_cntrl(11) = phystep ! time step in the physics |
|---|
| 70 | tab_cntrl(12) = 0. |
|---|
| 71 | tab_cntrl(13) = 0. |
|---|
| 72 | |
|---|
| 73 | ! Informations about Mars, only for physics |
|---|
| 74 | tab_cntrl(14) = year_day ! length of year (sols) ~668.6 |
|---|
| 75 | tab_cntrl(15) = periastr ! min. star-planet distance (AU) |
|---|
| 76 | tab_cntrl(16) = apoastr ! max. star-planet distance (AU) |
|---|
| 77 | tab_cntrl(17) = peri_day ! date of periastron (sols since N. spring) |
|---|
| 78 | tab_cntrl(18) = obliquit ! Obliquity of the planet (deg) ~23.98 |
|---|
| 79 | |
|---|
| 80 | ! Boundary layer and turbulence |
|---|
| 81 | tab_cntrl(19) = z0 ! surface roughness (m) ~0.01 |
|---|
| 82 | tab_cntrl(20) = lmixmin ! mixing length ~100 |
|---|
| 83 | tab_cntrl(21) = emin_turb ! minimal energy ~1.e-8 |
|---|
| 84 | |
|---|
| 85 | ! Optical properties of polar caps and ground emissivity |
|---|
| 86 | tab_cntrl(24) = emisice(1) ! Emissivity of northern cap ~0.95 |
|---|
| 87 | tab_cntrl(25) = emisice(2) ! Emissivity of southern cap ~0.95 |
|---|
| 88 | tab_cntrl(26) = emissiv ! Emissivity of martian soil ~.95 |
|---|
| 89 | tab_cntrl(31) = iceradius(1) ! mean scat radius of CO2 snow (north) |
|---|
| 90 | tab_cntrl(32) = iceradius(2) ! mean scat radius of CO2 snow (south) |
|---|
| 91 | tab_cntrl(33) = dtemisice(1) ! time scale for snow metamorphism (north) |
|---|
| 92 | tab_cntrl(34) = dtemisice(2) ! time scale for snow metamorphism (south) |
|---|
| 93 | |
|---|
| 94 | tab_cntrl(28) = 0. |
|---|
| 95 | tab_cntrl(29) = 0. |
|---|
| 96 | tab_cntrl(30) = 0. |
|---|
| 97 | ! Soil properties: |
|---|
| 98 | tab_cntrl(35) = volcapa ! soil volumetric heat capacity |
|---|
| 99 | |
|---|
| 100 | call put_var("controle","Control parameters",tab_cntrl) |
|---|
| 101 | |
|---|
| 102 | ! Write the mid-layer depths |
|---|
| 103 | call put_var("soildepth","Soil mid-layer depth",mlayer) |
|---|
| 104 | |
|---|
| 105 | ! Write longitudes |
|---|
| 106 | call put_field("longitude","Longitudes of physics grid",lonfi) |
|---|
| 107 | |
|---|
| 108 | ! Write latitudes |
|---|
| 109 | call put_field("latitude","Latitudes of physics grid",latfi) |
|---|
| 110 | |
|---|
| 111 | ! Write mesh areas |
|---|
| 112 | call put_field("area","Mesh area",cell_area) |
|---|
| 113 | |
|---|
| 114 | ! Write surface geopotential |
|---|
| 115 | call put_field("phisfi","Geopotential at the surface",phisfi) |
|---|
| 116 | |
|---|
| 117 | ! Write surface albedo |
|---|
| 118 | call put_field("albedodat","Albedo of bare ground",alb) |
|---|
| 119 | |
|---|
| 120 | ! Subgrid topogaphy variables |
|---|
| 121 | call put_field("ZMEA","Relief: mean relief",zmea) |
|---|
| 122 | call put_field("ZSTD","Relief: standard deviation",zstd) |
|---|
| 123 | call put_field("ZSIG","Relief: sigma parameter",zsig) |
|---|
| 124 | call put_field("ZGAM","Relief: gamma parameter",zgam) |
|---|
| 125 | call put_field("ZTHE","Relief: theta parameter",zthe) |
|---|
| 126 | |
|---|
| 127 | ! Soil thermal inertia |
|---|
| 128 | call put_field("inertiedat","Soil thermal inertia",ith) |
|---|
| 129 | |
|---|
| 130 | ! Close file |
|---|
| 131 | call close_restartphy |
|---|
| 132 | |
|---|
| 133 | end subroutine physdem0 |
|---|
| 134 | |
|---|
| 135 | subroutine physdem1(filename,nsoil,ngrid,nlay,nq, & |
|---|
| 136 | phystep,time,tsurf,tsoil,emis,albedo,q2,qsurf, & |
|---|
| 137 | cloudfrac,totcloudfrac,hice, & |
|---|
| 138 | rnat,pctsrf_sic,tslab,tsea_ice,tice,sea_ice) |
|---|
| 139 | ! write time-dependent variable to restart file |
|---|
| 140 | use iostart, only : open_restartphy, close_restartphy, & |
|---|
| 141 | put_var, put_field |
|---|
| 142 | use tracer_h, only: noms |
|---|
| 143 | USE radinc_h, ONLY: L_NSPECTV |
|---|
| 144 | ! use slab_ice_h, only: noceanmx |
|---|
| 145 | USE ocean_slab_mod, ONLY: nslay |
|---|
| 146 | use callkeys_mod, only: ok_slab_ocean, calllott_nonoro |
|---|
| 147 | use nonoro_gwd_ran_mod, only: du_nonoro_gwd, dv_nonoro_gwd, & |
|---|
| 148 | east_gwstress, west_gwstress |
|---|
| 149 | |
|---|
| 150 | implicit none |
|---|
| 151 | |
|---|
| 152 | character(len=*),intent(in) :: filename |
|---|
| 153 | integer,intent(in) :: nsoil |
|---|
| 154 | integer,intent(in) :: ngrid |
|---|
| 155 | integer,intent(in) :: nlay |
|---|
| 156 | integer,intent(in) :: nq |
|---|
| 157 | real,intent(in) :: phystep |
|---|
| 158 | real,intent(in) :: time |
|---|
| 159 | real,intent(in) :: tsurf(ngrid) |
|---|
| 160 | real,intent(in) :: tsoil(ngrid,nsoil) |
|---|
| 161 | real,intent(in) :: emis(ngrid) |
|---|
| 162 | real,intent(in) :: albedo(ngrid,L_NSPECTV) |
|---|
| 163 | real,intent(in) :: q2(ngrid,nlay+1) |
|---|
| 164 | real,intent(in) :: qsurf(ngrid,nq) |
|---|
| 165 | real,intent(in) :: cloudfrac(ngrid,nlay) |
|---|
| 166 | real,intent(in) :: totcloudfrac(ngrid) |
|---|
| 167 | real,intent(in) :: hice(ngrid) |
|---|
| 168 | real,intent(in) :: rnat(ngrid) |
|---|
| 169 | real,intent(in) :: pctsrf_sic(ngrid) |
|---|
| 170 | real,intent(in) :: tslab(ngrid,nslay) |
|---|
| 171 | real,intent(in) :: tsea_ice(ngrid) |
|---|
| 172 | real,intent(in) :: sea_ice(ngrid) |
|---|
| 173 | real,intent(in) :: tice(ngrid) |
|---|
| 174 | |
|---|
| 175 | integer :: iq |
|---|
| 176 | |
|---|
| 177 | ! Open file |
|---|
| 178 | call open_restartphy(filename) |
|---|
| 179 | |
|---|
| 180 | ! First variable to write must be "Time", in order to correctly |
|---|
| 181 | ! set time counter in file |
|---|
| 182 | !call put_var("Time","Temps de simulation",time) |
|---|
| 183 | |
|---|
| 184 | ! Surface temperature |
|---|
| 185 | call put_field("tsurf","Surface temperature",tsurf) |
|---|
| 186 | |
|---|
| 187 | ! Soil temperature |
|---|
| 188 | call put_field("tsoil","Soil temperature",tsoil) |
|---|
| 189 | |
|---|
| 190 | ! Emissivity of the surface |
|---|
| 191 | call put_field("emis","Surface emissivity",emis) |
|---|
| 192 | |
|---|
| 193 | ! Surface albedo (assume homegeneous spectral albedo for now) |
|---|
| 194 | call put_field("albedo","Surface albedo",albedo(:,1)) |
|---|
| 195 | |
|---|
| 196 | ! Planetary Boundary Layer |
|---|
| 197 | call put_field("q2","pbl wind variance",q2) |
|---|
| 198 | |
|---|
| 199 | ! cloud fraction and sea ice (NB: these should be optional... to be improved) |
|---|
| 200 | call put_field("cloudfrac","Cloud fraction",cloudfrac) |
|---|
| 201 | call put_field("totcloudfrac","Total fraction",totcloudfrac) |
|---|
| 202 | call put_field("hice","Height of oceanic ice",hice) |
|---|
| 203 | |
|---|
| 204 | !Slab ocean |
|---|
| 205 | if(ok_slab_ocean) then |
|---|
| 206 | call put_field("rnat","Nature of surface",rnat) |
|---|
| 207 | call put_field("pctsrf_sic","Pourcentage sea ice",pctsrf_sic) |
|---|
| 208 | call put_field("tslab","Temperature slab ocean",tslab) |
|---|
| 209 | call put_field("tsea_ice","Temperature of top layer (seaice or snow)",tsea_ice) |
|---|
| 210 | call put_field("tice","Temperature of seaice",tice) |
|---|
| 211 | call put_field("sea_ice","Sea ice mass",sea_ice) |
|---|
| 212 | endif!(ok_slab_ocean) |
|---|
| 213 | |
|---|
| 214 | ! tracers |
|---|
| 215 | if (nq>0) then |
|---|
| 216 | do iq=1,nq |
|---|
| 217 | call put_field(noms(iq),"tracer on surface",qsurf(:,iq)) |
|---|
| 218 | enddo |
|---|
| 219 | endif ! of if (nq>0) |
|---|
| 220 | |
|---|
| 221 | ! Non-orographic gavity waves |
|---|
| 222 | if (calllott_nonoro) then |
|---|
| 223 | call put_field("du_nonoro_gwd","Zonal wind tendency due to GW",du_nonoro_gwd) |
|---|
| 224 | call put_field("dv_nonoro_gwd","Meridional wind tendency due to GW",dv_nonoro_gwd) |
|---|
| 225 | call put_field("east_gwstress","Eastward stress profile due to GW",east_gwstress) |
|---|
| 226 | call put_field("west_gwstress","Westward stress profile due to GW",west_gwstress) |
|---|
| 227 | endif |
|---|
| 228 | |
|---|
| 229 | ! close file |
|---|
| 230 | CALL close_restartphy |
|---|
| 231 | !$OMP BARRIER |
|---|
| 232 | |
|---|
| 233 | end subroutine physdem1 |
|---|
| 234 | |
|---|
| 235 | end module phyredem |
|---|
