Changeset 6035 for CONFIG/UNIFORM/v6/IPSLCM5A2CHT.2/GENERAL/DRIVER

Timestamp:

01/07/22 11:35:28 (2 years ago)

Author:

acosce

Message:

in IPSLCM5A2CHT configuration :

• add parameters files need by paleo Experiment
• Add Addnoise functionality in configuration (see rev[5106] on IPSLCM5A2 configuration)
• remove .resol file not need
• add modification in driver for paleo experiment
• add iom.F90 to allow flexible handling of calendar (see rev[4930] in IPSLCM5A2CHT configuration)

Location:

CONFIG/UNIFORM/v6/IPSLCM5A2CHT.2/GENERAL/DRIVER

Files:

• lmdz.driver (modified) (1 diff)
• oasis.driver (modified) (2 diffs)
• opa9.driver (modified) (3 diffs)
• orchidee.driver (modified) (1 diff)

CONFIG/UNIFORM/v6/IPSLCM5A2CHT.2/GENERAL/DRIVER/lmdz.driver

@@ -279 +279 @@
     ## gcm.def parameters : 
     # Modification of variable iphysiq depending on the physics
-    if [ X${LMDZ_Physics} = X"AP" ] ; then
+    ISPALEO=${config_UserChoices_ExpType##*/}
+
+    if [ X${LMDZ_Physics} = X"AP" ] && [ ${ISPALEO} != "paleo" ]; then
         if [ ${RESOL_ATM_3D} = 96x95x39 ]; then
             value=10

CONFIG/UNIFORM/v6/IPSLCM5A2CHT.2/GENERAL/DRIVER/oasis.driver

@@ -25 +25 @@
         esac
     done
+
+    # Initialisation required for AddNoise done during one Period
+    Test_AddNoise_First_Only=y
 
     IGCM_debug_PopStack "CPL_Initialize"
@@ -95 +98 @@
     IGCM_comp_modifyXmlFile force iodef.xml using_oasis NONE true
 
+    if [ X"${Test_AddNoise_First_Only}" = X"y" ]; then
+
+    Test_AddNoise_First_Only=n
+
+    # AddNoise if required but only for one Period
+
+    # Read ByPass_addnoise_sst option in oasis.card --> apply addnoise for sst in CPL restart file if [ $ByPass_addnoise_sst = y ]
+    ByPass_addnoise_sst=${oasis_UserChoices_ByPass_addnoise_sst}
+    if [ X"${ByPass_addnoise_sst}" = X"y" ] ; then
+
+      ByPass_PerturbExe=${oasis_UserChoices_ByPass_PerturbExe}
+      ByPass_FileOut=${oasis_UserChoices_ByPass_FileOut}
+      ByPass_PerturbVar=${oasis_UserChoices_ByPass_PerturbVar}
+      ByPass_PerturbAmp=${oasis_UserChoices_ByPass_PerturbAmp}
+
+      IGCM_sys_Chmod 644 ${RUN_DIR}/${ByPass_FileOut}.nc
+      IGCM_debug_Print 1 "${ByPass_PerturbExe} ${RUN_DIR}/${ByPass_FileOut}.nc ${ByPass_PerturbVar} ${ByPass_PerturbAmp}"
+
+      ${ByPass_PerturbExe} ${RUN_DIR}/${ByPass_FileOut}.nc ${ByPass_PerturbVar} ${ByPass_PerturbAmp}
+      if [ $? -ne 0 ] ; then
+        IGCM_debug_Exit "Error with $( basename ${ByPass_PerturbExe} )"
+        IGCM_debug_Verif_Exit
+      fi
+
+      echo
+      IGCM_debug_Print 1 "ByPass_addnoise_sst : applied addnoise for sst in CPL restart file"
+      echo
+      ByPass_addnoise_sst=n
+      # necessary to reset ByPass in next Period during the same Job (PeriodNb<>1)
+      export oasis_UserChoicesByPass_addnoise_sst=n
+      IGCM_card_WriteOption ${SUBMIT_DIR}/COMP/oasis.card UserChoices ByPass_addnoise_sst "${ByPass_addnoise_sst}"
+    fi
+
+    fi
+
     IGCM_debug_PopStack "CPL_Update"
 }

