Changeset 5857 for BOL/LMDZ_Setup

Timestamp:

Oct 22, 2025, 1:51:30 PM (4 weeks ago)

Author:

idelkadi

Message:

Added the possibility to use read_climoz=1 (in config.def) to read a single ozone climatology that will be used for both day and night during the simulation
To enable it, set ozone=1 in setup.sh (by default ozone=0)

Location:

BOL/LMDZ_Setup

Files:

• script_SIMU (modified) (2 diffs)
• setup.sh (modified) (9 diffs)

BOL/LMDZ_Setup/script_SIMU

@@ -64 +64 @@
 ok_guide=y   # y/n guidage ou non
 climato=1
+ozone=0
 
 echo '##############################################################'
@@ -197 +198 @@
 if [ $climato = 1 ] ; then
    ${GET} $LIMITDIR/limit.nc limit.nc
+   if [ $ozone = 1 ] ; then ${GET} $LIMITDIR/climoz_LMDZ_1979_2008_clim.nc climoz_LMDZ.nc ; fi
 else
    ${GET} $LIMITDIR/limit.$year.nc limit.nc
+   if [ $ozone = 1 ] ; then ${GET} $LIMITDIR/climoz_$year.nc climoz_LMDZ.nc ; fi
 fi
 

BOL/LMDZ_Setup/setup.sh

@@ -35 +35 @@
 
   aerosols="clim"
+
+  # OZONE : 
+  # AI et ZE oct 2025
+  #         0 => analytic formula fonction of altitude and latitude
+  #         1 => ozonread a single ozone climatology that will be used day and night
+  ozone=0
 
   # SURFACE/VEGETATION SCHEME
@@ -42 +48 @@
   #   If you need other orch versions, and also require XIOS, you'll need
   #                    to create the appropriate files in DEF/XMLfilesOR...
-  veget="CMIP6"
+  veget="none"
 
   # New snow scheme INLANDSIS! "y" / "n"
@@ -139 +145 @@
              -netcdf  0/1/DIR        Netcdf installation
              -aerosols n/clim/spla   Aerosols : none / climatological / interactive with SPLA
+             -ozone 0/1              Ozone : not read / read an ozone climatology
              -veget CMIP6            CMIP6 version of rchidee
                     7983             orch2.2 
@@ -169 +176 @@
       "-nudging") ok_guide=y; shift;;
       "-climato") climato=$2; shift 2;;
+      "-ozone") ozone=$2; shift 2;;
       "-mthini") mthini=$2; shift 2;;
       "-mthend") mthend=$2; shift 2;;
@@ -390 +398 @@
   pwd
 
+  # OZONE : modifier la cles logique pour ozone
+  if [[ $ozone -eq 1 ]]; then
+     sed -i'' -e 's/read_climoz=.*.$/read_climoz='$ozone'/' DEF/config.def
+  fi
+
   # COSP : ok_cosp desactive COSP si on a compile avec; il ne l'active pas si on a compile sans
   sed -i'' -e 's/ok_cosp.*.$/ok_cosp='$cosp'/' DEF/config.def
@@ -533 +546 @@
   -e "s/^isotopes=.*.$/isotopes=$isotopes/" \
   -e "s/^climato=.*.$/climato=$climato/" \
+  -e "s/^ozone=.*.$/ozone=$ozone/" \
   -e "s/^ok_guide=.*.$/ok_guide=$ok_guide/" \
   "$local_d/script_SIMU" >| "$SIMRUNDIR/tmp_$SIM"
@@ -664 +678 @@
        ln_from_pub  3DInputData/Limit $file
     done
+
+    #  Ozone
+    if [[ $ozone -eq 1 ]]; then
+       if [[ $climato -eq 1 ]]; then
+        echo LISTEGET tro3_1979_2008_clim.nc
+           wget_pub 3DInputData/Ozone/original tro3_1979_2008_clim.nc
+           ln_from_pub  3DInputData/Ozone/original tro3_1979_2008_clim.nc
+       elif [[ $climato -eq 0 ]]; then
+           yr=$yrini
+           #while [[ $yr -le $yrend ]]; do
+           echo LISTEGET tro3_${yr}.nc
+           wget_pub 3DInputData/Ozone/original tro3_$(( $yr - 1 )).nc
+           ln_from_pub 3DInputData/Ozone/original tro3_$(( $yr - 1 )).nc
+           wget_pub 3DInputData/Ozone/original tro3_${yr}.nc
+           ln_from_pub 3DInputData/Ozone/original tro3_$yr.nc
+           wget_pub 3DInputData/Ozone/original tro3_$(( $yr + 1 )).nc
+           ln_from_pub 3DInputData/Ozone/original tro3_$(( $yr + 1 )).nc
+           #done
+       fi
+    fi
     #-------------------------------------------------------------------
     # Initial state ECDYN
@@ -717 +751 @@
 echo "Executable : $ce0l"
 for yr in $yrs; do
+  # sst and sic
   if [ $climato = 0 ]; then sufyr=$suf\$yr; else sufyr=$suf; fi
   ln -sf amipbc_sic_\$sufyr.nc amipbc_sic_1x1.nc
   ln -sf amipbc_sst_\$sufyr.nc amipbc_sst_1x1.nc
+  # ozone
+  if [[ $ozone -eq 1 ]]; then
+     if [[ $climato -eq 1 ]]; then
+        ln -s tro3_1979_2008_clim.nc climoz.nc
+     else
+        ln -s tro3_${yr}.nc  climoz.nc
+        ln -s tro3_$(( $yr - 1 )).nc climoz_m.nc
+        ln -s tro3_$(( $yr + 1 )).nc climoz_p.nc
+     fi
+   fi
+  # modif annee dans run.def
   sed -e 's/anneeref=.*.$/anneeref='\$yr'/' DEF/run.def >| run.def
   echo Starting initialisation
   # Runing ce0l.e
   $mpicmd $ce0l  # ce0l requires MPI=OMP=1
+  # Renomer limit et/ou climoz_LMDZ.nc si climato=0
   if [ $climato = 0 ]; then mv limit.nc limit.\$yr.nc; fi
+  if [[ $ozone -eq 1 ]]; then
+     if [ $climato = 0 ]; then mv climoz_LMDZ.nc climoz_LMDZ_\$yr.nc; fi
+  fi
 done
 # ANCIEN MULTI STEP   ;;
@@ -740 +790 @@
 for f in sta* gri*nc; do cp \$f $local_d/$INIT/\$f; done
 if [[ $climato = 1 ]]; then cp limit.nc $local_d/$INIT/limit.nc; fi
+if [[ $climato = 1 ]]; then
+   if [[ $ozone -eq 1 ]]; then cp climoz_LMDZ.nc $local_d/$INIT/climoz_LMDZ_1979_2008_clim.nc; fi
+fi
 mkdir -p $local_d/$LIMIT
 for f in limit*.nc ; do cp \$f $local_d/$LIMIT/\$f; done
+# ozone
+if [[ $ozone -eq 1 ]]; then
+   if [[ $climato -eq 1 ]]; then cp climoz_LMDZ.nc $local_d/$LIMIT/climoz_LMDZ_1979_2008_clim.nc; fi
+   if [[ $climato -eq 0 ]]; then cp climoz_LMDZ_${yr}.nc $local_d/$LIMIT/.; fi
+fi
 if [ $aerosols = clim ]; then  for f in aerosols[.0-9]*nc; do cp \$f $local_d/$LIMIT/\$f; done; fi
 #