Changeset 6058 for LMDZ6/trunk/libf/phylmd/cv_driver_mod.F90

Timestamp:

Jan 29, 2026, 12:52:39 AM (6 weeks ago)

Author:

fhourdin

Message:

Travail pour la replayisation de la convection

Reunion de tous les anciens common devenus modules, dans lmdz_cv_ini.
Déplacement de presque toutes les routines d'initialisation dans lmdz_cv_ini.
Encapsulage de certains sous-programmes dans des modules.
Suppression de programmes inutilisés (cv3_crit et cv3_incp)
Reste :

• à sortir des routines d'initialisation "_pre" de cv_driver et

cva_driver

• à passer le variables argunement en intent(in/out/inout).

La convergence numérique a été testée pour
iflag_con=3/30/4
en 3D parallèle.
La compilation de la version isotopique fonctionne.

File:

• LMDZ6/trunk/libf/phylmd/cv_driver_mod.F90 (moved) (moved from LMDZ6/trunk/libf/phylmd/cv_driver.F90) (4 diffs)

LMDZ6/trunk/libf/phylmd/cv_driver_mod.F90

@@ -1 @@
-
 ! $Header$
+MODULE cv_driver_mod
+
+CONTAINS
 
 SUBROUTINE cv_driver(len, nd, ndp1, ntra, iflag_con, t1, q1, qs1, u1, v1, &
@@ -12 +14 @@
 
   USE dimphy
-  USE cv30_routines_mod, ONLY: cv30_param, cv30_prelim, cv30_feed, cv30_undilute1, cv30_trigger, cv30_compress, cv30_undilute2, &
+  USE cv30_routines_mod, ONLY: cv30_prelim, cv30_feed, cv30_undilute1, cv30_trigger, cv30_compress, cv30_undilute2, &
           cv30_closure, cv30_epmax_fn_cape, cv30_mixing, cv30_unsat, cv30_yield, cv30_tracer, cv30_uncompress
-  USE cv_routines_mod
+  USE cv_routines_mod, ONLY : cv_prelim, cv_feed, cv_undilute1, cv_trigger, cv_compress, cv_undilute2, cv_closure, cv_mixing, cv_unsat, cv_yield, cv_uncompress
+
+  USE lmdz_cv_ini, ONLY : cv_param, cv_flag, cv_thermo ! A retirer apres replayisation
+
   IMPLICIT NONE
 
@@ -352 +357 @@
   ! -------------------------------------------------------------------
   ! --- SET CONSTANTS AND PARAMETERS
-  ! -------------------------------------------------------------------
-  ! print *, '-> cv_driver'      !jyg
-  ! -- set simulation flags:
-  ! (common cvflag)
-
-  CALL cv_flag(0)
-
-  ! -- set thermodynamical constants:
-  ! (common cvthermo)
-
-  CALL cv_thermo(iflag_con)
-
-  ! -- set convect parameters
-
-  ! includes microphysical parameters and parameters that
-  ! control the rate of approach to quasi-equilibrium)
-  ! (common cvparam)
-
-
-  IF (iflag_con==30) THEN
-    CALL cv30_param(nd, delt)
-  END IF
-
-  IF (iflag_con==4) THEN
-    CALL cv_param(nd)
-  END IF
+! Replay deplace vers lmdz_cv_ini  ! -------------------------------------------------------------------
+! Replay deplace vers lmdz_cv_ini  ! print *, '-> cv_driver'      !jyg
+! Replay deplace vers lmdz_cv_ini  ! -- set simulation flags:
+! Replay deplace vers lmdz_cv_ini  ! (common cvflag)
+! Replay deplace vers lmdz_cv_ini
+! Replay deplace vers lmdz_cv_ini  CALL cv_flag(0)
+! Replay deplace vers lmdz_cv_ini
+! Replay deplace vers lmdz_cv_ini  ! -- set thermodynamical constants:
+! Replay deplace vers lmdz_cv_ini  ! (common cvthermo)
+! Replay deplace vers lmdz_cv_ini
+! Replay deplace vers lmdz_cv_ini  CALL cv_thermo(iflag_con)
+! Replay deplace vers lmdz_cv_ini
+! Replay deplace vers lmdz_cv_ini  ! -- set convect parameters
+! Replay deplace vers lmdz_cv_ini
+! Replay deplace vers lmdz_cv_ini  ! includes microphysical parameters and parameters that
+! Replay deplace vers lmdz_cv_ini  ! control the rate of approach to quasi-equilibrium)
+! Replay deplace vers lmdz_cv_ini  ! (common cvparam)
+! Replay deplace vers lmdz_cv_ini
+! Replay deplace vers lmdz_cv_ini
+! Replay deplace vers lmdz_cv_ini  IF (iflag_con==30) THEN
+! Replay deplace vers lmdz_cv_ini    CALL cv30_param(nd, delt)
+! Replay deplace vers lmdz_cv_ini  END IF
+! Replay deplace vers lmdz_cv_ini
+! Replay deplace vers lmdz_cv_ini  IF (iflag_con==4) THEN
+! Replay deplace vers lmdz_cv_ini    CALL cv_param(nd)
+! Replay deplace vers lmdz_cv_ini  END IF
 
   ! ---------------------------------------------------------------------
