Changeset 5276 for LMDZ6/trunk/libf/phylmd
- Timestamp:
- Oct 25, 2024, 5:59:01 PM (9 months ago)
- Location:
- LMDZ6/trunk/libf/phylmd
- Files:
-
- 11 edited
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/cv30_routines.f90
r5275 r5276 5 5 6 6 SUBROUTINE cv30_param(nd, delt) 7 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 8 , clmci, eps, epsi, epsim1, ginv, hrd, grav 7 9 IMPLICIT NONE 8 10 … … 85 87 SUBROUTINE cv30_prelim(len, nd, ndp1, t, q, p, ph, lv, cpn, tv, gz, h, hm, & 86 88 th) 89 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 90 , clmci, eps, epsi, epsim1, ginv, hrd, grav 87 91 IMPLICIT NONE 88 92 … … 108 112 REAL cpx(len, nd) 109 113 110 include "cvthermo.h"111 114 include "cv30param.h" 112 115 … … 343 346 SUBROUTINE cv30_undilute1(len, nd, t, q, qs, gz, plcl, p, nk, icb, tp, tvp, & 344 347 clw, icbs) 348 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 349 , clmci, eps, epsi, epsim1, ginv, hrd, grav 345 350 IMPLICIT NONE 346 351 … … 357 362 ! - if icbs=icb, compute also tp(icb+1),tvp(icb+1) & clw(icb+1) 358 363 ! ---------------------------------------------------------------- 359 360 include "cvthermo.h"361 364 include "cv30param.h" 362 365 … … 839 842 ep, sigp, buoy) 840 843 ! epmax_cape: ajout arguments 844 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 845 , clmci, eps, epsi, epsim1, ginv, hrd, grav 841 846 IMPLICIT NONE 842 847 … … 859 864 ! --------------------------------------------------------------------- 860 865 861 include "cvthermo.h"862 866 include "cv30param.h" 863 867 include "conema3.h" … … 1212 1216 SUBROUTINE cv30_closure(nloc, ncum, nd, icb, inb, pbase, p, ph, tv, buoy, & 1213 1217 sig, w0, cape, m) 1218 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 1219 , clmci, eps, epsi, epsim1, ginv, hrd, grav 1214 1220 IMPLICIT NONE 1215 1221 … … 1220 1226 ! =================================================================== 1221 1227 1222 include "cvthermo.h"1223 1228 include "cv30param.h" 1224 1229 … … 1418 1423 u, v, tra, h, lv, qnk, hp, tv, tvp, ep, clw, m, sig, ment, qent, uent, & 1419 1424 vent, sij, elij, ments, qents, traent) 1425 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 1426 , clmci, eps, epsi, epsim1, ginv, hrd, grav 1420 1427 IMPLICIT NONE 1421 1428 … … 1426 1433 ! --------------------------------------------------------------------- 1427 1434 1428 include "cvthermo.h"1429 1435 include "cv30param.h" 1430 1436 … … 1816 1822 mp, rp, up, vp, trap, wt, water, evap, b & ! RomP-jyg 1817 1823 , wdtraina, wdtrainm) ! 26/08/10 RomP-jyg 1824 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 1825 , clmci, eps, epsi, epsim1, ginv, hrd, grav 1818 1826 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 1819 1827 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac 1820 1828 IMPLICIT NONE 1821 include "cvthermo.h"1822 1829 include "cv30param.h" 1823 1830 … … 2183 2190 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 2184 2191 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac 2192 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 2193 , clmci, eps, epsi, epsim1, ginv, hrd, grav 2185 2194 IMPLICIT NONE 2186 2195 2187 include "cvthermo.h"2188 2196 include "cv30param.h" 2189 2197 include "conema3.h" … … 3274 3282 ,cape,ep,hp,icb,inb,clw,nk,t,h,lv & 3275 3283 ,epmax_diag) 3284 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 3285 , clmci, eps, epsi, epsim1, ginv, hrd, grav 3276 3286 implicit none 3277 3287 … … 3281 3291 ! Toutes les autres variables fn de ep sont calcul�es plus bas. 3282 3292 3283 INCLUDE "cvthermo.h"3284 3293 INCLUDE "cv30param.h" 3285 3294 INCLUDE "conema3.h" -
LMDZ6/trunk/libf/phylmd/cv3_buoy.f90
r5268 r5276 10 10 ! ************************************************************** 11 11 12 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 13 , clmci, eps, epsi, epsim1, ginv, hrd, grav 12 14 IMPLICIT NONE 13 15 14 include "cvthermo.h"15 16 include "cv3param.h" 16 17 include "YOMCST2.h" -
LMDZ6/trunk/libf/phylmd/cv3_cine.f90
r5274 r5276 15 15 ! ************************************************************** 16 16 17 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 18 , clmci, eps, epsi, epsim1, ginv, hrd, grav 17 19 USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO & 18 20 , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA & … … 30 32 31 33 32 include "cvthermo.h"33 34 include "cv3param.h" 34 35 ! input: -
LMDZ6/trunk/libf/phylmd/cv3_enthalpmix.f90
r5274 r5276 11 11 ! ************************************************************** 12 12 13 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 14 , clmci, eps, epsi, epsim1, ginv, hrd, grav 13 15 USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO & 14 16 , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA & … … 34 36 ! =============================================================== 35 37 36 include "cvthermo.h"37 38 include "YOETHF.h" 38 39 -
LMDZ6/trunk/libf/phylmd/cv3_estatmix.f90
r5274 r5276 12 12 ! **************************************************************** 13 13 14 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 15 , clmci, eps, epsi, epsim1, ginv, hrd, grav 14 16 USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO & 15 17 , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA & … … 35 37 ! =============================================================== 36 38 37 include "cvthermo.h"38 39 include "YOETHF.h" 39 40 -
LMDZ6/trunk/libf/phylmd/cv3_routines.f90
r5275 r5276 184 184 SUBROUTINE cv3_incrcount(len, nd, delt, sig) 185 185 186 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 186 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 187 , clmci, eps, epsi, epsim1, ginv, hrd, grav 188 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 187 189 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac 188 190 IMPLICIT NONE … … 225 227 SUBROUTINE cv3_prelim(len, nd, ndp1, t, q, p, ph, & 226 228 lv, lf, cpn, tv, gz, h, hm, th) 229 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 230 , clmci, eps, epsi, epsim1, ginv, hrd, grav 227 231 IMPLICIT NONE 228 232 … … 248 252 REAL cpx(len, nd) 249 253 250 include "cvthermo.h"251 254 include "cv3param.h" 252 255 … … 314 317 USE add_phys_tend_mod, ONLY: fl_cor_ebil 315 318 USE print_control_mod, ONLY: prt_level 319 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 320 , clmci, eps, epsi, epsim1, ginv, hrd, grav 316 321 IMPLICIT NONE 317 322 … … 332 337 333 338 include "cv3param.h" 334 include "cvthermo.h"335 339 336 340 !inputs: … … 605 609 SUBROUTINE cv3_undilute1(len, nd, t, qs, gz, plcl, p, icb, tnk, qnk, gznk, & 606 610 tp, tvp, clw, icbs) 611 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 612 , clmci, eps, epsi, epsim1, ginv, hrd, grav 607 613 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 608 614 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac … … 622 628 ! ---------------------------------------------------------------- 623 629 624 include "cvthermo.h"625 630 include "cv3param.h" 626 631 … … 1141 1146 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 1142 1147 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac 1148 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 1149 , clmci, eps, epsi, epsim1, ginv, hrd, grav 1143 1150 IMPLICIT NONE 1144 1151 … … 1161 1168 ! --------------------------------------------------------------------- 1162 1169 1163 include "cvthermo.h"1164 1170 include "cv3param.h" 1165 1171 include "conema3.h" … … 2050 2056 pbase, p, ph, tv, buoy, & 2051 2057 sig, w0, cape, m, iflag) 2058 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 2059 , clmci, eps, epsi, epsim1, ginv, hrd, grav 2052 2060 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 2053 2061 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac … … 2060 2068 ! =================================================================== 2061 2069 2062 include "cvthermo.h"2063 2070 include "cv3param.h" 2064 2071 … … 2296 2303 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 2297 2304 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac 2305 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 2306 , clmci, eps, epsi, epsim1, ginv, hrd, grav 2298 2307 IMPLICIT NONE 2299 2308 … … 2303 2312 ! --------------------------------------------------------------------- 2304 2313 2305 include "cvthermo.h"2306 2314 include "cv3param.h" 2307 2315 … … 2717 2725 faci, b, sigd, & 2718 2726 wdtrainA, wdtrainS, wdtrainM) ! RomP 2727 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 2728 , clmci, eps, epsi, epsim1, ginv, hrd, grav 2719 2729 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 2720 2730 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac … … 2722 2732 IMPLICIT NONE 2723 2733 2724 include "cvthermo.h"2725 2734 include "cv3param.h" 2726 2735 include "nuage.h" … … 3466 3475 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 3467 3476 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac 3468 3477 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 3478 , clmci, eps, epsi, epsim1, ginv, hrd, grav 3469 3479 IMPLICIT NONE 3470 3480 3471 include "cvthermo.h"3472 3481 include "cv3param.h" 3473 3482 include "conema3.h" … … 4961 4970 ep, Vprecip, elij, clw, epmlmMm, eplaMm, & 4962 4971 icb, inb) 4972 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 4973 , clmci, eps, epsi, epsim1, ginv, hrd, grav 4963 4974 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 4964 4975 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac … … 5151 5162 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 5152 5163 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac 5164 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 5165 , clmci, eps, epsi, epsim1, ginv, hrd, grav 5153 5166 implicit none 5154 5167 … … 5158 5171 ! Toutes les autres variables fn de ep sont calcul�es plus bas. 5159 5172 5160 include "cvthermo.h"5161 5173 include "cv3param.h" 5162 5174 include "conema3.h" -
LMDZ6/trunk/libf/phylmd/cv3p1_closure.f90
r5274 r5276 19 19 ! ************************************************************** 20 20 21 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 22 , clmci, eps, epsi, epsim1, ginv, hrd, grav 21 23 USE print_control_mod, ONLY: prt_level, lunout 22 24 USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO & … … 34 36 IMPLICIT NONE 35 37 36 include "cvthermo.h"37 38 include "cv3param.h" 38 39 include "YOMCST2.h" -
LMDZ6/trunk/libf/phylmd/cv3p2_closure.f90
r5275 r5276 18 18 ! ************************************************************** 19 19 20 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 21 , clmci, eps, epsi, epsim1, ginv, hrd, grav 20 22 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 21 23 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac … … 35 37 IMPLICIT NONE 36 38 37 include "cvthermo.h"38 39 include "cv3param.h" 39 40 include "YOMCST2.h" -
LMDZ6/trunk/libf/phylmd/cv3p_mixing.f90
r5275 r5276 13 13 ! ************************************************************** 14 14 15 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 16 , clmci, eps, epsi, epsim1, ginv, hrd, grav 15 17 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 16 18 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac … … 21 23 IMPLICIT NONE 22 24 23 include "cvthermo.h"24 25 include "cv3param.h" 25 26 include "YOMCST2.h" -
LMDZ6/trunk/libf/phylmd/cv_driver.F90
r5275 r5276 682 682 SUBROUTINE cv_flag(iflag_ice_thermo) 683 683 684 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 685 , clmci, eps, epsi, epsim1, ginv, hrd, grav 684 686 USE cvflag_mod_h, ONLY: icvflag_Tpa, cvflag_grav, cvflag_ice, ok_optim_yield, ok_entrain, ok_homo_tend, & 685 687 ok_convstop, ok_intermittent, cvflag_prec_eject, qsat_depends_on_qt, adiab_ascent_mass_flux_depends_on_ejectliq, keepbug_ice_frac … … 699 701 ! 700 702 ! si icvflag_Tpa=0, alors la fraction de glace dans l'ascendance adiabatique est 701 ! fonction de la temperature de l'environnement et la temperature de l'ascendance est 702 ! calculee en deux it�rations, une en supposant qu'il n'y a pas de glace et l'autre 703 ! fonction de la temperature de l'environnement et la temperature de l'ascendance est 704 ! calculee en deux it�rations, une en supposant qu'il n'y a pas de glace et l'autre 703 705 ! en ajoutant la glace (ancien sch�ma d'Arnaud Jam). 704 706 ! si icvflag_Tpa=1, alors la fraction de glace dans l'ascendance adiabatique est 705 ! fonction de la temperature de l'environnement et la temperature de l'ascendance est 707 ! fonction de la temperature de l'environnement et la temperature de l'ascendance est 706 708 ! calculee en une seule iteration. 707 709 ! si icvflag_Tpa=2, alors la fraction de glace dans l'ascendance adiabatique est 708 ! fonction de la temperature de l'ascendance et la temperature de l'ascendance est 710 ! fonction de la temperature de l'ascendance et la temperature de l'ascendance est 709 711 ! calculee en une seule iteration. 710 icvflag_Tpa=0 712 icvflag_Tpa=0 711 713 call getin_p('icvflag_Tpa', icvflag_Tpa) 712 714 … … 728 730 , RESTT, RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS & 729 731 , RALPD, RBETD, RGAMD 732 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 733 , clmci, eps, epsi, epsim1, ginv, hrd, grav 730 734 IMPLICIT NONE 731 735 … … 735 739 736 740 737 include "cvthermo.h"738 741 739 742 INTEGER iflag_con -
LMDZ6/trunk/libf/phylmd/cv_routines.f90
r5268 r5276 3 3 4 4 SUBROUTINE cv_param(nd) 5 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 6 , clmci, eps, epsi, epsim1, ginv, hrd, grav 5 7 IMPLICIT NONE 6 8 … … 75 77 76 78 SUBROUTINE cv_prelim(len, nd, ndp1, t, q, p, ph, lv, cpn, tv, gz, h, hm) 79 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 80 , clmci, eps, epsi, epsim1, ginv, hrd, grav 77 81 IMPLICIT NONE 78 82 … … 93 97 REAL cpx(len, nd) 94 98 95 include "cvthermo.h"96 99 include "cvparam.h" 97 100 … … 254 257 SUBROUTINE cv_undilute1(len, nd, t, q, qs, gz, p, nk, icb, icbmax, tp, tvp, & 255 258 clw) 259 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 260 , clmci, eps, epsi, epsim1, ginv, hrd, grav 256 261 IMPLICIT NONE 257 262 258 include "cvthermo.h"259 263 include "cvparam.h" 260 264 … … 481 485 SUBROUTINE cv_undilute2(nloc, ncum, nd, icb, nk, tnk, qnk, gznk, t, q, qs, & 482 486 gz, p, dph, h, tv, lv, inb, inb1, tp, tvp, clw, hp, ep, sigp, frac) 487 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 488 , clmci, eps, epsi, epsim1, ginv, hrd, grav 483 489 IMPLICIT NONE 484 490 … … 492 498 ! FIND THE LEVEL OF NEUTRAL BUOYANCY 493 499 ! --------------------------------------------------------------------- 494 495 include "cvthermo.h"496 500 include "cvparam.h" 497 501 … … 764 768 SUBROUTINE cv_closure(nloc, ncum, nd, nk, icb, tv, tvp, p, ph, dph, plcl, & 765 769 cpn, iflag, cbmf) 770 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 771 , clmci, eps, epsi, epsim1, ginv, hrd, grav 766 772 IMPLICIT NONE 767 773 … … 781 787 REAL dtpbl(nloc), dtmin(nloc), tvpplcl(nloc), tvaplcl(nloc) 782 788 REAL work(nloc) 783 784 include "cvthermo.h"785 789 include "cvparam.h" 786 790 … … 847 851 h, lv, qnk, hp, tv, tvp, ep, clw, cbmf, m, ment, qent, uent, vent, nent, & 848 852 sij, elij) 853 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 854 , clmci, eps, epsi, epsim1, ginv, hrd, grav 849 855 IMPLICIT NONE 850 856 851 include "cvthermo.h"852 857 include "cvparam.h" 853 858 … … 1098 1103 SUBROUTINE cv_unsat(nloc, ncum, nd, inb, t, q, qs, gz, u, v, p, ph, h, lv, & 1099 1104 ep, sigp, clw, m, ment, elij, iflag, mp, qp, up, vp, wt, water, evap) 1105 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 1106 , clmci, eps, epsi, epsim1, ginv, hrd, grav 1100 1107 IMPLICIT NONE 1101 1102 1103 include "cvthermo.h"1104 1108 include "cvparam.h" 1105 1109 … … 1301 1305 ment, qent, uent, vent, nent, elij, tv, tvp, iflag, wd, qprime, tprime, & 1302 1306 precip, cbmf, ft, fq, fu, fv, ma, qcondc) 1307 USE cvthermo_mod_h, ONLY: cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 1308 , clmci, eps, epsi, epsim1, ginv, hrd, grav 1303 1309 IMPLICIT NONE 1304 1310 1305 include "cvthermo.h"1306 1311 include "cvparam.h" 1307 1312 -
LMDZ6/trunk/libf/phylmd/cvthermo_mod_h.f90
r5275 r5276 1 ! 2 ! $Header$ 3 ! 4 ! Thermodynamical constants for convectL: 1 ! Replaces cvthermo.h 5 2 6 real cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0 7 real clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl, clmci 8 real eps, epsi, epsim1 9 real ginv, hrd 10 real grav 3 MODULE cvthermo_mod_h 4 IMPLICIT NONE; PRIVATE 5 PUBLIC cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 6 , clmci, eps, epsi, epsim1, ginv, hrd, grav 11 7 12 COMMON /cvthermo/ cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl & 13 ,t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 14 ,clmci, eps, epsi, epsim1, ginv, hrd, grav 8 ! Thermodynamical constants for convectL: 9 REAL cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0 10 REAL clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl, clmci 11 REAL eps, epsi, epsim1 12 REAL ginv, hrd 13 REAL grav 15 14 16 !$OMP THREADPRIVATE(/cvthermo/) 15 !$OMP THREADPRIVATE(cpd, cpv, cl, ci, rrv, rrd, lv0, lf0, g, rowl, t0, clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl & 16 !$OMP , clmci, eps, epsi, epsim1, ginv, hrd, grav) 17 END MODULE cvthermo_mod_h
Note: See TracChangeset
for help on using the changeset viewer.