Changeset 6545

Timestamp:

07/11/23 15:41:04 (5 months ago)

Author:

acosce

Message:

update driver and def to be use with inca model

Location:

CONFIG/UNIFORM/v7/IPSLCM7/GENERAL

Files:

• DRIVER/dynamico.driver (modified) (1 diff)
• DRIVER/inca.driver (modified) (4 diffs)
• PARAM/inca.def (modified) (2 diffs)

CONFIG/UNIFORM/v7/IPSLCM7/GENERAL/DRIVER/dynamico.driver

@@ -24 +24 @@
 
 
-    if [ -f ${SUBMIT_DIR}/../../../bin/inca_${OptMode}_${ConfChem}.dat ];  then 
-        eval $(grep nbtrac ${SUBMIT_DIR}/../../../bin/inca_${OptMode}_${ConfChem}.dat) || nbtrac=1
+    if [ -f ${SUBMIT_DIR}/../../../bin/inca_${ConfChem}.dat ];  then 
+        eval $(grep nbtrac ${SUBMIT_DIR}/../../../bin/inca_${ConfChem}.dat) || nbtrac=1
     else
         nbtrac=1

CONFIG/UNIFORM/v7/IPSLCM7/GENERAL/DRIVER/inca.driver

@@ -20 +20 @@
     IGCM_debug_PushStack "CHM_Initialize"
 
-    if [ X$ResolAtm != X ] ; then
-        # ResolAtm is set in config.card
-        RESOL_ATM_3D=${ResolAtm}
-        RESOL_ATM_Z=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $3}' )
-        RESOL_ATM_X=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $1}' )
-        RESOL_ATM_Y=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $2}' )
-
-        RESOL_CHM=${RESOL_ATM_X}${RESOL_ATM_Y}-L${RESOL_ATM_Z}
-    elif [ -f ${SUBMIT_DIR}/../.resol ] ; then
-        [ -f ${SUBMIT_DIR}/../.resol ] && eval $(grep RESOL_ATM_3D ${SUBMIT_DIR}/../.resol) || RESOL_ATM_3D=96x95x19
-        RESOL_CHM_Z=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $3}' | awk "-F-" '{print $1}')
-    else
+    if [ X${config_ListOfComponents_ICO} = X ] ; then
+
+        if [ X$ResolAtm != X ] ; then
+            # ResolAtm is set in config.card
+            RESOL_ATM_3D=${ResolAtm}
+            RESOL_ATM_Z=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $3}' )
+            RESOL_ATM_X=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $1}' )
+            RESOL_ATM_Y=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $2}' )
+            
+            RESOL_CHM=${RESOL_ATM_X}${RESOL_ATM_Y}-L${RESOL_ATM_Z}
+        elif [ -f ${SUBMIT_DIR}/../.resol ] ; then
+            [ -f ${SUBMIT_DIR}/../.resol ] && eval $(grep RESOL_ATM_3D ${SUBMIT_DIR}/../.resol) || RESOL_ATM_3D=96x95x19
+            RESOL_CHM_Z=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $3}' | awk "-F-" '{print $1}')
+        else
         IGCM_debug_Exit "ResolAtm is not set in config.card and the .resol file does not exist."
         IGCM_debug_Verif_Exit
+        fi
+    else
+        RESOL_ATM_Z=${dynamico_UserChoices_RESOL_ATM_Z}
+        RESOL_NBP=${dynamico_UserChoices_RESOL_NBP}
+
     fi
 
@@ -332 +339 @@
 
 
+    # inca.def : regrid output files on regular grid or not 
+    if [ ! X${inca_UserChoices_xios_remap_output} = X ] ; then
+        IGCM_comp_modifyDefFile nonblocker inca.def xios_remap_output  ${inca_UserChoices_xios_remap_output}
+    else
+        IGCM_comp_modifyDefFile nonblocker inca.def xios_remap_output  y
+    fi
+
+    # inca.def : choose if we read input file with xios or netcdf 
+    if [ ! X${inca_UserChoices_xios_read} = X ] ; then
+        IGCM_comp_modifyDefFile nonblocker inca.def xios_read  ${inca_UserChoices_xios_read}
+    else
+        IGCM_comp_modifyDefFile nonblocker inca.def xios_read  y
+    fi
+
     # inca.def : Activate reading of 10m wind in files if LMDZ_10m_winds=n in inca.card
     if [ ! X${inca_UserChoices_LMDZ_10m_winds} = X ] ; then
@@ -419 +440 @@
     #  <context id="inca" src="./context_inca.xml"/>
     echo '<context id="inca" src="./context_inca.xml"/>' > add.tmp
