Changeset 4690 for LMDZ6/trunk/libf/phylmdiso

Timestamp:

Sep 17, 2023, 1:16:26 PM (16 months ago)

Author:

fhourdin

Message:

Preparation/test de convergence num en temps

File:

• LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -83 +83 @@
     USE yamada_ini_mod, ONLY : yamada_ini
     USE lmdz_atke_turbulence_ini, ONLY : atke_ini
-    USE lmdz_thermcell_ini, ONLY : thermcell_ini
+    USE lmdz_thermcell_ini, ONLY : thermcell_ini, iflag_thermals_tenv
     USE lmdz_thermcell_dtke, ONLY : thermcell_dtke
     USE blowing_snow_ini_mod, ONLY : blowing_snow_ini , qbst_bs 
@@ -4465 +4465 @@
 #endif
           !>jyg
+          ! WARNING !!! introduction d'une subtilite avec la distinction
+          ! entre deux types de profiles pour les thermiques.
+          ! Ne change pas les resultats avec les isotopes tant qu'on
+          ! tourne avec le defaut : iflag_thermals_tenv=0
+          if ( iflag_thermals_tenv == 0 ) then 
+              abort_message ='iflag_thermals_env/=0 : non prevu avec les isotopes'
+              CALL abort_physic (modname,abort_message,1)
+          endif
+
           CALL calltherm(pdtphys &
                ,pplay,paprs,pphi,weak_inversion &
                         ! ,u_seri,v_seri,t_seri,q_seri,zqsat,debut & !jyg
-               ,u_therm,v_therm,t_therm,q_therm,zqsat,debut &  !jyg
+               ,u_therm,v_therm,t_therm,q_therm,t_therm,q_therm,zqsat,debut &  !jyg
                ,d_u_ajs,d_v_ajs,d_t_ajs,d_q_ajs &
                ,fm_therm,entr_therm,detr_therm &