Changeset 4743 for LMDZ6/branches/Portage_acc/libf/phylmd/conf_phys_m.F90

Timestamp:

Oct 30, 2023, 5:37:00 PM (7 months ago)

Author:

Laurent Fairhead

Message:

Merge of ACC branch with 4740 revision from trunk

Location:

LMDZ6/branches/Portage_acc

Files:

• . (modified) (1 prop)
• libf/phylmd/conf_phys_m.F90 (modified) (16 diffs)

LMDZ6/branches/Portage_acc/libf/phylmd/conf_phys_m.F90

@@ -16 +16 @@
        fact_cldcon, facttemps,ok_newmicro,iflag_radia,&
        iflag_cld_th, &
-       iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
+       ratqsbas,ratqshaut,tau_ratqs, &
        ok_ade, ok_aie, ok_alw, ok_cdnc, ok_volcan, flag_volc_surfstrat, aerosol_couple, & 
        chemistry_couple, flag_aerosol, flag_aerosol_strat, flag_aer_feedback, &
@@ -29 +29 @@
     USE carbon_cycle_mod,  ONLY: read_fco2_ocean_cor,var_fco2_ocean_cor
     USE carbon_cycle_mod,  ONLY: read_fco2_land_cor,var_fco2_land_cor
+    USE chemistry_cycle_mod, ONLY: dms_cycle_cpl
     USE mod_grid_phy_lmdz, ONLY: klon_glo
     USE print_control_mod, ONLY: lunout
     use config_ocean_skin_m, only: config_ocean_skin
     USE phys_state_var_mod, ONLY: phys_tstep
+    USE infotrac_phy, ONLY: type_trac
 
     INCLUDE "conema3.h"
@@ -85 +87 @@
     REAL                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut,tau_ratqs
     INTEGER              :: iflag_cld_th
-    INTEGER              :: iflag_ratqs
+     
 
     CHARACTER (len = 6), SAVE  :: type_ocean_omp, version_ocean_omp, ocean_omp
@@ -173 +175 @@
     LOGICAL,SAVE :: ok_ice_sursat_omp
     LOGICAL,SAVE :: ok_plane_h2o_omp, ok_plane_contrail_omp
-    INTEGER,SAVE :: iflag_rei_omp
     REAL,SAVE :: rad_froid_omp, rad_chau1_omp, rad_chau2_omp
-    REAL,SAVE :: rei_min_omp, rei_max_omp
     INTEGER,SAVE :: iflag_sic_omp, iflag_inertie_omp
     REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_sic_omp
@@ -233 +233 @@
     LOGICAL, SAVE :: carbon_cycle_cpl_omp
     LOGICAL, SAVE :: carbon_cycle_rad_omp
+    LOGICAL, SAVE :: dms_cycle_cpl_omp
     INTEGER, SAVE :: level_coupling_esm_omp
     LOGICAL, SAVE :: read_fco2_ocean_cor_omp
@@ -242 +243 @@
     LOGICAL, SAVE :: ok_new_lscp_omp
     LOGICAL, SAVE :: ok_icefra_lscp_omp
+    !rajout de choix_bulk et nit_bulk par Olivier Torres
+    INTEGER,SAVE  :: choix_bulk_omp
+    INTEGER,SAVE  :: nit_bulk_omp
+    INTEGER,SAVE  :: kz0_omp
     LOGICAL, SAVE :: ok_bs_omp, ok_rad_bs_omp
 
@@ -935 +940 @@
     nbapp_rad_omp = 12
     CALL getin('nbapp_rad',nbapp_rad_omp)
+
+    !rajout Olivier Torres
+    !Config  Key  = choix_bulk
+    !Config  Desc = choix de la formulation bulk a prendre dans clcdrag au-dessus de l'ocean
+    !Config  Def  = 0
+    !Config         0 -> originale (lmdz/Louis 79)
+    !Config         1 -> COARE
+    !Config         2 -> CORE-"pure" (cf. Large)
+    !Config         3 -> CORE-"mixte" (avec z_0 et C_T^N donnees par Smith 88)
+    choix_bulk_omp = 0
+    CALL getin('choix_bulk',choix_bulk_omp)
+
+    !Config  Key  = nit_bulk
+    !Config  Desc = choix du nombre d'it de pt fixe dans la bulk
+    !Config  Def  = 5
+    nit_bulk_omp = 1
+    CALL getin('nit_bulk',nit_bulk_omp)
+
+    !Config  Key  = kz0
+    !Config  Desc = choix de la formulation z0 pour la bulk ECUME
+    !Config  Def  = 1
+    !Config         0 -> ARPEGE formulation
+    !Config         1 -> Smith Formulation
+    !Config         2 -> Direct computation using the stability functions
+    kz0_omp = 0
+    CALL getin('kz0',kz0_omp)
+
 
     !Config  Key  = iflag_con
