[909] | 1 | MODULE phys_state_var_mod |
---|
| 2 | ! Variables sauvegardees pour le startphy.nc |
---|
| 3 | !====================================================================== |
---|
| 4 | ! |
---|
| 5 | ! |
---|
| 6 | !====================================================================== |
---|
| 7 | ! Declaration des variables |
---|
| 8 | USE dimphy |
---|
[952] | 9 | INTEGER, PARAMETER :: nlevSTD=17 |
---|
| 10 | INTEGER, PARAMETER :: nout=3 |
---|
| 11 | INTEGER, PARAMETER :: napisccp=1 |
---|
[967] | 12 | INTEGER, SAVE :: radpas |
---|
[1227] | 13 | REAL, SAVE :: dtime, solaire_etat0 |
---|
[967] | 14 | !$OMP THREADPRIVATE(radpas) |
---|
[1227] | 15 | !$OMP THREADPRIVATE(dtime, solaire_etat0) |
---|
[967] | 16 | |
---|
[909] | 17 | REAL, ALLOCATABLE, SAVE :: rlat(:), rlon(:), pctsrf(:,:) |
---|
[967] | 18 | !$OMP THREADPRIVATE(rlat, rlon, pctsrf) |
---|
[909] | 19 | REAL, ALLOCATABLE, SAVE :: ftsol(:,:) |
---|
[967] | 20 | !$OMP THREADPRIVATE(ftsol) |
---|
| 21 | ! character(len=6), SAVE :: ocean |
---|
| 22 | !!!!!!$OMP THREADPRIVATE(ocean) |
---|
| 23 | ! logical, SAVE :: ok_veget |
---|
| 24 | !!!!!!$OMP THREADPRIVATE(ok_veget) |
---|
[909] | 25 | REAL, ALLOCATABLE, SAVE :: falb1(:,:), falb2(:,:) |
---|
[967] | 26 | !$OMP THREADPRIVATE(falb1, falb2) |
---|
[909] | 27 | REAL, ALLOCATABLE, SAVE :: rain_fall(:), snow_fall(:) |
---|
[967] | 28 | !$OMP THREADPRIVATE( rain_fall, snow_fall) |
---|
[909] | 29 | REAL, ALLOCATABLE, SAVE :: solsw(:), sollw(:) |
---|
[967] | 30 | !$OMP THREADPRIVATE(solsw, sollw) |
---|
[909] | 31 | REAL, ALLOCATABLE, SAVE :: radsol(:) |
---|
[967] | 32 | !$OMP THREADPRIVATE(radsol) |
---|
| 33 | |
---|
[909] | 34 | !clesphy0 param physiq |
---|
| 35 | ! |
---|
| 36 | ! Parametres de l'Orographie a l'Echelle Sous-Maille (OESM): |
---|
| 37 | ! |
---|
| 38 | REAL, ALLOCATABLE, SAVE :: zmea(:), zstd(:), zsig(:), zgam(:) |
---|
[967] | 39 | !$OMP THREADPRIVATE(zmea, zstd, zsig, zgam) |
---|
[909] | 40 | REAL, ALLOCATABLE, SAVE :: zthe(:), zpic(:), zval(:) |
---|
[967] | 41 | !$OMP THREADPRIVATE(zthe, zpic, zval) |
---|
[909] | 42 | ! REAL tabcntr0(100) |
---|
| 43 | REAL, ALLOCATABLE, SAVE :: rugoro(:) |
---|
[967] | 44 | !$OMP THREADPRIVATE(rugoro) |
---|
[909] | 45 | REAL, ALLOCATABLE, SAVE :: t_ancien(:,:), q_ancien(:,:) |
---|
[967] | 46 | !$OMP THREADPRIVATE(t_ancien, q_ancien) |
---|
[1054] | 47 | REAL, ALLOCATABLE, SAVE :: u_ancien(:,:), v_ancien(:,:) |
---|
| 48 | !$OMP THREADPRIVATE(u_ancien, v_ancien) |
---|
[909] | 49 | LOGICAL, SAVE :: ancien_ok |
---|
[967] | 50 | !$OMP THREADPRIVATE(ancien_ok) |
---|
[909] | 51 | REAL, ALLOCATABLE, SAVE :: clwcon(:,:),rnebcon(:,:) |
---|
[967] | 52 | !$OMP THREADPRIVATE(clwcon,rnebcon) |
---|
[909] | 53 | REAL, ALLOCATABLE, SAVE :: ratqs(:,:) |
---|
[967] | 54 | !$OMP THREADPRIVATE(ratqs) |
---|
[909] | 55 | REAL, ALLOCATABLE, SAVE :: pbl_tke(:,:,:) ! turb kinetic energy |
---|
[967] | 56 | !$OMP THREADPRIVATE(pbl_tke) |
---|
[927] | 57 | REAL, ALLOCATABLE, SAVE :: zmax0(:), f0(:) ! |
---|
[967] | 58 | !$OMP THREADPRIVATE(zmax0,f0) |
---|
[937] | 59 | REAL, ALLOCATABLE, SAVE :: ema_work1(:,:), ema_work2(:,:) |
---|
[967] | 60 | !$OMP THREADPRIVATE(ema_work1,ema_work2) |
---|
[942] | 61 | REAL, ALLOCATABLE, SAVE :: entr_therm(:,:), fm_therm(:,:) |
---|
[967] | 62 | !$OMP THREADPRIVATE(entr_therm,fm_therm) |
---|
[952] | 63 | REAL, ALLOCATABLE, SAVE :: detr_therm(:,:) |
---|
[967] | 64 | !$OMP THREADPRIVATE(detr_therm) |
---|
[952] | 65 | !IM 150408 |
---|
| 66 | ! pour phsystoke avec thermiques |
---|
| 67 | REAL,ALLOCATABLE,SAVE :: clwcon0th(:,:),rnebcon0th(:,:) |
---|
[967] | 68 | !$OMP THREADPRIVATE(clwcon0th,rnebcon0th) |
---|
[952] | 69 | ! radiation outputs |
---|
| 70 | REAL,ALLOCATABLE,SAVE :: swdn0(:,:), swdn(:,:) |
---|
[967] | 71 | !$OMP THREADPRIVATE(swdn0,swdn) |
---|
[952] | 72 | REAL,ALLOCATABLE,SAVE :: swup0(:,:), swup(:,:) |
---|
[967] | 73 | !$OMP THREADPRIVATE(swup0,swup) |
---|
[952] | 74 | REAL,ALLOCATABLE,SAVE :: SWdn200clr(:), SWdn200(:) |
---|
[967] | 75 | !$OMP THREADPRIVATE(SWdn200clr,SWdn200) |
---|
[952] | 76 | REAL,ALLOCATABLE,SAVE :: SWup200clr(:), SWup200(:) |
---|
[967] | 77 | !$OMP THREADPRIVATE(SWup200clr,SWup200) |
---|
[952] | 78 | REAL,ALLOCATABLE,SAVE :: lwdn0(:,:), lwdn(:,:) |
---|
[967] | 79 | !$OMP THREADPRIVATE(lwdn0,lwdn) |
---|
[952] | 80 | REAL,ALLOCATABLE,SAVE :: lwup0(:,:), lwup(:,:) |
---|
[967] | 81 | !$OMP THREADPRIVATE(lwup0,lwup) |
---|
[952] | 82 | REAL,ALLOCATABLE,SAVE :: LWdn200clr(:), LWdn200(:) |
---|
[967] | 83 | !$OMP THREADPRIVATE(LWdn200clr,LWdn200) |
---|
[952] | 84 | REAL,ALLOCATABLE,SAVE :: LWup200clr(:), LWup200(:) |
---|
[967] | 85 | !$OMP THREADPRIVATE(LWup200clr,LWup200) |
---|
[952] | 86 | REAL,ALLOCATABLE,SAVE :: LWdnTOA(:), LWdnTOAclr(:) |
---|
[967] | 87 | !$OMP THREADPRIVATE(LWdnTOA,LWdnTOAclr) |
---|
[952] | 88 | ! pressure level |
---|
| 89 | REAL,ALLOCATABLE,SAVE :: tsumSTD(:,:,:) |
---|
[967] | 90 | !$OMP THREADPRIVATE(tsumSTD) |
---|
[952] | 91 | REAL,ALLOCATABLE,SAVE :: usumSTD(:,:,:), vsumSTD(:,:,:) |
---|
[967] | 92 | !$OMP THREADPRIVATE(usumSTD,vsumSTD) |
---|
[952] | 93 | REAL,ALLOCATABLE,SAVE :: wsumSTD(:,:,:), phisumSTD(:,:,:) |
---|
[967] | 94 | !$OMP THREADPRIVATE(wsumSTD,phisumSTD) |
---|
[952] | 95 | REAL,ALLOCATABLE,SAVE :: qsumSTD(:,:,:), rhsumSTD(:,:,:) |
---|
[967] | 96 | !$OMP THREADPRIVATE(qsumSTD,rhsumSTD) |
---|
[952] | 97 | REAL,ALLOCATABLE,SAVE :: tnondef(:,:,:) |
---|
[967] | 98 | !$OMP THREADPRIVATE(tnondef) |
---|
[952] | 99 | REAL,ALLOCATABLE,SAVE :: uvsumSTD(:,:,:) |
---|
[967] | 100 | !$OMP THREADPRIVATE(uvsumSTD) |
---|
[952] | 101 | REAL,ALLOCATABLE,SAVE :: vqsumSTD(:,:,:) |
---|
[967] | 102 | !$OMP THREADPRIVATE(vqsumSTD) |
---|
[952] | 103 | REAL,ALLOCATABLE,SAVE :: vTsumSTD(:,:,:) |
---|
[967] | 104 | !$OMP THREADPRIVATE(vTsumSTD) |
---|
[952] | 105 | REAL,ALLOCATABLE,SAVE :: wqsumSTD(:,:,:) |
---|
[967] | 106 | !$OMP THREADPRIVATE(wqsumSTD) |
---|
[952] | 107 | REAL,ALLOCATABLE,SAVE :: vphisumSTD(:,:,:) |
---|
[967] | 108 | !$OMP THREADPRIVATE(vphisumSTD) |
---|
[952] | 109 | REAL,ALLOCATABLE,SAVE :: wTsumSTD(:,:,:) |
---|
[967] | 110 | !$OMP THREADPRIVATE(wTsumSTD) |
---|
[952] | 111 | REAL,ALLOCATABLE,SAVE :: u2sumSTD(:,:,:) |
---|
[967] | 112 | !$OMP THREADPRIVATE(u2sumSTD) |
---|
[952] | 113 | REAL,ALLOCATABLE,SAVE :: v2sumSTD(:,:,:) |
---|
[967] | 114 | !$OMP THREADPRIVATE(v2sumSTD) |
---|
[952] | 115 | REAL,ALLOCATABLE,SAVE :: T2sumSTD(:,:,:) |
---|
[967] | 116 | !$OMP THREADPRIVATE(T2sumSTD) |
---|
[952] | 117 | INTEGER,ALLOCATABLE,SAVE :: seed_old(:,:) |
---|
[967] | 118 | !$OMP THREADPRIVATE(seed_old) |
---|
[952] | 119 | REAL,ALLOCATABLE,SAVE :: zuthe(:),zvthe(:) |
---|
[967] | 120 | !$OMP THREADPRIVATE(zuthe,zvthe) |
---|
[952] | 121 | REAL,ALLOCATABLE,SAVE :: alb_neig(:) |
---|
[967] | 122 | !$OMP THREADPRIVATE(alb_neig) |
---|
[952] | 123 | !cloud base mass flux |
---|
| 124 | REAL,ALLOCATABLE,SAVE :: ema_workcbmf(:), ema_cbmf(:) |
---|
[967] | 125 | !$OMP THREADPRIVATE(ema_workcbmf,ema_cbmf) |
---|
[952] | 126 | !cloud base pressure & cloud top pressure |
---|
| 127 | REAL,ALLOCATABLE,SAVE :: ema_pcb(:), ema_pct(:) |
---|
[967] | 128 | !$OMP THREADPRIVATE(ema_pcb,ema_pct) |
---|
[952] | 129 | REAL,ALLOCATABLE,SAVE :: Ma(:,:) ! undilute upward mass flux |
---|
[967] | 130 | !$OMP THREADPRIVATE(Ma) |
---|
[952] | 131 | REAL,ALLOCATABLE,SAVE :: qcondc(:,:) ! in-cld water content from convect |
---|
[967] | 132 | !$OMP THREADPRIVATE(qcondc) |
---|
[952] | 133 | REAL,ALLOCATABLE,SAVE :: wd(:) ! sb |
---|
[967] | 134 | !$OMP THREADPRIVATE(wd) |
---|
[952] | 135 | REAL,ALLOCATABLE,SAVE :: sigd(:) |
---|
[967] | 136 | !$OMP THREADPRIVATE(sigd) |
---|
[952] | 137 | ! |
---|
| 138 | REAL,ALLOCATABLE,SAVE :: cin(:) |
---|
[967] | 139 | !$OMP THREADPRIVATE(cin) |
---|
[952] | 140 | ! ftd : differential heating between wake and environment |
---|
| 141 | REAL,ALLOCATABLE,SAVE :: ftd(:,:) |
---|
[967] | 142 | !$OMP THREADPRIVATE(ftd) |
---|
[952] | 143 | ! fqd : differential moistening between wake and environment |
---|
| 144 | REAL,ALLOCATABLE,SAVE :: fqd(:,:) |
---|
[967] | 145 | !$OMP THREADPRIVATE(fqd) |
---|
[952] | 146 | !34EK |
---|
| 147 | ! -- Variables de controle de ALE et ALP |
---|
| 148 | !ALE : Energie disponible pour soulevement : utilisee par la |
---|
| 149 | ! convection d'Emanuel pour le declenchement et la regulation |
---|
| 150 | REAL,ALLOCATABLE,SAVE :: ALE(:) |
---|
[967] | 151 | !$OMP THREADPRIVATE(ALE) |
---|
[952] | 152 | !ALP : Puissance disponible pour soulevement |
---|
| 153 | REAL,ALLOCATABLE,SAVE :: ALP(:) |
---|
[967] | 154 | !$OMP THREADPRIVATE(ALP) |
---|
[952] | 155 | ! |
---|
| 156 | ! nouvelles variables pour le couplage convection-couche limite |
---|
| 157 | REAL,ALLOCATABLE,SAVE :: Ale_bl(:) |
---|
[967] | 158 | !$OMP THREADPRIVATE(Ale_bl) |
---|
[952] | 159 | REAL,ALLOCATABLE,SAVE :: Alp_bl(:) |
---|
[967] | 160 | !$OMP THREADPRIVATE(Alp_bl) |
---|
[952] | 161 | INTEGER,ALLOCATABLE,SAVE :: lalim_conv(:) |
---|
[967] | 162 | !$OMP THREADPRIVATE(lalim_conv) |
---|
[952] | 163 | REAL,ALLOCATABLE,SAVE :: wght_th(:,:) |
---|
[967] | 164 | !$OMP THREADPRIVATE(wght_th) |
---|
[952] | 165 | ! |
---|
| 166 | ! variables de la wake |
---|
| 167 | ! wake_deltat : ecart de temperature avec la zone non perturbee |
---|
| 168 | ! wake_deltaq : ecart d'humidite avec la zone non perturbee |
---|
| 169 | ! wake_Cstar : vitesse d'etalement de la poche |
---|
| 170 | ! wake_s : fraction surfacique occupee par la poche froide |
---|
| 171 | ! wake_fip : Gust Front Impinging power - ALP |
---|
| 172 | ! dt_wake, dq_wake: LS tendencies due to wake |
---|
| 173 | REAL,ALLOCATABLE,SAVE :: wake_deltat(:,:) |
---|
[967] | 174 | !$OMP THREADPRIVATE(wake_deltat) |
---|
[952] | 175 | REAL,ALLOCATABLE,SAVE :: wake_deltaq(:,:) |
---|
[967] | 176 | !$OMP THREADPRIVATE(wake_deltaq) |
---|
[952] | 177 | REAL,ALLOCATABLE,SAVE :: wake_Cstar(:) |
---|
[967] | 178 | !$OMP THREADPRIVATE(wake_Cstar) |
---|
[952] | 179 | REAL,ALLOCATABLE,SAVE :: wake_s(:) |
---|
[967] | 180 | !$OMP THREADPRIVATE(wake_s) |
---|
[952] | 181 | REAL,ALLOCATABLE,SAVE :: wake_fip(:) |
---|
[967] | 182 | !$OMP THREADPRIVATE(wake_fip) |
---|
[952] | 183 | REAL,ALLOCATABLE,SAVE :: dt_wake(:,:) |
---|
[967] | 184 | !$OMP THREADPRIVATE(dt_wake) |
---|
[952] | 185 | REAL,ALLOCATABLE,SAVE :: dq_wake(:,:) |
---|
[967] | 186 | !$OMP THREADPRIVATE(dq_wake) |
---|
[952] | 187 | ! |
---|
| 188 | ! pfrac_impa : Produits des coefs lessivage impaction |
---|
| 189 | ! pfrac_nucl : Produits des coefs lessivage nucleation |
---|
| 190 | ! pfrac_1nucl: Produits des coefs lessi nucl (alpha = 1) |
---|
| 191 | REAL,ALLOCATABLE,SAVE :: pfrac_impa(:,:), pfrac_nucl(:,:) |
---|
[967] | 192 | !$OMP THREADPRIVATE(pfrac_impa,pfrac_nucl) |
---|
[952] | 193 | REAL,ALLOCATABLE,SAVE :: pfrac_1nucl(:,:) |
---|
[967] | 194 | !$OMP THREADPRIVATE(pfrac_1nucl) |
---|
[952] | 195 | ! |
---|
| 196 | REAL,ALLOCATABLE,SAVE :: total_rain(:), nday_rain(:) |
---|
[967] | 197 | !$OMP THREADPRIVATE(total_rain,nday_rain) |
---|
[952] | 198 | REAL,ALLOCATABLE,SAVE :: paire_ter(:) |
---|
[967] | 199 | !$OMP THREADPRIVATE(paire_ter) |
---|
[952] | 200 | ! albsol1: albedo du sol total pour SW visible |
---|
| 201 | ! albsol2: albedo du sol total pour SW proche IR |
---|
| 202 | REAL,ALLOCATABLE,SAVE :: albsol1(:), albsol2(:) |
---|
[967] | 203 | !$OMP THREADPRIVATE(albsol1,albsol2) |
---|
[1215] | 204 | |
---|
| 205 | REAL, ALLOCATABLE, SAVE:: wo(:,:) |
---|
| 206 | ! column-density of ozone in a layer, in kilo-Dobsons |
---|
| 207 | !$OMP THREADPRIVATE(wo) |
---|
| 208 | |
---|
[952] | 209 | ! heat : chauffage solaire |
---|
| 210 | ! heat0: chauffage solaire ciel clair |
---|
| 211 | ! cool : refroidissement infrarouge |
---|
| 212 | ! cool0 : refroidissement infrarouge ciel clair |
---|
| 213 | ! sollwdown : downward LW flux at surface |
---|
| 214 | ! sollwdownclr : downward CS LW flux at surface |
---|
| 215 | ! toplwdown : downward CS LW flux at TOA |
---|
| 216 | ! toplwdownclr : downward CS LW flux at TOA |
---|
| 217 | REAL,ALLOCATABLE,SAVE :: clwcon0(:,:),rnebcon0(:,:) |
---|
[967] | 218 | !$OMP THREADPRIVATE(clwcon0,rnebcon0) |
---|
[952] | 219 | REAL,ALLOCATABLE,SAVE :: heat(:,:) |
---|
[967] | 220 | !$OMP THREADPRIVATE(heat) |
---|
[952] | 221 | REAL,ALLOCATABLE,SAVE :: heat0(:,:) |
---|
[967] | 222 | !$OMP THREADPRIVATE(heat0) |
---|
[952] | 223 | REAL,ALLOCATABLE,SAVE :: cool(:,:) |
---|
[967] | 224 | !$OMP THREADPRIVATE(cool) |
---|
[952] | 225 | REAL,ALLOCATABLE,SAVE :: cool0(:,:) |
---|
[967] | 226 | !$OMP THREADPRIVATE(cool0) |
---|
[952] | 227 | REAL,ALLOCATABLE,SAVE :: topsw(:), toplw(:) |
---|
[967] | 228 | !$OMP THREADPRIVATE(topsw,toplw) |
---|
[952] | 229 | REAL,ALLOCATABLE,SAVE :: sollwdown(:) |
---|
[967] | 230 | !$OMP THREADPRIVATE(sollwdown) |
---|
[952] | 231 | REAL,ALLOCATABLE,SAVE :: sollwdownclr(:) |
---|
[967] | 232 | !$OMP THREADPRIVATE(sollwdownclr) |
---|
[952] | 233 | REAL,ALLOCATABLE,SAVE :: toplwdown(:) |
---|
[967] | 234 | !$OMP THREADPRIVATE(toplwdown) |
---|
[952] | 235 | REAL,ALLOCATABLE,SAVE :: toplwdownclr(:) |
---|
[967] | 236 | !$OMP THREADPRIVATE(toplwdownclr) |
---|
[952] | 237 | REAL,ALLOCATABLE,SAVE :: topsw0(:),toplw0(:),solsw0(:),sollw0(:) |
---|
[967] | 238 | !$OMP THREADPRIVATE(topsw0,toplw0,solsw0,sollw0) |
---|
[952] | 239 | REAL,ALLOCATABLE,SAVE :: albpla(:) |
---|
[967] | 240 | !$OMP THREADPRIVATE(albpla) |
---|
[952] | 241 | ! pbase : cloud base pressure |
---|
| 242 | ! bbase : cloud base buoyancy |
---|
| 243 | REAL,ALLOCATABLE,SAVE :: cape(:) |
---|
[967] | 244 | !$OMP THREADPRIVATE(cape) |
---|
[952] | 245 | REAL,ALLOCATABLE,SAVE :: pbase(:) |
---|
[967] | 246 | !$OMP THREADPRIVATE(pbase) |
---|
[952] | 247 | REAL,ALLOCATABLE,SAVE :: bbase(:) |
---|
[967] | 248 | !$OMP THREADPRIVATE(bbase) |
---|
[952] | 249 | ! |
---|
| 250 | REAL,SAVE,ALLOCATABLE :: zqasc(:,:) |
---|
[967] | 251 | !$OMP THREADPRIVATE( zqasc) |
---|
[952] | 252 | INTEGER,ALLOCATABLE,SAVE :: ibas_con(:), itop_con(:) |
---|
[967] | 253 | !$OMP THREADPRIVATE(ibas_con,itop_con) |
---|
[952] | 254 | REAL,SAVE,ALLOCATABLE :: rain_con(:) |
---|
[967] | 255 | !$OMP THREADPRIVATE(rain_con) |
---|
[952] | 256 | REAL,SAVE,ALLOCATABLE :: snow_con(:) |
---|
[967] | 257 | !$OMP THREADPRIVATE(snow_con) |
---|
[952] | 258 | ! |
---|
| 259 | REAL,SAVE,ALLOCATABLE :: rlonPOS(:) |
---|
[967] | 260 | !$OMP THREADPRIVATE(rlonPOS) |
---|
[952] | 261 | REAL,SAVE,ALLOCATABLE :: newsst(:) |
---|
[967] | 262 | !$OMP THREADPRIVATE(newsst) |
---|
[952] | 263 | REAL,SAVE,ALLOCATABLE :: u10m(:,:), v10m(:,:) |
---|
[967] | 264 | !$OMP THREADPRIVATE(u10m,v10m) |
---|
[952] | 265 | ! |
---|
| 266 | ! ok_ade=T -ADE=topswad-topsw |
---|
| 267 | ! ok_aie=T -> |
---|
| 268 | ! ok_ade=T -AIE=topswai-topswad |
---|
[1150] | 269 | ! ok_ade=F -AIE=topswai-topsw |
---|
[952] | 270 | ! |
---|
| 271 | !topswad, solswad : Aerosol direct effect |
---|
| 272 | REAL,SAVE,ALLOCATABLE :: topswad(:), solswad(:) |
---|
[967] | 273 | !$OMP THREADPRIVATE(topswad,solswad) |
---|
[952] | 274 | !topswai, solswai : Aerosol indirect effect |
---|
| 275 | REAL,SAVE,ALLOCATABLE :: topswai(:), solswai(:) |
---|
[967] | 276 | !$OMP THREADPRIVATE(topswai,solswai) |
---|
[959] | 277 | |
---|
[1179] | 278 | REAL,SAVE,ALLOCATABLE :: tau_aero(:,:,:,:), piz_aero(:,:,:,:), cg_aero(:,:,:,:) |
---|
| 279 | !$OMP THREADPRIVATE(tau_aero, piz_aero, cg_aero) |
---|
[952] | 280 | REAL,SAVE,ALLOCATABLE :: ccm(:,:,:) |
---|
[967] | 281 | !$OMP THREADPRIVATE(ccm) |
---|
[952] | 282 | |
---|
[909] | 283 | CONTAINS |
---|
| 284 | |
---|
| 285 | !====================================================================== |
---|
| 286 | SUBROUTINE phys_state_var_init |
---|
| 287 | use dimphy |
---|
[1183] | 288 | use aero_mod |
---|
[909] | 289 | IMPLICIT NONE |
---|
| 290 | #include "indicesol.h" |
---|
[959] | 291 | #include "control.h" |
---|
[952] | 292 | ALLOCATE(rlat(klon), rlon(klon)) |
---|
| 293 | ALLOCATE(pctsrf(klon,nbsrf)) |
---|
| 294 | ALLOCATE(ftsol(klon,nbsrf)) |
---|
| 295 | ALLOCATE(falb1(klon,nbsrf)) |
---|
| 296 | ALLOCATE(falb2(klon,nbsrf)) |
---|
| 297 | ALLOCATE(rain_fall(klon)) |
---|
| 298 | ALLOCATE(snow_fall(klon)) |
---|
[909] | 299 | ALLOCATE(solsw(klon), sollw(klon)) |
---|
[952] | 300 | ALLOCATE(radsol(klon)) |
---|
| 301 | ALLOCATE(zmea(klon), zstd(klon), zsig(klon), zgam(klon)) |
---|
| 302 | ALLOCATE(zthe(klon), zpic(klon), zval(klon)) |
---|
[909] | 303 | |
---|
[952] | 304 | ALLOCATE(rugoro(klon)) |
---|
| 305 | ALLOCATE(t_ancien(klon,klev), q_ancien(klon,klev)) |
---|
[1054] | 306 | ALLOCATE(u_ancien(klon,klev), v_ancien(klon,klev)) |
---|
[952] | 307 | ALLOCATE(clwcon(klon,klev),rnebcon(klon,klev)) |
---|
| 308 | ALLOCATE(ratqs(klon,klev)) |
---|
[909] | 309 | ALLOCATE(pbl_tke(klon,klev+1,nbsrf)) |
---|
[927] | 310 | ALLOCATE(zmax0(klon), f0(klon)) |
---|
[937] | 311 | ALLOCATE(ema_work1(klon,klev), ema_work2(klon,klev)) |
---|
[942] | 312 | ALLOCATE(entr_therm(klon,klev), fm_therm(klon,klev+1)) |
---|
[952] | 313 | ALLOCATE(detr_therm(klon,klev)) |
---|
| 314 | ! pour phsystoke avec thermiques |
---|
| 315 | ALLOCATE(clwcon0th(klon,klev),rnebcon0th(klon,klev)) |
---|
| 316 | ! radiation outputs |
---|
| 317 | ALLOCATE(swdn0(klon,klevp1), swdn(klon,klevp1)) |
---|
| 318 | ALLOCATE(swup0(klon,klevp1), swup(klon,klevp1)) |
---|
| 319 | ALLOCATE(lwdn0(klon,klevp1), lwdn(klon,klevp1)) |
---|
| 320 | ALLOCATE(lwup0(klon,klevp1), lwup(klon,klevp1)) |
---|
| 321 | ALLOCATE(SWdn200clr(klon), SWdn200(klon)) |
---|
| 322 | ALLOCATE(SWup200clr(klon), SWup200(klon)) |
---|
| 323 | ALLOCATE(LWdn200clr(klon), LWdn200(klon)) |
---|
| 324 | ALLOCATE(LWup200clr(klon), LWup200(klon)) |
---|
| 325 | ALLOCATE(LWdnTOA(klon), LWdnTOAclr(klon)) |
---|
| 326 | ! pressure level |
---|
| 327 | ALLOCATE(tsumSTD(klon,nlevSTD,nout)) |
---|
| 328 | ALLOCATE(usumSTD(klon,nlevSTD,nout), vsumSTD(klon,nlevSTD,nout)) |
---|
| 329 | ALLOCATE(wsumSTD(klon,nlevSTD,nout), phisumSTD(klon,nlevSTD,nout)) |
---|
| 330 | ALLOCATE(qsumSTD(klon,nlevSTD,nout), rhsumSTD(klon,nlevSTD,nout)) |
---|
| 331 | ALLOCATE(tnondef(klon,nlevSTD,nout)) |
---|
| 332 | ALLOCATE(uvsumSTD(klon,nlevSTD,nout)) |
---|
| 333 | ALLOCATE(vqsumSTD(klon,nlevSTD,nout)) |
---|
| 334 | ALLOCATE(vTsumSTD(klon,nlevSTD,nout)) |
---|
| 335 | ALLOCATE(wqsumSTD(klon,nlevSTD,nout)) |
---|
| 336 | ALLOCATE(vphisumSTD(klon,nlevSTD,nout)) |
---|
| 337 | ALLOCATE(wTsumSTD(klon,nlevSTD,nout)) |
---|
| 338 | ALLOCATE(u2sumSTD(klon,nlevSTD,nout)) |
---|
| 339 | ALLOCATE(v2sumSTD(klon,nlevSTD,nout)) |
---|
| 340 | ALLOCATE(T2sumSTD(klon,nlevSTD,nout)) |
---|
| 341 | ALLOCATE(seed_old(klon,napisccp)) |
---|
| 342 | ALLOCATE(zuthe(klon),zvthe(klon)) |
---|
| 343 | ALLOCATE(alb_neig(klon)) |
---|
| 344 | !cloud base mass flux |
---|
| 345 | ALLOCATE(ema_workcbmf(klon), ema_cbmf(klon)) |
---|
| 346 | !cloud base pressure & cloud top pressure |
---|
| 347 | ALLOCATE(ema_pcb(klon), ema_pct(klon)) |
---|
| 348 | ! |
---|
| 349 | ALLOCATE(Ma(klon,klev)) |
---|
| 350 | ALLOCATE(qcondc(klon,klev)) |
---|
| 351 | ALLOCATE(wd(klon)) |
---|
| 352 | ALLOCATE(sigd(klon)) |
---|
| 353 | ALLOCATE(cin(klon), ALE(klon), ALP(klon)) |
---|
| 354 | ALLOCATE(ftd(klon,klev), fqd(klon,klev)) |
---|
| 355 | ALLOCATE(Ale_bl(klon)) |
---|
| 356 | ALLOCATE(Alp_bl(klon)) |
---|
| 357 | ALLOCATE(lalim_conv(klon)) |
---|
| 358 | ALLOCATE(wght_th(klon,klev)) |
---|
| 359 | ALLOCATE(wake_deltat(klon,klev), wake_deltaq(klon,klev)) |
---|
| 360 | ALLOCATE(wake_Cstar(klon), wake_s(klon), wake_fip(klon)) |
---|
| 361 | ALLOCATE(dt_wake(klon,klev), dq_wake(klon,klev)) |
---|
| 362 | ALLOCATE(pfrac_impa(klon,klev), pfrac_nucl(klon,klev)) |
---|
| 363 | ALLOCATE(pfrac_1nucl(klon,klev)) |
---|
| 364 | ALLOCATE(total_rain(klon), nday_rain(klon)) |
---|
| 365 | ALLOCATE(paire_ter(klon)) |
---|
| 366 | ALLOCATE(albsol1(klon), albsol2(klon)) |
---|
| 367 | ALLOCATE(wo(klon,klev)) |
---|
| 368 | ALLOCATE(clwcon0(klon,klev),rnebcon0(klon,klev)) |
---|
| 369 | ALLOCATE(heat(klon,klev), heat0(klon,klev)) |
---|
| 370 | ALLOCATE(cool(klon,klev), cool0(klon,klev)) |
---|
| 371 | ALLOCATE(topsw(klon), toplw(klon)) |
---|
| 372 | ALLOCATE(sollwdown(klon), sollwdownclr(klon)) |
---|
| 373 | ALLOCATE(toplwdown(klon), toplwdownclr(klon)) |
---|
| 374 | ALLOCATE(topsw0(klon),toplw0(klon),solsw0(klon),sollw0(klon)) |
---|
| 375 | ALLOCATE(albpla(klon)) |
---|
| 376 | ALLOCATE(cape(klon)) |
---|
| 377 | ALLOCATE(pbase(klon),bbase(klon)) |
---|
| 378 | ALLOCATE(zqasc(klon,klev)) |
---|
| 379 | ALLOCATE(ibas_con(klon), itop_con(klon)) |
---|
| 380 | ALLOCATE(rain_con(klon), snow_con(klon)) |
---|
| 381 | ALLOCATE(rlonPOS(klon)) |
---|
| 382 | ALLOCATE(newsst(klon)) |
---|
| 383 | ALLOCATE(u10m(klon,nbsrf), v10m(klon,nbsrf)) |
---|
| 384 | ALLOCATE(topswad(klon), solswad(klon)) |
---|
| 385 | ALLOCATE(topswai(klon), solswai(klon)) |
---|
[1181] | 386 | ALLOCATE(tau_aero(klon,klev,naero_grp,nbands),piz_aero(klon,klev,naero_grp,nbands),cg_aero(klon,klev,naero_grp,nbands)) |
---|
| 387 | ALLOCATE(ccm(klon,klev,nbands)) |
---|
[909] | 388 | |
---|
| 389 | END SUBROUTINE phys_state_var_init |
---|
| 390 | |
---|
| 391 | !====================================================================== |
---|
| 392 | SUBROUTINE phys_state_var_end |
---|
| 393 | use dimphy |
---|
| 394 | IMPLICIT NONE |
---|
| 395 | #include "indicesol.h" |
---|
[959] | 396 | #include "control.h" |
---|
[909] | 397 | |
---|
[952] | 398 | deallocate(rlat, rlon, pctsrf, ftsol, falb1, falb2) |
---|
| 399 | deallocate(rain_fall, snow_fall, solsw, sollw, radsol) |
---|
| 400 | deallocate(zmea, zstd, zsig, zgam) |
---|
| 401 | deallocate(zthe, zpic, zval) |
---|
| 402 | deallocate(rugoro, t_ancien, q_ancien, clwcon, rnebcon) |
---|
[1054] | 403 | deallocate( u_ancien, v_ancien ) |
---|
[952] | 404 | deallocate(ratqs, pbl_tke) |
---|
| 405 | deallocate(zmax0, f0) |
---|
| 406 | deallocate(ema_work1, ema_work2) |
---|
| 407 | deallocate(entr_therm, fm_therm) |
---|
| 408 | deallocate(detr_therm) |
---|
| 409 | deallocate(clwcon0th, rnebcon0th) |
---|
| 410 | ! radiation outputs |
---|
| 411 | deallocate(swdn0, swdn) |
---|
| 412 | deallocate(swup0, swup) |
---|
| 413 | deallocate(lwdn0, lwdn) |
---|
| 414 | deallocate(lwup0, lwup) |
---|
| 415 | deallocate(SWdn200clr, SWdn200) |
---|
| 416 | deallocate(SWup200clr, SWup200) |
---|
| 417 | deallocate(LWdn200clr, LWdn200) |
---|
| 418 | deallocate(LWup200clr, LWup200) |
---|
| 419 | deallocate(LWdnTOA, LWdnTOAclr) |
---|
| 420 | ! pressure level |
---|
| 421 | deallocate(tsumSTD) |
---|
| 422 | deallocate(usumSTD, vsumSTD) |
---|
| 423 | deallocate(wsumSTD, phisumSTD) |
---|
| 424 | deallocate(tnondef) |
---|
[970] | 425 | deallocate(qsumSTD, rhsumSTD) |
---|
[952] | 426 | deallocate(uvsumSTD) |
---|
| 427 | deallocate(vqsumSTD) |
---|
| 428 | deallocate(vTsumSTD) |
---|
| 429 | deallocate(wqsumSTD) |
---|
| 430 | deallocate(vphisumSTD) |
---|
| 431 | deallocate(wTsumSTD) |
---|
| 432 | deallocate(u2sumSTD) |
---|
| 433 | deallocate(v2sumSTD) |
---|
| 434 | deallocate(T2sumSTD) |
---|
| 435 | deallocate(seed_old) |
---|
| 436 | deallocate(zuthe, zvthe) |
---|
| 437 | deallocate(alb_neig) |
---|
| 438 | deallocate(ema_workcbmf, ema_cbmf) |
---|
| 439 | deallocate(ema_pcb, ema_pct) |
---|
| 440 | deallocate(Ma, qcondc) |
---|
| 441 | deallocate(wd, sigd) |
---|
| 442 | deallocate(cin, ALE, ALP) |
---|
| 443 | deallocate(ftd, fqd) |
---|
| 444 | deallocate(Ale_bl, Alp_bl) |
---|
| 445 | deallocate(lalim_conv, wght_th) |
---|
| 446 | deallocate(wake_deltat, wake_deltaq) |
---|
| 447 | deallocate(wake_Cstar, wake_s, wake_fip) |
---|
| 448 | deallocate(dt_wake, dq_wake) |
---|
| 449 | deallocate(pfrac_impa, pfrac_nucl) |
---|
| 450 | deallocate(pfrac_1nucl) |
---|
| 451 | deallocate(total_rain, nday_rain) |
---|
| 452 | deallocate(paire_ter) |
---|
| 453 | deallocate(albsol1, albsol2) |
---|
| 454 | deallocate(wo) |
---|
| 455 | deallocate(clwcon0,rnebcon0) |
---|
| 456 | deallocate(heat, heat0) |
---|
| 457 | deallocate(cool, cool0) |
---|
| 458 | deallocate(topsw, toplw) |
---|
| 459 | deallocate(sollwdown, sollwdownclr) |
---|
| 460 | deallocate(toplwdown, toplwdownclr) |
---|
| 461 | deallocate(topsw0,toplw0,solsw0,sollw0) |
---|
| 462 | deallocate(albpla) |
---|
| 463 | deallocate(cape) |
---|
| 464 | deallocate(pbase,bbase) |
---|
| 465 | deallocate(zqasc) |
---|
| 466 | deallocate(ibas_con, itop_con) |
---|
| 467 | deallocate(rain_con, snow_con) |
---|
| 468 | deallocate(rlonPOS) |
---|
| 469 | deallocate(newsst) |
---|
| 470 | deallocate(u10m, v10m) |
---|
| 471 | deallocate(topswad, solswad) |
---|
| 472 | deallocate(topswai, solswai) |
---|
[1179] | 473 | deallocate(tau_aero,piz_aero,cg_aero) |
---|
[952] | 474 | deallocate(ccm) |
---|
[959] | 475 | |
---|
[909] | 476 | END SUBROUTINE phys_state_var_end |
---|
| 477 | |
---|
| 478 | END MODULE phys_state_var_mod |
---|