Changeset 6052 for CONFIG/UNIFORM

Timestamp:

01/28/22 15:46:28 (2 years ago)

Author:

acosce

Message:

Add possibility to nudge winds in dynamico

Location:

CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7

Files:

• EXPERIMENTS/ICOLMDZORINCA/GES/COMP/dynamico.card (modified) (3 diffs)
• GENERAL/DRIVER/dynamico.driver (modified) (4 diffs)
• GENERAL/PARAM/context_dynamico.xml (added)
• GENERAL/PARAM/nudging_dynamico.xml (added)
• GENERAL/PARAM/run_dynamico.def_nbp40 (modified) (1 diff)

CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/EXPERIMENTS/ICOLMDZORINCA/GES/COMP/dynamico.card

@@ -13 +13 @@
 iflag_phys=1
 
+ok_guide=n
 
 [InitialStateFiles]
@@ -19 +20 @@
 [BoundaryFiles]
 List=()
+
+#List=($CCCSTOREDIR/../../subipsl/subipsl/ERA5/NETCDF/GLOBAL_025/4xdaily/AN_PL/${year}/u.${year}${month}.aphe5.GLOBAL_025.nc, u.nc),\
+#     ($CCCSTOREDIR/../../subipsl/subipsl/ERA5/NETCDF/GLOBAL_025/4xdaily/AN_PL/${year}/v.${year}${month}.aphe5.GLOBAL_025.nc, v.nc)
+
 ListNonDel= ()
 
@@ -29 +34 @@
         (${SUBMIT_DIR}/PARAM/run_dynamico.def_nbp${RESOL_NBP}, run_dynamico.def),\
         (${SUBMIT_DIR}/PARAM/iodef.xml, .),\
-        (${MODIPSL}/modeles/DYNAMICO/XML/context_dynamico.xml, .),\
-        (${MODIPSL}/modeles/DYNAMICO/XML/context_input_dynamico.xml, .),\
-        (${MODIPSL}/modeles/DYNAMICO/XML/field_def_dynamico.xml, .),\
+        (${SUBMIT_DIR}/PARAM/context_dynamico.xml, .),\
+        (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\
+        (${MODIPSL}/modeles/DYNAMICO/xml/field_def_dynamico.xml, .),\
         (${MODIPSL}/modeles/DYNAMICO/xml/file_def_dynamico.xml, .)
+#       (${SUBMIT_DIR}/PARAM/nudging_dynamico.xml, .),\
+
+
+
+#       (${MODIPSL}/modeles/DYNAMICO/xml/context_dynamico.xml, .),\
+#       (${MODIPSL}/modeles/DYNAMICO/xml/context_input_dynamico.xml, .),\
+#       (${MODIPSL}/modeles/DYNAMICO/xml/field_def_dynamico.xml, .),\
+#       (${MODIPSL}/modeles/DYNAMICO/xml/file_def_dynamico.xml, .)
 
 

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

@@ -32 +32 @@
 
 
+    if [ ! -f ${SUBMIT_DIR}/../../../DYNAMICO/xml/file_def_dynamico.xml ];  then 
+        module load python3
+        python3 ${SUBMIT_DIR}/../../../DYNAMICO/xml/file_def_dynamico.xml.py > ${SUBMIT_DIR}/../../../DYNAMICO/xml/file_def_dynamico.xml
+    fi
+
     IGCM_debug_PopStack "ICO_Initialize"
 }
@@ -54 +59 @@
     # Set lenght for execution 
     IGCM_comp_modifyDefFile blocker run_dynamico.def run_length $(( ${PeriodLengthInDays} * 86400 ))
+
+    #  Activate nudging if ok_guide=y set in dynamico.card
+    if [  X${dynamico_UserChoices_ok_guide} = Xy ] ; then
+        IGCM_comp_modifyDefFile nonblocker run_dynamico.def guided_type  nudging
+    else
+        IGCM_comp_modifyDefFile nonblocker run_dynamico.def guided_type  none
+    fi
 
     # Change run_dynamico.def according to UserChoices options set in dynmaico.card
@@ -109 +121 @@
     fi
 
+    if [ -f nudging_dynamico.xml ] ; then
+        echo '<context id="icosagcm" src="./nudging_dynamico.xml" />'  >> add.tmp
+    fi
+
     if [ -f file_def_dynamico.xml ] ; then
         echo '<context id="icosagcm" src="./file_def_dynamico.xml" />' >> add.tmp
@@ -137 +153 @@
       rm add.tmp
     fi
+
+    #remove unlimited time axis 
+    if [ -f u.nc ] ; then 
+        nccopy -u u.nc tmp.nc
+        mv tmp.nc u.nc
+    fi 
+    if [ -f v.nc ] ; then 
+        nccopy -u v.nc tmp.nc
+        mv tmp.nc v.nc
+    fi 
+
 
     if [ X${config_ListOfComponents_CPL} = Xoasis ] ; then

CONFIG/UNIFORM/v7/ICOLMDZORINCA_v7/GENERAL/PARAM/run_dynamico.def_nbp40

@@ +1 @@
+metric_type = icosa_global
+
+#---------------- guided ----------------
+#guided_type = none
+guided_type = _AUTO_
+#nudging
+#frequence de lecture dans fichier forcage / coller a frequence du fichier par rapport à dt
+# toutes les 6h 
+itau_nudging = 48
+
+nudging_zone = global
+
+
+#temps relaxation en secondes = force de rappel (voir dissip/nudging.f90)  >= dt
+U_relax_in=21600
+PS_relax_in=1e30
+T_relax_in=1e30
+Q_relax_in=1e30
+
+#guide_U = u, v 
+guide_U=y
+
+# -------------------------------  ---------------------------------
+
 mpi_threading_mode=funneled
 # ------------------------------- Mesh ---------------------------------