CONFIG/UNIFORM/v6/IPSLCM5A2CHT.2/GENERAL/DRIVER/opa9.driver

@@ -26 +26 @@
     supergrep () { grep "^ *$1 *=" $2 | sed -e "s% *\!.*%%" -e "s%^ *$1 *=%%" ; }
     ##--Variables used by OPA --
-    NAMELIST_OPA_CFG=${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE}_cfg
+    ISPALEO=${config_UserChoices_ExpType##*/}
+    if [ "${ISPALEO}" = "paleo" ] ; then
+        NAMELIST_OPA_CFG=${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE}_cfg_paleo
+    else
+        NAMELIST_OPA_CFG=${SUBMIT_DIR}/PARAM/namelist_${RESOL_OCE}_cfg
+    fi
     IGCM_debug_Print 1 "namelist_cfg : ${NAMELIST_OPA_CFG}"
     if [ ! -r ${NAMELIST_OPA_CFG} ] ; then
@@ -38 +43 @@
     # OPA_NN_FSBC=$(   supergrep nn_fsbc       namelist_cfg )
     (( OPA_NPDT_DAY   =  86400 / OPA_RDT        ))
-    (( OPA_NPDT_YEAR  =  OPA_NPDT_DAY * 365     ))
     (( OPA_RDT_SBC    =  OPA_RDT * OPA_NN_FSBC ))
+    # nleapy configuration + OPA_NPDT_YEAR
+    case ${config_UserChoices_CalendarType} in
+        ( leap|gregorian) OPA_NLEAPY=1
+                          (( OPA_NPDT_YEAR  =  OPA_NPDT_DAY * 365     )) ;;
+        ( noleap        ) OPA_NLEAPY=0
+                          (( OPA_NPDT_YEAR  =  OPA_NPDT_DAY * 365     )) ;;
+        ( 360d          ) OPA_NLEAPY=30
+                          (( OPA_NPDT_YEAR  =  OPA_NPDT_DAY * 360     )) ;;
+        ( *             ) OPA_NLEAPY=30
+                          (( OPA_NPDT_YEAR  =  OPA_NPDT_DAY * 360     )) ;;
+    esac
  
     # Period Length In Days between DateBegin and first day of calendar 0001 01 01
@@ -145 +160 @@
     fi
 
-    # nleapy configuration
-    case ${config_UserChoices_CalendarType} in
-        ( leap|gregorian) OPA_NLEAPY=1  ;;
-        ( noleap        ) OPA_NLEAPY=0  ;;
-        ( 360d          ) OPA_NLEAPY=30 ;;
-        ( *             ) OPA_NLEAPY=30 ;;
-    esac
 
     if ( [ "${opa9_UserChoices_Reproducibility_after_restart}" = "y" ] ) ; then

CONFIG/UNIFORM/v6/IPSLCM5A2CHT.2/GENERAL/DRIVER/orchidee.driver

@@ -55 +55 @@
     # Activate creation of river description file only for the first period
     if [ ${CumulPeriod} -eq 1 ] ; then
-        IGCM_comp_modifyDefFile nonblocker orchidee.def RIVER_DESC y
+        ISPALEO=${config_UserChoices_ExpType##*/}
+        if [ "${ISPALEO}" = "paleo" ] ; then # change to y if needed for paleo
+            IGCM_comp_modifyDefFile nonblocker orchidee.def RIVER_DESC n
+        else
+            IGCM_comp_modifyDefFile nonblocker orchidee.def RIVER_DESC y
+        fi
     else
         IGCM_comp_modifyDefFile nonblocker orchidee.def RIVER_DESC n