@@ -683 +688 @@
 END SUBROUTINE cv_driver
 
-! ==================================================================
-SUBROUTINE cv_flag(iflag_ice_thermo)
-
-  USE cvflag_mod_h
-  USE ioipsl_getin_p_mod, ONLY : getin_p
-
-  IMPLICIT NONE
-
-  ! Argument : iflag_ice_thermo : ice thermodynamics is taken into account if
-  ! iflag_ice_thermo >=1
-  INTEGER iflag_ice_thermo
-
-
-  ! -- si .TRUE., on rend la gravite plus explicite et eventuellement
-  ! differente de 10.0 dans convect3:
-  cvflag_grav = .TRUE.
-  cvflag_ice = iflag_ice_thermo >= 1
-  !
-! si icvflag_Tpa=0, alors la fraction de glace dans l'ascendance adiabatique est
-  ! fonction de la temperature de l'environnement et la temperature de l'ascendance est
-  ! calculee en deux it�rations, une en supposant qu'il n'y a pas de glace et l'autre
-  ! en ajoutant la glace (ancien sch�ma d'Arnaud Jam).
-! si icvflag_Tpa=1, alors la fraction de glace dans l'ascendance adiabatique est
-  ! fonction de la temperature de l'environnement et la temperature de l'ascendance est
-  ! calculee en une seule iteration.
-! si icvflag_Tpa=2, alors la fraction de glace dans l'ascendance adiabatique est
-  ! fonction de la temperature de l'ascendance et la temperature de l'ascendance est
-  ! calculee en une seule iteration.
-  icvflag_Tpa=0
-  call getin_p('icvflag_Tpa', icvflag_Tpa)
-
-  RETURN
-END SUBROUTINE cv_flag
-
-! ==================================================================
-SUBROUTINE cv_thermo(iflag_con)
-  USE yomcst_mod_h
-  USE lmdz_cv_ini, ONLY : ci,cl,clmci,clmcpd,clmcpv,cpd,cpdmcp,cpv,cpvmcl,cpvmcpd,eps,epsi,epsim1,g,ginv,grav,hrd,lf0,lv0,rowl,rrd,rrv,t0
-
-IMPLICIT NONE
-
-  ! -------------------------------------------------------------
-  ! Set thermodynamical constants for convectL
-  ! -------------------------------------------------------------
-
-
-
-  INTEGER iflag_con
-
-
-  ! original set from convect:
-  IF (iflag_con==4) THEN
-    cpd = 1005.7
-    cpv = 1870.0
-    cl = 4190.0
-    rrv = 461.5
-    rrd = 287.04
-    lv0 = 2.501E6
-    g = 9.8
-    t0 = 273.15
-    grav = g
-  ELSE
-
-    ! constants consistent with LMDZ:
-    cpd = rcpd
-    cpv = rcpv
-    cl = rcw
-    ci = rcs
-    rrv = rv
-    rrd = rd
-    lv0 = rlvtt
-    lf0 = rlstt - rlvtt
-    g = rg ! not used in convect3
-    ! ori      t0  = RTT
-    t0 = 273.15 ! convect3 (RTT=273.16)
-    ! maf       grav= 10.    ! implicitely or explicitely used in convect3
-    grav = g ! implicitely or explicitely used in convect3
-  END IF
-
-  rowl = 1000.0 !(a quelle variable de YOMCST cela correspond-il?)
-
-  clmcpv = cl - cpv
-  clmcpd = cl - cpd
-  clmci = cl - ci
-  cpdmcp = cpd - cpv
-  cpvmcpd = cpv - cpd
-  cpvmcl = cl - cpv ! for convect3
-  eps = rrd/rrv
-  epsi = 1.0/eps
-  epsim1 = epsi - 1.0
-  ! ginv=1.0/g
-  ginv = 1.0/grav
-  hrd = 0.5*rrd
-
-  RETURN
-END SUBROUTINE cv_thermo
+END MODULE cv_driver_mod