+    if [ -f context_input_inca.xml ] ; then
+        echo '<context id="inca" src="./context_input_inca.xml"/>' >> add.tmp
+    fi
+    if [ -f context_sflx_inca.xml ] ; then
+        echo '<context id="inca" src="./context_sflx_inca.xml"/>' >> add.tmp
+    fi
+    if [ -f context_co2_inca.xml ] ; then
+        echo '<context id="inca" src="./context_co2_inca.xml"/>' >> add.tmp
+    fi
+    if [ -f context_oxydants_inca.xml ] ; then
+        echo '<context id="inca" src="./context_oxydants_inca.xml"/>' >> add.tmp
+    fi
+    if [ -f context_start_inca.xml ] ; then
+        if  [ -f chem_start.nc ] ; then 
+            echo '<context id="inca" src="./context_start_inca.xml"/>' >> add.tmp
+        fi 
+    fi
+    if [ -f context_dust_inca.xml ] ; then
+        echo '<context id="inca" src="./context_dust_inca.xml"/>' >> add.tmp
+    fi
+    if [ -f context_winds_inca.xml ] ; then
+        echo '<context id="inca" src="./context_winds_inca.xml"/>' >> add.tmp
+    fi
+
+    if [ -f context_sad_inca.xml ] ; then
+        echo '<context id="inca" src="./context_sad_inca.xml"/>' >> add.tmp
+    fi
+
+
     cp iodef.xml iodef.xml.tmp
     sed -e "/COMPONENT CONTEXT/r add.tmp" \
@@ -439 +489 @@
     fi
 
+    # remove unlimited from time axis for allowed xios to read all values in one time 
+    echo "retire UNLIMITED sur axe temps "
+    nccopy -u u10xios.nc tmp.nc
+    mv tmp.nc u10xios.nc
+    nccopy -u v10xios.nc tmp.nc
+    mv tmp.nc v10xios.nc
+    nccopy -u oxydants.nc tmp.nc
+    mv tmp.nc oxydants.nc
+    nccopy -u sflx_ANT_LR.nc tmp.nc
+    mv tmp.nc sflx_ANT_LR.nc
+    nccopy -u sflx_BBG_LR.nc tmp.nc
+    mv tmp.nc sflx_BBG_LR.nc
+    nccopy -u sflx_GHG_LR.nc tmp.nc
+    mv tmp.nc sflx_GHG_LR.nc
+    nccopy -u sflx_NAT_LR.nc tmp.nc
+    mv tmp.nc sflx_NAT_LR.nc
+    nccopy -u sflx_MISC_LR.nc tmp.nc
+    mv tmp.nc sflx_MISC_LR.nc
+    nccopy -u sad.nc tmp.nc
+    mv tmp.nc sad.nc
+
+
+    # Set interpolated output grid depending on resoultion, only if ICO component is set in config.card
+    if [ ! "X${config_ListOfComponents_ICO}" = "X" ] ; then   
+        if [ ${RESOL_NBP} == 40  ] ; then
+            valuei=144
+            valuej=143
+        elif [ ${RESOL_NBP} == 80  ] ; then
+            valuei=360
+            valuej=180
+        elif [ ${RESOL_NBP} == 160  ] ; then
+            valuei=720
+            valuej=360
+    elif [ ${RESOL_NBP} == 320  ] ; then
+        valuei=1440
+        valuej=720
+        else
+            IGCM_debug_Exit "RESOL_NBP=${RESOL_NBP} is not yet implemented in drivers"
+            IGCM_debug_Verif_Exit
+        fi
+        
+        IGCM_debug_Print 1 "Modify output grid parameters for current resolution in context_inca.xml"
+        cp context_inca.xml context_inca.xml.init
+        sed -e "s/ni_glo.*=.*\"144\"/ni_glo=\"${valuei}\"/" \
+            -e "s/nj_glo.*=.*\"143\"/nj_glo=\"${valuej}\"/" \
+            context_inca.xml.init > context_inca.xml
+        
+    fi
+
+
     IGCM_debug_PopStack "CHM_Update"
 }

CONFIG/UNIFORM/v7/IPSLCM7/GENERAL/PARAM/inca.def

@@ -2 +2 @@
 # parametres pour INCA
 #######################################################################
+#choose if we regrid or not output files on regular grid if we are running with ICO grid 
+xios_remap_output=_AUTO_: DEFAULT=y
+#choose if we read input file with xios or not (mandatory for unstructured grid)
+xios_read=_AUTO_:DEFAULT=y
 #LMDZ_10m_winds=y --> use 10m wind from LMDZ, =n --> use 10m wind from output file
 LMDZ_10m_winds=_AUTO_
-10m_winds_method=_AUTO_
 #choose if there is a feedback with the areosol effect
 feedb=_AUTO_
-#indicate if we are running with xios and/or ioipsl
-XIOS_INCA_OK=_AUTO_
-IOIPSL_INCA_OK=_AUTO_
-#weibull 10m wind parameter / def 4
-kref=_AUTO_: DEFAULT =4
 #threshold 10m wind parameter / def 0.78
 tref=_AUTO_: DEFAULT =0.78
 #seasalt_correctif / def = 1 
 ss_corr=_AUTO_: DEFAULT =1 
-#sedimentation 
-multilayer_sediment=_AUTO_
 #coupled model with orchidee
 CoupOrchInca=_AUTO_
@@ -34 +30 @@
 #choose which time interpolation you want to interpolate sflx file. ( 0 - 1 - 2) (no interpolation / point to point interpolation / Sheng & Zwiers interpolation) 
 emi_interp_time=_AUTO_
-
+#activate the flag cmip_xml if you are running a CMIP experiment
+cmip_xml=_AUTO_: DEFAULT =n