- Timestamp:
- Mar 12, 2025, 10:44:50 AM (13 months ago)
- Location:
- LMDZ6/branches/contrails
- Files:
-
- 9 edited
-
DefLists/context_input_lmdz.xml (modified) (4 diffs)
-
DefLists/file_def_histday_lmdz.xml (modified) (1 diff)
-
DefLists/file_def_histhf_lmdz.xml (modified) (1 diff)
-
DefLists/file_def_histmth_lmdz.xml (modified) (1 diff)
-
libf/phylmd/lmdz_aviation.f90 (modified) (3 diffs)
-
libf/phylmd/lmdz_lscp.f90 (modified) (5 diffs)
-
libf/phylmd/phys_output_ctrlout_mod.F90 (modified) (1 diff)
-
libf/phylmd/phys_state_var_mod.F90 (modified) (2 diffs)
-
libf/phylmd/physiq_mod.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/contrails/DefLists/context_input_lmdz.xml
r5546 r5573 399 399 <file id="aviation_file" name="aviation" mode="read" output_freq="1mo" type="one_file" time_counter_name="toto" > 400 400 401 <field id="KMFLOWN_id" name=" KMFLOWN" operation="instant" grid_ref="aviation_grid" freq_offset="1ts" />402 <field id="levaviation_id" name=" LEV" axis_ref="aviation_lev" operation="instant" freq_offset="1ts" />403 <field id="timeaviation_id" name=" TIME" axis_ref="aviation_time" operation="instant" freq_offset="1ts" />401 <field id="KMFLOWN_id" name="seg_length_km" operation="instant" grid_ref="aviation_grid" freq_offset="1ts" /> 402 <field id="levaviation_id" name="pressure_Pa" axis_ref="aviation_lev" operation="instant" freq_offset="1ts" /> 403 <field id="timeaviation_id" name="time" axis_ref="aviation_time" operation="instant" freq_offset="1ts" /> 404 404 405 405 </file> … … 427 427 <axis_definition> 428 428 429 <axis id="aviation_lev" unit="hPa" axis_type="Z"/>430 <axis id="aviation_time" name=" TIME" />429 <axis id="aviation_lev" axis_type="Z" /> 430 <axis id="aviation_time" name="time" /> 431 431 432 432 </axis_definition> … … 435 435 436 436 <field id ="KMFLOWN_read" field_ref="KMFLOWN_id" read_access="true" /> 437 <field id ="levaviation_read" field_ref="levaviation_id" read_access="true" />438 <field id ="timeaviation_read" field_ref="timeaviation_id" read_access="true" />439 440 437 <field id ="KMFLOWN_interp" field_ref="KMFLOWN_read" read_access="true" grid_ref="grid_from_aviation" /> 441 438 … … 454 451 <grid id="grid_from_aviation"> 455 452 <domain domain_ref="dom_from_aviation" /> 456 <axis axis_ref=" klev" />453 <axis axis_ref="aviation_lev" /> 457 454 <axis axis_ref="aviation_time" /> 458 455 </grid> -
LMDZ6/branches/contrails/DefLists/file_def_histday_lmdz.xml
r5343 r5573 504 504 <!-- VARS 3D --> 505 505 <field_group operation="average" grid_ref="grid_out_presnivs"> 506 <field field_ref="flightdist" level="1" /> 507 <field field_ref="cfseri" level="1" /> 508 <field field_ref="qissr" level="1" /> 509 <field field_ref="qcld" level="1" /> 510 <field field_ref="subfra" level="1" /> 511 <field field_ref="issrfra" level="1" /> 512 <field field_ref="dqvcmix" level="1" /> 513 <field field_ref="dqimix" level="1" /> 514 <field field_ref="dcfmix" level="1" /> 506 515 <field field_ref="CO2" level="5" /> <!-- Added PC --> 507 516 <field field_ref="dCO2_vdf" level="5" /> -
LMDZ6/branches/contrails/DefLists/file_def_histhf_lmdz.xml
r4846 r5573 1 <!-- $Id$ -->2 <file_definition>3 <file_group id="defile">4 5 6 7 <!-- <field_group operation="once"> -->8 <!-- <field field_ref="Ap" level="1" /> -->9 <!-- <field field_ref="Bp" level="1" /> -->10 <!-- <field field_ref="Alt" level="1" /> -->11 <!-- </field_group> -->12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 <field field_ref="snow" level="5" />70 <field field_ref="bsfall" level="5" />71 <field field_ref="evap" level="5" />72 <field field_ref="snowerosion" level="5" />73 <field field_ref="rhosnow_lic" level="5" />74 75 <field field_ref="evap_ter" level="10" />76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 <field field_ref="prsw" level="10" />223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 1 <!-- $Id$ --> 2 <file_definition> 3 <file_group id="defile"> 4 <file id="histhf" name="histhf" output_freq="1h" output_level="_AUTO_" enabled="_AUTO_" compression_level="2" sync_freq="1h"> 5 6 <!-- VARS 1D --> 7 <!-- <field_group operation="once"> --> 8 <!-- <field field_ref="Ap" level="1" /> --> 9 <!-- <field field_ref="Bp" level="1" /> --> 10 <!-- <field field_ref="Alt" level="1" /> --> 11 <!-- </field_group> --> 12 13 <!-- VARS 2D --> 14 <field_group group_ref="remap_1h" > 15 <field_group grid_ref="grid_out" > 16 <field field_ref="phis" level="5" /> 17 <field field_ref="aire" level="10" operation="once" /> 18 <field field_ref="contfracATM" level="5" operation="once" /> 19 <field field_ref="contfracOR" level="5" /> 20 <field field_ref="flat" level="5" /> 21 <field field_ref="slp" level="5" /> 22 <field field_ref="sza" level="10" /> 23 <field field_ref="ptstar" level="10" /> 24 <field field_ref="pt0" level="10" /> 25 <field field_ref="tsol" level="5" /> 26 <field field_ref="t2m" level="5" /> 27 <field field_ref="t2m_min" level="10" operation="minimum" /> 28 <field field_ref="t2m_max" level="10" operation="maximum" /> 29 <field field_ref="t2m_ter" level="10" /> 30 <field field_ref="t2m_lic" level="10" /> 31 <field field_ref="t2m_oce" level="10" /> 32 <field field_ref="t2m_sic" level="10" /> 33 <field field_ref="gusts" level="10" /> 34 <field field_ref="wind10m" level="5" /> 35 <field field_ref="wind100m" level="10" /> 36 <field field_ref="wind10max" level="10" operation="maximum" /> 37 <field field_ref="sicf" level="10" /> 38 <field field_ref="q2m" level="5" /> 39 <field field_ref="ustar" level="5" /> 40 <field field_ref="u10m" level="5" /> 41 <field field_ref="v10m" level="5" /> 42 <field field_ref="psol" level="5" /> 43 <field field_ref="qsurf" level="5" /> 44 <field field_ref="ustar_ter" level="10" /> 45 <field field_ref="ustar_lic" level="10" /> 46 <field field_ref="ustar_oce" level="10" /> 47 <field field_ref="ustar_sic" level="10" /> 48 <field field_ref="wstar_ter" level="10" /> 49 <field field_ref="wstar_lic" level="10" /> 50 <field field_ref="wstar_oce" level="10" /> 51 <field field_ref="wstar_sic" level="10" /> 52 <field field_ref="wstar" level="10" /> 53 <field field_ref="u10m_ter" level="10" /> 54 <field field_ref="u10m_lic" level="10" /> 55 <field field_ref="u10m_oce" level="10" /> 56 <field field_ref="u10m_sic" level="10" /> 57 <field field_ref="v10m_ter" level="10" /> 58 <field field_ref="v10m_lic" level="10" /> 59 <field field_ref="v10m_oce" level="10" /> 60 <field field_ref="v10m_sic" level="10" /> 61 <field field_ref="qsol" level="10" /> 62 <field field_ref="ndayrain" level="10" operation="instant" /> 63 <field field_ref="precip" level="5" /> 64 <field field_ref="rain_fall" level="10" /> 65 <field field_ref="rain_con" level="10" /> 66 <field field_ref="plul" level="10" /> 67 <field field_ref="plun" level="10" /> 68 <field field_ref="pluc" level="10" /> 69 <field field_ref="snow" level="5" /> 70 <field field_ref="bsfall" level="5" /> 71 <field field_ref="evap" level="5" /> 72 <field field_ref="snowerosion" level="5" /> 73 <field field_ref="rhosnow_lic" level="5" /> 74 <field field_ref="ustart_lic" level="5" /> 75 <field field_ref="evap_ter" level="10" /> 76 <field field_ref="evap_lic" level="10" /> 77 <field field_ref="evap_oce" level="10" /> 78 <field field_ref="evap_sic" level="10" /> 79 <field field_ref="msnow" level="5" /> 80 <field field_ref="fsnow" level="10" /> 81 <field field_ref="tops" level="10" /> 82 <field field_ref="tops0" level="10" /> 83 <field field_ref="topl" level="10" /> 84 <field field_ref="topl0" level="10" /> 85 <field field_ref="SWupTOA" level="10" /> 86 <field field_ref="SWupTOAclr" level="10" /> 87 <field field_ref="SWdnTOA" level="10" /> 88 <field field_ref="SWdnTOAclr" level="10" /> 89 <field field_ref="nettop" level="10" /> 90 <field field_ref="SWup200" level="15" /> 91 <field field_ref="SWup200clr" level="15" /> 92 <field field_ref="SWdn200" level="15" /> 93 <field field_ref="SWdn200clr" level="15" /> 94 <field field_ref="LWup200" level="15" /> 95 <field field_ref="LWup200clr" level="15" /> 96 <field field_ref="LWdn200" level="15" /> 97 <field field_ref="LWdn200clr" level="15" /> 98 <field field_ref="sols" level="10" /> 99 <field field_ref="sols0" level="10" /> 100 <field field_ref="soll" level="10" /> 101 <field field_ref="soll0" level="10" /> 102 <field field_ref="radsol" level="10" /> 103 <field field_ref="SWupSFC" level="10" /> 104 <field field_ref="SWupSFCclr" level="10" /> 105 <field field_ref="SWdnSFC" level="5" /> 106 <field field_ref="SWdnSFCclr" level="10" /> 107 <field field_ref="LWupSFC" level="10" /> 108 <field field_ref="LWupSFCclr" level="10" /> 109 <field field_ref="LWdnSFC" level="5" /> 110 <field field_ref="LWdnSFCclr" level="10" /> 111 <field field_ref="bils" level="10" /> 112 <field field_ref="bils_tke" level="10" /> 113 <field field_ref="bils_diss" level="10" /> 114 <field field_ref="bils_ec" level="10" /> 115 <field field_ref="bils_kinetic" level="10" /> 116 <field field_ref="bils_enthalp" level="10" /> 117 <field field_ref="bils_latent" level="10" /> 118 <field field_ref="sens" level="5" /> 119 <field field_ref="sens_x" level="11" /> 120 <field field_ref="sens_w" level="11" /> 121 <field field_ref="flat_x" level="11" /> 122 <field field_ref="flat_w" level="11" /> 123 <field field_ref="delta_tsurf" level="11" /> 124 <field field_ref="cdragh_x" level="11" /> 125 <field field_ref="cdragh_w" level="11" /> 126 <field field_ref="cdragm_x" level="11" /> 127 <field field_ref="cdragm_w" level="11" /> 128 <field field_ref="kh" level="11" /> 129 <field field_ref="kh_x" level="11" /> 130 <field field_ref="kh_w" level="11" /> 131 <field field_ref="fder" level="10" /> 132 <field field_ref="ffonte" level="10" /> 133 <field field_ref="fqcalving" level="10" /> 134 <field field_ref="fqfonte" level="10" /> 135 <field field_ref="mrroli" level="10" /> 136 <field field_ref="runofflic" level="10" /> 137 <field field_ref="taux" level="10" /> 138 <field field_ref="tauy" level="10" /> 139 <field field_ref="taux_ter" level="10" /> 140 <field field_ref="taux_lic" level="10" /> 141 <field field_ref="taux_oce" level="10" /> 142 <field field_ref="taux_sic" level="10" /> 143 <field field_ref="tauy_ter" level="10" /> 144 <field field_ref="tauy_lic" level="10" /> 145 <field field_ref="tauy_oce" level="10" /> 146 <field field_ref="tauy_sic" level="10" /> 147 <field field_ref="pourc_ter" level="10" /> 148 <field field_ref="pourc_lic" level="10" /> 149 <field field_ref="pourc_oce" level="10" /> 150 <field field_ref="pourc_sic" level="10" /> 151 <field field_ref="fract_ter" level="10" /> 152 <field field_ref="fract_lic" level="10" /> 153 <field field_ref="fract_oce" level="10" /> 154 <field field_ref="fract_sic" level="10" /> 155 <field field_ref="tsol_ter" level="10" /> 156 <field field_ref="tsol_lic" level="10" /> 157 <field field_ref="tsol_oce" level="10" /> 158 <field field_ref="tsol_sic" level="10" /> 159 <field field_ref="evappot_ter" level="10" /> 160 <field field_ref="evappot_lic" level="10" /> 161 <field field_ref="evappot_oce" level="10" /> 162 <field field_ref="evappot_sic" level="10" /> 163 <field field_ref="sens_ter" level="7" /> 164 <field field_ref="sens_lic" level="7" /> 165 <field field_ref="sens_oce" level="7" /> 166 <field field_ref="sens_sic" level="7" /> 167 <field field_ref="lat_ter" level="7" /> 168 <field field_ref="lat_lic" level="7" /> 169 <field field_ref="lat_oce" level="7" /> 170 <field field_ref="lat_sic" level="7" /> 171 <field field_ref="flw_ter" level="10" /> 172 <field field_ref="flw_lic" level="10" /> 173 <field field_ref="flw_oce" level="10" /> 174 <field field_ref="flw_sic" level="10" /> 175 <field field_ref="fsw_ter" level="10" /> 176 <field field_ref="fsw_lic" level="10" /> 177 <field field_ref="fsw_oce" level="10" /> 178 <field field_ref="fsw_sic" level="10" /> 179 <field field_ref="wbils_ter" level="10" /> 180 <field field_ref="wbils_lic" level="10" /> 181 <field field_ref="wbils_oce" level="10" /> 182 <field field_ref="wbils_sic" level="10" /> 183 <field field_ref="wbilo_ter" level="10" /> 184 <field field_ref="wbilo_lic" level="10" /> 185 <field field_ref="wbilo_oce" level="10" /> 186 <field field_ref="wbilo_sic" level="10" /> 187 <field field_ref="wevap_ter" level="10" /> 188 <field field_ref="wevap_lic" level="10" /> 189 <field field_ref="wevap_oce" level="10" /> 190 <field field_ref="wevap_sic" level="10" /> 191 <field field_ref="wrain_ter" level="10" /> 192 <field field_ref="wrain_lic" level="10" /> 193 <field field_ref="wrain_oce" level="10" /> 194 <field field_ref="wrain_sic" level="10" /> 195 <field field_ref="wsnow_ter" level="10" /> 196 <field field_ref="wsnow_lic" level="10" /> 197 <field field_ref="wsnow_oce" level="10" /> 198 <field field_ref="wsnow_sic" level="10" /> 199 <field field_ref="cdrm" level="10" /> 200 <field field_ref="cdrh" level="7" /> 201 <field field_ref="cldl" level="10" /> 202 <field field_ref="cldm" level="10" /> 203 <field field_ref="cldh" level="10" /> 204 <field field_ref="cldt" level="5" /> 205 <field field_ref="cldq" level="10" /> 206 <field field_ref="lwp" level="10" /> 207 <field field_ref="iwp" level="10" /> 208 <field field_ref="ue" level="10" /> 209 <field field_ref="ve" level="10" /> 210 <field field_ref="uq" level="10" /> 211 <field field_ref="vq" level="10" /> 212 <field field_ref="cape" level="10" /> 213 <field field_ref="pbase" level="10" /> 214 <field field_ref="ptop" level="10" /> 215 <field field_ref="fbase" level="10" /> 216 <field field_ref="plcl" level="10" /> 217 <field field_ref="plfc" level="10" /> 218 <field field_ref="wbeff" level="10" /> 219 <field field_ref="convoccur" level="10" /> 220 <field field_ref="prw" level="5" /> 221 <field field_ref="prlw" level="10" /> 222 <field field_ref="prsw" level="10" /> 223 <field field_ref="prbsw" level="5" /> 224 <field field_ref="s_pblh" level="10" /> 225 <field field_ref="s_pblt" level="10" /> 226 <field field_ref="s_lcl" level="10" /> 227 <field field_ref="s_therm" level="10" /> 228 <field field_ref="slab_bils_oce" level="10" /> 229 <field field_ref="ale_bl" level="10" /> 230 <field field_ref="alp_bl" level="10" /> 231 <field field_ref="ale_wk" level="10" /> 232 <field field_ref="alp_wk" level="10" /> 233 <field field_ref="ale" level="10" /> 234 <field field_ref="alp" level="10" /> 235 <field field_ref="cin" level="10" /> 236 <field field_ref="wape" level="10" /> 237 <field field_ref="n2" level="6" /> 238 <field field_ref="s2" level="6" /> 239 <field field_ref="strig" level="6" /> 240 <field field_ref="zcong" level="6" /> 241 <field field_ref="zlcl_th" level="6" /> 242 <field field_ref="proba_notrig" level="6" /> 243 <field field_ref="random_notrig" level="6" /> 244 <field field_ref="cv_gen" level="6" /> 245 <field field_ref="ale_bl_stat" level="6" /> 246 <field field_ref="ale_bl_trig" level="6" /> 247 <field field_ref="alp_bl_det" level="10" /> 248 <field field_ref="alp_bl_fluct_m" level="10" /> 249 <field field_ref="alp_bl_fluct_tke" level="10" /> 250 <field field_ref="alp_bl_conv" level="10" /> 251 <field field_ref="alp_bl_stat" level="10" /> 252 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 253 <field_group operation="average" freq_op="1ts" detect_missing_value=".true."> 254 <field field_ref="u850" level="10" /> 255 <field field_ref="u700" level="10" /> 256 <field field_ref="u500" level="5" /> 257 <field field_ref="u200" level="10" /> 258 <field field_ref="u100" level="10" /> 259 <field field_ref="u50" level="10" /> 260 <field field_ref="u10" level="10" /> 261 <field field_ref="v850" level="10" /> 262 <field field_ref="v700" level="10" /> 263 <field field_ref="v500" level="5" /> 264 <field field_ref="v200" level="10" /> 265 <field field_ref="v100" level="10" /> 266 <field field_ref="v50" level="10" /> 267 <field field_ref="v10" level="10" /> 268 <field field_ref="w850" level="10" /> 269 <field field_ref="w700" level="10" /> 270 <field field_ref="w500" level="10" /> 271 <field field_ref="w200" level="10" /> 272 <field field_ref="w100" level="10" /> 273 <field field_ref="w50" level="10" /> 274 <field field_ref="w10" level="10" /> 275 <field field_ref="t850" level="10" /> 276 <field field_ref="t700" level="10" /> 277 <field field_ref="t500" level="5" /> 278 <field field_ref="t200" level="10" /> 279 <field field_ref="t100" level="10" /> 280 <field field_ref="t50" level="10" /> 281 <field field_ref="t10" level="10" /> 282 <field field_ref="q850" level="10" /> 283 <field field_ref="q700" level="10" /> 284 <field field_ref="q500" level="10" /> 285 <field field_ref="q200" level="10" /> 286 <field field_ref="q100" level="10" /> 287 <field field_ref="q50" level="10" /> 288 <field field_ref="q10" level="10" /> 289 <field field_ref="z850" level="10" /> 290 <field field_ref="z700" level="10" /> 291 <field field_ref="z500" level="5" /> 292 <field field_ref="z200" level="10" /> 293 <field field_ref="z100" level="10" /> 294 <field field_ref="z50" level="10" /> 295 <field field_ref="z10" level="10" /> 296 </field_group> 297 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 298 <field field_ref="t_oce_sic" level="10" /> 299 <field field_ref="weakinv" level="10" /> 300 <field field_ref="dthmin" level="10" /> 301 <field field_ref="u10_ter" level="10" /> 302 <field field_ref="u10_lic" level="10" /> 303 <field field_ref="u10_oce" level="10" /> 304 <field field_ref="u10_sic" level="10" /> 305 <field field_ref="v10_ter" level="10" /> 306 <field field_ref="v10_lic" level="10" /> 307 <field field_ref="v10_oce" level="10" /> 308 <field field_ref="v10_sic" level="10" /> 309 <field field_ref="rh2m" level="5" /> 310 <field field_ref="rh2m_min" level="10" operation="minimum" /> 311 <field field_ref="rh2m_max" level="10" operation="maximum" /> 312 <field field_ref="qsat2m" level="10" /> 313 <field field_ref="tpot" level="10" /> 314 <field field_ref="tpote" level="10" /> 315 <field field_ref="SWnetOR" level="10" /> 316 <field field_ref="SWdownOR" level="5" /> 317 <field field_ref="LWdownOR" level="5" /> 318 <field field_ref="snowl" level="10" /> 319 <field field_ref="cape_max" level="10" operation="maximum" /> 320 <field field_ref="solldown" level="10" /> 321 <field field_ref="dtsvdfo" level="10" /> 322 <field field_ref="dtsvdft" level="10" /> 323 <field field_ref="dtsvdfg" level="10" /> 324 <field field_ref="dtsvdfi" level="10" /> 325 <field field_ref="rugs" level="10" /> 326 <field field_ref="z0m" level="10" /> 327 <field field_ref="z0h" level="10" /> 328 <field field_ref="topswad" level="10" /> 329 <field field_ref="topswad0" level="10" /> 330 <field field_ref="topswai" level="10" /> 331 <field field_ref="solswad" level="10" /> 332 <field field_ref="solswad0" level="10" /> 333 <field field_ref="solswai" level="10" /> 334 <field field_ref="od550_ASBCM" level="10" /> 335 <field field_ref="od550_ASPOMM" level="10" /> 336 <field field_ref="od550_ASSO4M" level="10" /> 337 <field field_ref="od550_CSSO4M" level="10" /> 338 <field field_ref="od550_SSSSM" level="10" /> 339 <field field_ref="od550_ASSSM" level="10" /> 340 <field field_ref="od550_CSSSM" level="10" /> 341 <field field_ref="od550_CIDUSTM" level="10" /> 342 <field field_ref="od550_AIBCM" level="10" /> 343 <field field_ref="od550_AIPOMM" level="10" /> 344 <field field_ref="od550_ASNO3M" level="10" /> 345 <field field_ref="od550_CSNO3M" level="10" /> 346 <field field_ref="od550_CINO3M" level="10" /> 347 <field field_ref="od550_STRAT" level="10" /> 348 <field field_ref="od550aer" level="10" /> 349 <field field_ref="od865aer" level="10" /> 350 <field field_ref="abs550aer" level="10" /> 351 <field field_ref="od550lt1aer" level="10" /> 352 <field field_ref="sconcso4" level="10" /> 353 <field field_ref="sconcoa" level="10" /> 354 <field field_ref="sconcbc" level="10" /> 355 <field field_ref="sconcss" level="10" /> 356 <field field_ref="sconcdust" level="10" /> 357 <field field_ref="sconcno3" level="10" /> 358 <field field_ref="loadso4" level="10" /> 359 <field field_ref="loadoa" level="10" /> 360 <field field_ref="loadbc" level="10" /> 361 <field field_ref="loadss" level="10" /> 362 <field field_ref="loaddust" level="10" /> 363 <field field_ref="loadno3" level="10" /> 364 <field field_ref="swtoaas_nat" level="10" /> 365 <field field_ref="swsrfas_nat" level="10" /> 366 <field field_ref="swtoacs_nat" level="10" /> 367 <field field_ref="swsrfcs_nat" level="10" /> 368 <field field_ref="swtoaas_ant" level="10" /> 369 <field field_ref="swsrfas_ant" level="10" /> 370 <field field_ref="swtoacs_ant" level="10" /> 371 <field field_ref="swsrfcs_ant" level="10" /> 372 <field field_ref="swtoacf_nat" level="10" /> 373 <field field_ref="swsrfcf_nat" level="10" /> 374 <field field_ref="swtoacf_ant" level="10" /> 375 <field field_ref="swsrfcf_ant" level="10" /> 376 <field field_ref="swtoacf_zero" level="10" /> 377 <field field_ref="swsrfcf_zero" level="10" /> 378 <field field_ref="cldncl" level="10" /> 379 <field field_ref="reffclwtop" level="10" /> 380 <field field_ref="cldnvi" level="10" /> 381 <field field_ref="lcc" level="10" /> 382 <field field_ref="wvapp" level="10" /> 383 <field field_ref="albe_ter" level="7" /> 384 <field field_ref="albe_lic" level="7" /> 385 <field field_ref="albe_oce" level="7" /> 386 <field field_ref="albe_sic" level="7" /> 387 <field field_ref="ages_ter" level="10" /> 388 <field field_ref="ages_lic" level="10" /> 389 <field field_ref="ages_oce" level="10" /> 390 <field field_ref="ages_sic" level="10" /> 391 <field field_ref="snow_ter" level="10" /> 392 <field field_ref="snow_lic" level="10" /> 393 <field field_ref="snow_oce" level="10" /> 394 <field field_ref="snow_sic" level="10" /> 395 <field field_ref="rugs_ter" level="10" /> 396 <field field_ref="rugs_lic" level="10" /> 397 <field field_ref="rugs_oce" level="10" /> 398 <field field_ref="rugs_sic" level="10" /> 399 <field field_ref="z0m_ter" level="10" /> 400 <field field_ref="z0m_lic" level="10" /> 401 <field field_ref="z0m_oce" level="10" /> 402 <field field_ref="z0m_sic" level="10" /> 403 <field field_ref="z0h_ter" level="10" /> 404 <field field_ref="z0h_lic" level="10" /> 405 <field field_ref="z0h_oce" level="10" /> 406 <field field_ref="z0h_sic" level="10" /> 407 <field field_ref="alb1" level="10" /> 408 <field field_ref="alb2" level="10" /> 409 <field field_ref="ftime_con" level="10" operation="instant" /> 410 <field field_ref="wake_h" level="10" /> 411 <field field_ref="wake_dens" level="10" /> 412 <field field_ref="wake_s" level="10" /> 413 <field field_ref="epmax" level="10" /> 414 <field field_ref="plulth" level="10" /> 415 <field field_ref="plulst" level="10" /> 416 <field field_ref="lmaxth" level="10" /> 417 <field field_ref="ftime_th" level="10" /> 418 <field field_ref="ftime_deepcv" level="10" /> 419 <field field_ref="f0_th" level="10" /> 420 <field field_ref="zmax_th" level="10" /> 421 <field field_ref="rsut4co2" level="10" /> 422 <field field_ref="rlut4co2" level="10" /> 423 <field field_ref="rsutcs4co2" level="10" /> 424 <field field_ref="rlutcs4co2" level="10" /> 425 <field field_ref="dqphy2d" level="10" /> 426 <field field_ref="dqlphy2d" level="10" /> 427 <field field_ref="dqsphy2d" level="10" /> 428 <field field_ref="dqdyn2d" level="10" /> 429 <field field_ref="dqldyn2d" level="10" /> 430 <field field_ref="dqsdyn2d" level="10" /> 431 <field field_ref="dqcon2d" level="10" /> 432 <field field_ref="dqwak2d" level="10" /> 433 <field field_ref="dqlsc2d" level="10" /> 434 <field field_ref="dqvdf2d" level="10" /> 435 <field field_ref="dqeva2d" level="10" /> 436 <field field_ref="dqlscth2d" level="10" /> 437 <field field_ref="dqlscst2d" level="10" /> 438 <field field_ref="dqthe2d" level="10" /> 439 <field field_ref="dqajs2d" level="10" /> 440 <field field_ref="flx_co2_ocean" level="10" /> <!-- Added PC --> 441 <field field_ref="flx_co2_land" level="10" /> <!-- Added PC --> 442 <field field_ref="flx_co2_ocean_cor" level="10" /> <!-- Added PC --> 443 <field field_ref="flx_co2_land_cor" level="10" /> <!-- Added PC --> 444 <field field_ref="flx_co2_ff" level="10" /> <!-- Added PC --> 445 <field field_ref="flx_co2_bb" level="10" /> <!-- Added PC --> 446 446 447 447 </field_group> 448 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 449 <field_group operation="instant" freq_op="1ts" > 450 <field field_ref="alt_tropo" level="10" /> 451 <field field_ref="map_prop_hc" level="10" /> 452 <field field_ref="map_prop_hist" level="10" /> 453 <field field_ref="map_emis_hc" level="10" /> 454 <field field_ref="map_iwp_hc" level="10" /> 455 <field field_ref="map_deltaz_hc" level="10" /> 456 <field field_ref="map_pcld_hc" level="10" /> 457 <field field_ref="map_tcld_hc" level="10" /> 458 <field field_ref="map_emis_hist" level="10" /> 459 <field field_ref="map_iwp_hist" level="10" /> 460 <field field_ref="map_deltaz_hist" level="10" /> 461 <field field_ref="map_rad_hist" level="10" /> 462 <field field_ref="map_emis_Cb" level="10" /> 463 <field field_ref="map_pcld_Cb" level="10" /> 464 <field field_ref="map_tcld_Cb" level="10" /> 465 <field field_ref="map_emis_Anv" level="10" /> 466 <field field_ref="map_pcld_Anv" level="10" /> 467 <field field_ref="map_tcld_Anv" level="10" /> 468 <field field_ref="map_emis_ThCi" level="10" /> 469 <field field_ref="map_pcld_ThCi" level="10" /> 470 <field field_ref="map_tcld_ThCi" level="10" /> 471 <field field_ref="map_ntot" level="10" /> 472 <field field_ref="map_hc" level="10" /> 473 <field field_ref="map_hist" level="10" /> 474 <field field_ref="map_Cb" level="10" /> 475 <field field_ref="map_ThCi" level="10" /> 476 <field field_ref="map_Anv" level="10" /> 477 </field_group> 478 478 479 <!-- VARS 3D --> 480 <field_group operation="average" grid_ref="grid_out_presnivs"> 481 <field field_ref="CO2" level="10" /> <!-- Added PC --> 482 <field field_ref="tke" level="10" /> 483 <field field_ref="tke_ter" level="10" /> 484 <field field_ref="tke_lic" level="10" /> 485 <field field_ref="tke_oce" level="10" /> 486 <field field_ref="tke_sic" level="10" /> 487 <field field_ref="tke_max_ter" level="10" operation="maximum" /> 488 <field field_ref="tke_max_lic" level="10" operation="maximum" /> 489 <field field_ref="tke_max_oce" level="10" operation="maximum" /> 490 <field field_ref="tke_max_sic" level="10" operation="maximum" /> 491 <field field_ref="l_mix_ter" level="10" /> 492 <field field_ref="l_mix_lic" level="10" /> 493 <field field_ref="l_mix_oce" level="10" /> 494 <field field_ref="l_mix_sic" level="10" /> 495 <field field_ref="l_mixmin_ter" level="10" /> 496 <field field_ref="l_mixmin_lic" level="10" /> 497 <field field_ref="l_mixmin_oce" level="10" /> 498 <field field_ref="l_mixmin_sic" level="10" /> 499 <field field_ref="treedrg_ter" level="10" /> 500 <field field_ref="treedrg_lic" level="10" /> 501 <field field_ref="treedrg_oce" level="10" /> 502 <field field_ref="treedrg_sic" level="10" /> 503 <field field_ref="cldtau" level="10" /> 504 <field field_ref="cldemi" level="10" /> 505 <field field_ref="tke_max" level="10" operation="maximum" /> 506 <field field_ref="concso4" level="10" /> 507 <field field_ref="concoa" level="10" /> 508 <field field_ref="concbc" level="10" /> 509 <field field_ref="concss" level="10" /> 510 <field field_ref="concdust" level="10" /> 511 <field field_ref="concno3" level="10" /> 512 <field field_ref="ec550aer" level="10" /> 513 <field field_ref="lwcon" level="5" /> 514 <field field_ref="iwcon" level="5" /> 515 <field field_ref="temp" level="4" /> 516 <field field_ref="theta" level="10" /> 517 <field field_ref="ovap" level="5" /> 518 <field field_ref="ovapinit" level="10" /> 519 <field field_ref="oliq" level="5" /> 520 <field field_ref="ocond" level="5" /> 521 <field field_ref="qbs" level="5" /> 522 <field field_ref="geop" level="10" /> 523 <field field_ref="geoph" level="10" /> 524 <field field_ref="vitu" level="5" /> 525 <field field_ref="vitv" level="5" /> 526 <field field_ref="vitw" level="6" /> 527 <field field_ref="pres" level="5" /> 528 <field field_ref="paprs" level="10" /> 529 <field field_ref="mass" level="10" /> 530 <field field_ref="zfull" level="5" /> 531 <field field_ref="zhalf" level="10" /> 532 <field field_ref="rneb" level="10" /> 533 <field field_ref="rnebcon" level="10" /> 534 <field field_ref="rnebls" level="10" /> 535 <field field_ref="rneblsvol" level="10" /> 536 <field field_ref="rhum" level="10" /> 537 <field field_ref="rhi" level="5" /> 538 <field field_ref="ozone" level="10" /> 539 <field field_ref="upwd" level="10" /> 540 <field field_ref="ep" level="10" /> 541 <field field_ref="duphy" level="10" /> 542 <field field_ref="dtphy" level="5" /> 543 <field field_ref="dqphy" level="5" /> 544 <field field_ref="dqlphy" level="10" /> 545 <field field_ref="dqsphy" level="10" /> 546 <field field_ref="pr_con_l" level="10" /> 547 <field field_ref="pr_con_i" level="10" /> 548 <field field_ref="pr_lsc_l" level="10" /> 549 <field field_ref="pr_lsc_i" level="10" /> 550 <field field_ref="pr_bs" level="5" /> 551 <field field_ref="re" level="10" /> 552 <field field_ref="fl" level="10" /> 553 <field field_ref="scdnc" level="10" /> 554 <field field_ref="reffclws" level="10" /> 479 <!-- VARS 3D --> 480 <field_group operation="average" grid_ref="grid_out_presnivs"> 481 <field field_ref="flightdist" level="1" operation="instant" /> 482 <field field_ref="CO2" level="10" /> <!-- Added PC --> 483 <field field_ref="tke" level="10" /> 484 <field field_ref="tke_ter" level="10" /> 485 <field field_ref="tke_lic" level="10" /> 486 <field field_ref="tke_oce" level="10" /> 487 <field field_ref="tke_sic" level="10" /> 488 <field field_ref="tke_max_ter" level="10" operation="maximum" /> 489 <field field_ref="tke_max_lic" level="10" operation="maximum" /> 490 <field field_ref="tke_max_oce" level="10" operation="maximum" /> 491 <field field_ref="tke_max_sic" level="10" operation="maximum" /> 492 <field field_ref="l_mix_ter" level="10" /> 493 <field field_ref="l_mix_lic" level="10" /> 494 <field field_ref="l_mix_oce" level="10" /> 495 <field field_ref="l_mix_sic" level="10" /> 496 <field field_ref="l_mixmin_ter" level="10" /> 497 <field field_ref="l_mixmin_lic" level="10" /> 498 <field field_ref="l_mixmin_oce" level="10" /> 499 <field field_ref="l_mixmin_sic" level="10" /> 500 <field field_ref="treedrg_ter" level="10" /> 501 <field field_ref="treedrg_lic" level="10" /> 502 <field field_ref="treedrg_oce" level="10" /> 503 <field field_ref="treedrg_sic" level="10" /> 504 <field field_ref="cldtau" level="10" /> 505 <field field_ref="cldemi" level="10" /> 506 <field field_ref="tke_max" level="10" operation="maximum" /> 507 <field field_ref="concso4" level="10" /> 508 <field field_ref="concoa" level="10" /> 509 <field field_ref="concbc" level="10" /> 510 <field field_ref="concss" level="10" /> 511 <field field_ref="concdust" level="10" /> 512 <field field_ref="concno3" level="10" /> 513 <field field_ref="ec550aer" level="10" /> 514 <field field_ref="lwcon" level="5" /> 515 <field field_ref="iwcon" level="5" /> 516 <field field_ref="temp" level="4" /> 517 <field field_ref="theta" level="10" /> 518 <field field_ref="ovap" level="5" /> 519 <field field_ref="ovapinit" level="10" /> 520 <field field_ref="oliq" level="5" /> 521 <field field_ref="ocond" level="5" /> 522 <field field_ref="qbs" level="5" /> 523 <field field_ref="geop" level="10" /> 524 <field field_ref="geoph" level="10" /> 525 <field field_ref="vitu" level="5" /> 526 <field field_ref="vitv" level="5" /> 527 <field field_ref="vitw" level="6" /> 528 <field field_ref="pres" level="5" /> 529 <field field_ref="paprs" level="10" /> 530 <field field_ref="mass" level="10" /> 531 <field field_ref="zfull" level="5" /> 532 <field field_ref="zhalf" level="10" /> 533 <field field_ref="rneb" level="10" /> 534 <field field_ref="rnebcon" level="10" /> 535 <field field_ref="rnebls" level="10" /> 536 <field field_ref="rneblsvol" level="10" /> 537 <field field_ref="rhum" level="10" /> 538 <field field_ref="rhi" level="5" /> 539 <field field_ref="ozone" level="10" /> 540 <field field_ref="upwd" level="10" /> 541 <field field_ref="ep" level="10" /> 542 <field field_ref="duphy" level="10" /> 543 <field field_ref="dtphy" level="5" /> 544 <field field_ref="dqphy" level="5" /> 545 <field field_ref="dqlphy" level="10" /> 546 <field field_ref="dqsphy" level="10" /> 547 <field field_ref="pr_con_l" level="10" /> 548 <field field_ref="pr_con_i" level="10" /> 549 <field field_ref="pr_lsc_l" level="10" /> 550 <field field_ref="pr_lsc_i" level="10" /> 551 <field field_ref="pr_bs" level="5" /> 552 <field field_ref="re" level="10" /> 553 <field field_ref="fl" level="10" /> 554 <field field_ref="scdnc" level="10" /> 555 <field field_ref="reffclws" level="10" /> 555 556 <field field_ref="reffclwc" level="10" /> 556 557 <field field_ref="lcc3d" level="10" /> -
LMDZ6/branches/contrails/DefLists/file_def_histmth_lmdz.xml
r5343 r5573 553 553 <!-- VARS 3D --> 554 554 <field_group operation="average" grid_ref="grid_out_presnivs"> 555 <field field_ref="flightdist" level="1" /> 556 <field field_ref="cfseri" level="1" /> 557 <field field_ref="qissr" level="1" /> 558 <field field_ref="qcld" level="1" /> 559 <field field_ref="subfra" level="1" /> 560 <field field_ref="issrfra" level="1" /> 561 <field field_ref="dqvcmix" level="1" /> 562 <field field_ref="dqimix" level="1" /> 563 <field field_ref="dcfmix" level="1" /> 564 <field field_ref="fallice_velo" level="1" /> 555 565 <field field_ref="CO2" level="5" /> <!-- Added PC --> 556 566 <field field_ref="dCO2_vdf" level="5" /> <!-- Added PC --> -
LMDZ6/branches/contrails/libf/phylmd/lmdz_aviation.f90
r5551 r5573 622 622 END FUNCTION contrail_cross_section_onera 623 623 624 SUBROUTINE read_aviation_emissions(klon, klev, flight_dist , flight_h2o)624 SUBROUTINE read_aviation_emissions(klon, klev, flight_dist_read, flight_h2o_read, aviation_lev, nleva) 625 625 ! This subroutine allows to read the traffic density data read in the file aviation.nc 626 626 ! This file is defined in ./COMP/lmdz.card … … 635 635 636 636 INTEGER, INTENT(IN) :: klon, klev ! number of horizontal grid points and vertical levels 637 REAL, DIMENSION(klon,klev), INTENT(OUT) :: flight_dist ! Aviation distance flown concentration [m/s/m3] 638 REAL, DIMENSION(klon,klev), INTENT(OUT) :: flight_h2o ! Aviation emitted H2O [kgH2O/s/m3] 637 INTEGER, INTENT(out) :: nleva ! Size of the vertical axis in the file 638 !REAL, DIMENSION(klon,klev), INTENT(OUT) :: flight_dist ! Aviation distance flown within the mesh [m/s/mesh] 639 !REAL, DIMENSION(klon,klev), INTENT(OUT) :: flight_h2o ! Aviation H2O emitted within the mesh [kgH2O/s/mesh] 640 REAL, ALLOCATABLE, INTENT(OUT) :: flight_dist_read(:,:,:) ! Aviation distance flown within the mesh [m/s/mesh] 641 REAL, ALLOCATABLE, INTENT(OUT) :: flight_h2o_read(:,:,:) ! Aviation H2O emitted within the mesh [kgH2O/s/mesh] 642 REAL, ALLOCATABLE, INTENT(OUT) :: aviation_lev(:) ! Pressure in the middle of the layers [Pa] 639 643 640 644 !---------------------------------------------------- 641 645 ! Local variable 642 646 !---------------------------------------------------- 643 REAL, DIMENSION(klon_mpi,klev,1) :: flight_dist_mpi 647 !REAL, DIMENSION(klon_mpi,klev,1) :: flight_dist_mpi 648 REAL, ALLOCATABLE :: flight_dist_mpi(:,:,:) 649 INTEGER :: ierr 650 651 ! Get number of vertical levels and level values 652 CALL xios_get_axis_attr( "aviation_lev", n_glo=nleva ) 653 654 ! Allocation of arrays 655 !$OMP MASTER 656 ALLOCATE(aviation_lev(nleva), STAT=ierr) 657 IF (ierr /= 0) CALL abort_physic('read_aviation_emissions', 'problem to allocate aviation_lev',1) 658 ALLOCATE(flight_dist_read(klon, nleva,1), STAT=ierr) 659 IF (ierr /= 0) CALL abort_physic('read_aviation_emissions', 'problem to allocate flight_dist',1) 660 ALLOCATE(flight_h2o_read(klon, nleva,1), STAT=ierr) 661 IF (ierr /= 0) CALL abort_physic('read_aviation_emissions', 'problem to allocate flight_h2o',1) 662 ALLOCATE(flight_dist_mpi(klon_mpi, nleva,1), STAT=ierr) 663 IF (ierr /= 0) CALL abort_physic('read_aviation_emissions', 'problem to allocate flight_dist_mpi',1) 664 !$OMP END MASTER 665 666 !$OMP BARRIER ! Ensure all threads wait until the arrays are allocated 644 667 645 668 !--Initialisation 646 flight_dist(:,:) = 0. 647 flight_h2o(:,:) = 0. 669 aviation_lev(:) = 0. 670 flight_dist_read(:,:,1) = 0. 671 flight_h2o_read(:,:,1) = 0. 672 673 ! Get number of vertical levels and level values 674 CALL xios_get_axis_attr( "aviation_lev", value=aviation_lev(:)) 648 675 649 676 ! Read the data from the file … … 651 678 IF (is_omp_master) CALL xios_recv_field("KMFLOWN_interp", flight_dist_mpi(:,:,1)) 652 679 653 654 655 CALL scatter_omp(flight_dist_mpi(:,:,1), flight_dist)680 ! Propagate to other OMP threads: flight_dist_mpi(klon_mpi,klev) to flight_dist(klon,klev) 681 ! (klon_mpi,klon) = (200,50) avec 80 MPI, 4 OMP, nbp40 682 CALL scatter_omp(flight_dist_mpi(:,:,1), flight_dist_read(:,:,1)) 656 683 657 684 END SUBROUTINE read_aviation_emissions 658 685 686 SUBROUTINE vertical_interpolation_aviation(klon, klev, paprs, pplay, flight_dist_read, & 687 flight_h2o_read, aviation_lev, nleva, flight_dist, & 688 flight_h2o) 689 ! This subroutine performs the vertical interpolation from the read data in aviation.nc 690 ! where there are nleva vertical levels described in aviation_lev to the klev levels or 691 ! the model. 692 ! flight_dist_read(klon,nleva) -> flight_dist(klon, klev) 693 ! flight_h2o_read(klon,nleva) -> flight_h2o(klon, klev) 694 USE print_control_mod, ONLY: lunout 695 IMPLICIT NONE 696 697 INTEGER, INTENT(IN) :: klon, klev ! number of horizontal grid points and vertical levels 698 REAL, INTENT(IN) :: paprs(klon, klev+1) ! inter-layer pressure [Pa] 699 REAL, INTENT(IN) :: pplay(klon, klev) ! mid-layer pressure [Pa] 700 INTEGER, INTENT(IN) :: nleva ! Size of the vertical axis in the file 701 REAL, INTENT(OUT) :: flight_dist(klon,klev,1) ! Aviation distance flown within the mesh [m/s/mesh] 702 REAL, INTENT(OUT) :: flight_h2o(klon,klev,1) ! Aviation H2O emitted within the mesh [kgH2O/s/mesh] 703 REAL, INTENT(IN) :: flight_dist_read(klon,nleva,1) ! Aviation distance flown within the mesh in file [m/s/mesh] 704 REAL, INTENT(IN) :: flight_h2o_read(klon,nleva,1) ! Aviation H2O emitted within the mesh in file [kgH2O/s/mesh] 705 REAL, INTENT(IN) :: aviation_lev(nleva) ! Pressure in the middle of the layers [Pa] 706 707 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 708 ! Local variable 709 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 710 REAL :: aviation_interface(1:nleva+1) ! Pressure of aviation file interfaces [ Pa ] 711 INTEGER :: k, kori ! Loop index for vertical layers 712 INTEGER :: i ! Loop index for horizontal grid 713 REAL :: zfrac ! Fraction of layer kori in layer k 714 REAL :: width_read_layer(1:nleva) ! width of a given layer [ Pa ] 715 716 ! Initialisation 717 flight_dist(:,:,1) = 0. 718 flight_h2o(:,:,1) = 0. 719 720 ! Compute the array with the vertical interface 721 ! It starts at 1 and has length nleva + 1 722 ! Note that aviation_lev has nleva and gives the altitude in the middle of the layers 723 ! Surface pressure in standard atmosphere model [ Pa ] 724 aviation_interface(1) = 101325. 725 DO kori=2, nleva 726 aviation_interface(kori) = (aviation_lev(kori-1)+aviation_lev(kori))/2.0 ! [ Pa ] 727 ENDDO 728 ! Last interface - we assume the same spacing as the very last one 729 aviation_interface(nleva+1) = aviation_interface(nleva) - (aviation_lev(nleva-1) - aviation_lev(nleva)) 730 731 ! Vertical width of each layer of the read file 732 ! It is positive 733 DO kori=1, nleva 734 width_read_layer(kori) = aviation_interface(kori) - aviation_interface(kori+1) 735 ENDDO 736 737 ! Vertical reprojection 738 ! The loop over klon is induced since it is done by MPI threads 739 ! zfrac is the fraction of layer kori (read file) included in layer k (model) 740 DO i=1,klon 741 DO k=1, klev 742 DO kori=1,nleva 743 ! Which of the lower interfaces is the highest (<=> the lowest pressure) ? 744 zfrac = min(paprs(i,k), aviation_interface(kori)) 745 ! Which of the upper interfaces is the lowest (<=> the greatest pressure) ? 746 zfrac = zfrac - max(paprs(i,k+1), aviation_interface(kori+1)) 747 ! If zfrac is negative, the layers are not overlapping 748 ! Otherwise, we get the fraction of layer kori that overlap with layer k 749 ! after normalisation to the total kori layer width 750 zfrac = max(0.0, zfrac) / width_read_layer(kori) 751 752 ! Vertical reprojection for each desired array 753 flight_dist(i,k,1) = flight_dist(i,k,1) + zfrac * flight_dist_read(i,kori,1) 754 flight_h2o(i,k,1) = flight_h2o(i,k,1) + zfrac * flight_h2o(i,kori,1) 755 ENDDO 756 ENDDO 757 ENDDO 758 759 END SUBROUTINE vertical_interpolation_aviation 760 659 761 END MODULE lmdz_aviation -
LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp.f90
r5551 r5573 24 24 dqi_adj, dqi_sub, dqi_con, dqi_mix, dqvc_adj, & 25 25 dqvc_sub, dqvc_con, dqvc_mix, qsatl, qsati, & 26 rcont_seri, flight_dist, flight_h2o, contfra, & 27 Tcritcont, qcritcont, potcontfraP, potcontfraNP, & 28 dcontfra_cir, dcf_avi, dqi_avi, dqvc_avi, & 29 cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, & 26 rcont_seri, flight_dist, flight_h2o, & 27 flight_dist_read, flight_h2o_read, & 28 aviation_lev, nleva, contfra, Tcritcont, qcritcont,& 29 potcontfraP, potcontfraNP, dcontfra_cir, dcf_avi, & 30 dqi_avi, dqvc_avi, cloudth_sth,cloudth_senv, & 31 cloudth_sigmath,cloudth_sigmaenv, & 30 32 qraindiag, qsnowdiag, dqreva, dqssub, dqrauto, & 31 33 dqrcol, dqrmelt, dqrfreez, dqsauto, dqsagg, dqsrim,& … … 122 124 USE lmdz_lscp_ini, ONLY : ok_plane_contrail 123 125 126 ! aviation module 127 USE mod_phys_lmdz_para, ONLY : is_omp_master 128 USE lmdz_aviation, ONLY : vertical_interpolation_aviation 129 124 130 IMPLICIT NONE 125 131 … … 175 181 !-------------------------------------------------- 176 182 REAL, DIMENSION(klon,klev), INTENT(INOUT):: rcont_seri ! ratio of contrails fraction to total cloud fraction [-] 177 REAL, DIMENSION(klon,klev), INTENT(IN) :: flight_dist ! aviation distance flown concentration [m/s/m3] 178 REAL, DIMENSION(klon,klev), INTENT(IN) :: flight_h2o ! aviation emitted H2O concentration [kgH2O/s/m3] 183 REAL, DIMENSION(klon,klev,1), INTENT(OUT) :: flight_dist ! aviation distance flown within the mesh [m/s/mesh] 184 REAL, DIMENSION(klon,klev,1), INTENT(OUT) :: flight_h2o ! aviation H2O emitted within the mesh [kgH2O/s/mesh] 185 ! Read from file - nleva vertical levels - define in state_var 186 INTEGER, INTENT(IN) :: nleva 187 REAL, DIMENSION(nleva), INTENT(IN) :: aviation_lev ! vertical levels [km] 188 REAL, DIMENSION(klon,nleva,1), INTENT(IN) :: flight_dist_read ! aviation distance flown within the mesh [m/s/mesh] 189 REAL, DIMENSION(klon,nleva,1), INTENT(IN) :: flight_h2o_read ! aviation H2O emitted within the mesh [kgH2O/s/mesh] 179 190 180 191 ! OUTPUT variables … … 452 463 !c_iso: variable initialisation for iso 453 464 465 IF ( ok_plane_contrail ) THEN 466 ! Vertical interpolation is done at each physical timestep 467 ! 468 IF (is_omp_master) CALL vertical_interpolation_aviation(klon, klev, paprs, pplay, flight_dist_read, & 469 flight_h2o_read, aviation_lev, nleva, flight_dist, flight_h2o) 470 ENDIF 454 471 455 472 !=============================================================================== … … 731 748 dqi_adj(:,k), dqi_sub(:,k), dqi_con(:,k), dqi_mix(:,k), & 732 749 dqvc_adj(:,k), dqvc_sub(:,k), dqvc_con(:,k), dqvc_mix(:,k), & 733 rcont_seri(:,k), flight_dist(:,k ), flight_h2o(:,k), contfra(:,k), &750 rcont_seri(:,k), flight_dist(:,k,1), flight_h2o(:,k,1), contfra(:,k), & 734 751 Tcritcont(:,k), qcritcont(:,k), potcontfraP(:,k), potcontfraNP(:,k), & 735 752 dcontfra_cir(:,k), dcf_avi(:,k), dqi_avi(:,k), dqvc_avi(:,k)) -
LMDZ6/branches/contrails/libf/phylmd/phys_output_ctrlout_mod.F90
r5551 r5573 2193 2193 'dqvcavi', 'Aviation cloudy water vapor tendency', 'kg/kg/s', (/ ('', i=1, 10) /)) 2194 2194 TYPE(ctrl_out), SAVE :: o_flight_dist = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), & 2195 'flightdist', 'Aviation flown distance concentration', 'm/s/m3', (/ ('', i=1, 10) /))2195 'flightdist', 'Aviation flown distance', 'm/s/m^3', (/ ('', i=1, 10) /)) 2196 2196 TYPE(ctrl_out), SAVE :: o_flight_h2o = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), & 2197 'flighth2o', 'Aviation emitted H2O concentration', 'kg H2O/s/m3', (/ ('', i=1, 10) /))2197 'flighth2o', 'Aviation H2O flight emission', 'kg H2O/s/m^3', (/ ('', i=1, 10) /)) 2198 2198 2199 2199 !!!!!!!!!!!!! Sorties niveaux standards de pression NMC -
LMDZ6/branches/contrails/libf/phylmd/phys_state_var_mod.F90
r5452 r5573 528 528 !$OMP THREADPRIVATE(delta_sal, ds_ns, dt_ns, delta_sst, dter, dser, dt_ds) 529 529 530 ! Arrays for the lecture of aviation files 531 ! The allocation is done in the read_aviation module 532 ! The size is (klon, nleva, 1) where 533 ! nleva is the size of the vertical axis (read from file) 534 ! flight_dist_read is the number of km per second 535 ! flight_h2o_read is the water content added to the air 536 ! aviation_lev is the value of the levels 537 REAL, SAVE, ALLOCATABLE :: flight_dist_read(:,:,:), flight_h2o_read(:,:,:) 538 REAL, SAVE, ALLOCATABLE :: aviation_lev(:) 539 !$OMP THREADPRIVATE(flight_dist_read, flight_h2o_read, aviation_lev) 540 INTEGER, SAVE :: nleva 541 !$OMP THREADPRIVATE(nleva) 542 543 530 544 CONTAINS 531 545 … … 942 956 DEALLOCATE(ratqs_inter_,sigma_qtherm) 943 957 958 ! DEALLOCATE aviation arrays 959 DEALLOCATE(flight_dist_read, flight_h2o_read, aviation_lev) 960 944 961 if (activate_ocean_skin >= 1) then 945 962 deALLOCATE(delta_sal, ds_ns, dt_ns, delta_sst, dter, dser) -
LMDZ6/branches/contrails/libf/phylmd/physiq_mod.F90
r5551 r5573 2124 2124 !--Read the aviation emissions 2125 2125 IF ( ok_plane_h2o .OR. ok_plane_contrail ) THEN 2126 CALL read_aviation_emissions(klon, klev, flight_dist, flight_h2o) 2126 CALL read_aviation_emissions(klon, klev, flight_dist_read, flight_h2o_read, & 2127 aviation_lev, nleva) 2127 2128 ENDIF 2128 2129 ! … … 3914 3915 dcf_sub, dcf_con, dcf_mix, dqi_adj, dqi_sub, dqi_con, dqi_mix, & 3915 3916 dqvc_adj, dqvc_sub, dqvc_con, dqvc_mix, qsatliq, qsatice, & 3916 rcont_seri, flight_dist, flight_h2o, contfra, &3917 Tcritcont, qcritcont, potcontfraP, potcontfraNP, &3918 dcontfra_cir, dcf_avi, dqi_avi, dqvc_avi, &3917 rcont_seri, flight_dist, flight_h2o, flight_dist_read, flight_h2o_read, & 3918 aviation_lev, nleva, contfra, Tcritcont, qcritcont, potcontfraP, & 3919 potcontfraNP, dcontfra_cir, dcf_avi, dqi_avi, dqvc_avi, & 3919 3920 cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, & 3920 3921 qraindiag, qsnowdiag, dqreva, dqssub, dqrauto, dqrcol, dqrmelt, &
Note: See TracChangeset
for help on using the changeset viewer.