@@ -1305 +1337 @@
 
 
-    !Config Key  = iflag_rei
-    !Config Desc =  
-    !Config Def  = 0
-    !Config Help = 
-    !
-    iflag_rei_omp = 0
-    CALL getin('iflag_rei',iflag_rei_omp)
-
-
     !
     !Config Key  = iflag_ice_thermo
@@ -1347 +1370 @@
     ok_plane_contrail_omp = .FALSE.
     CALL getin('ok_plane_contrail',ok_plane_contrail_omp)
-
-    !
-    !Config Key  = rei_min
-    !Config Desc =  
-    !Config Def  = 3.5
-    !Config Help = 
-    !
-    rei_min_omp = 3.5
-    CALL getin('rei_min',rei_min_omp)
-
-    !
-    !Config Key  = rei_max
-    !Config Desc =  
-    !Config Def  = 61.29
-    !Config Help = 
-    !
-    rei_max_omp = 61.29
-    CALL getin('rei_max',rei_max_omp)
 
     !
@@ -2240 +2245 @@
     CALL getin('level_coupling_esm',level_coupling_esm_omp)
 
+
+    dms_cycle_cpl_omp=.FALSE.
+    CALL getin('transm_dms_oa',dms_cycle_cpl_omp)
+    !
+    !Config Key  = dms_cycle_cpl
+    !Config Desc = receive dms from pisces via oasis
+    !Config Def  = .FALSE.
+    !Config Help = can be activated only if lmdz use inca for the atmospheric chemistry 
+
+    
     !$OMP END MASTER
     !$OMP BARRIER
@@ -2305 +2320 @@
     rad_chau1 = rad_chau1_omp
     rad_chau2 = rad_chau2_omp
-    iflag_rei=iflag_rei_omp
     iflag_ice_thermo = iflag_ice_thermo_omp
     ok_ice_sursat = ok_ice_sursat_omp
     ok_plane_h2o = ok_plane_h2o_omp
     ok_plane_contrail = ok_plane_contrail_omp
-    rei_min = rei_min_omp
-    rei_max = rei_max_omp
     top_height = top_height_omp
     overlap = overlap_omp
@@ -2514 +2526 @@
     read_fco2_land_cor = read_fco2_land_cor_omp
     var_fco2_land_cor = var_fco2_land_cor_omp
+    dms_cycle_cpl = dms_cycle_cpl_omp
+    !rajout Olivier Torres
+    kz0=kz0_omp
+    choix_bulk = choix_bulk_omp
+    nit_bulk = nit_bulk_omp
 
     ! Test of coherence between type_ocean and version_ocean
@@ -2633 +2650 @@
     ENDIF
 
+
+    ! Test on chemistry cycle 
+    IF ((type_trac .ne. "inca" .AND. type_trac .ne. "inco") .AND. dms_cycle_cpl ) THEN 
+       CALL abort_physic('conf_phys', 'dms_cycle_cpl has to be TRUE only with INCA coupling model',1)
+    ENDIF
+    
     ! ORCHIDEE must be activated for ifl_pbltree=1
     IF (.NOT. ok_veget .AND. ifl_pbltree==1) THEN
@@ -2719 +2742 @@
     WRITE(lunout,*) ' rad_chau1 = ',rad_chau1
     WRITE(lunout,*) ' rad_chau2 = ',rad_chau2
-    WRITE(lunout,*) ' iflag_rei = ',iflag_rei
     WRITE(lunout,*) ' iflag_ice_thermo = ',iflag_ice_thermo
     WRITE(lunout,*) ' ok_ice_sursat = ',ok_ice_sursat
     WRITE(lunout,*) ' ok_plane_h2o = ',ok_plane_h2o
     WRITE(lunout,*) ' ok_plane_contrail = ',ok_plane_contrail
-    WRITE(lunout,*) ' rei_min = ',rei_min
-    WRITE(lunout,*) ' rei_max = ',rei_max
     WRITE(lunout,*) ' overlap = ',overlap 
     WRITE(lunout,*) ' cdmmax = ',cdmmax 
@@ -2844 +2864 @@
     WRITE(lunout,*) ' read_fco2_land_cor = ', read_fco2_land_cor
     WRITE(lunout,*) ' var_fco2_land_cor = ', var_fco2_land_cor
+    WRITE(lunout,*) ' dms_cycle_cpl = ', dms_cycle_cpl
     WRITE(lunout,*) ' iflag_tsurf_inlandsis = ', iflag_tsurf_inlandsis
     WRITE(lunout,*) ' iflag_temp_inlandsis = ', iflag_temp_inlandsis
@@ -2861 +2882 @@
     WRITE(lunout,*) ' buf_sph_pol = ', buf_sph_pol
     WRITE(lunout,*) ' buf_siz_pol= ', buf_siz_pol
+    !rajout Olivier Torres
+    write(lunout,*) 'choix_bulk = ', choix_bulk
+    write(lunout,*) 'nit_bulk = ', nit_bulk
+    write(lunout,*) 'kz0 = ', kz0
 
     !$OMP END MASTER