Changeset 3026 for trunk/LMDZ.MARS/deftank

Timestamp:

Aug 9, 2023, 11:56:36 AM (17 months ago)

Author:

jbclement

Message:

Mars PCM/PEM 1D:
Small fixes to be able to run the Mars PCM 1D without "water" + Improvements/addition of scripts in deftank/pem to run the PEM 1D model according to Laskar orbital parameters.
JBC

Location:

trunk/LMDZ.MARS/deftank

Files:

• field_def_physics_mars.xml (modified) (2 diffs)
• file_def_physics_mars.xml (modified) (4 diffs)
• pem/modify_startfi_Ls.sh (added)
• pem/modify_startfi_orbit.sh (moved) (moved from trunk/LMDZ.MARS/deftank/pem/modify_startfi.sh) (3 diffs)
• pem/run_PEM1D.sh (moved) (moved from trunk/LMDZ.MARS/deftank/pem/run_pem1d_1)

trunk/LMDZ.MARS/deftank/field_def_physics_mars.xml

@@ -334 +334 @@
                    unit="W.m-2" />
             <field id="tauTES"
+                   long_name="Water ice cloud opacity at TES abs 825 cm-1"
+                   unit="" />
+            <field id="tauTESap"
                    long_name="Water ice cloud opacity at TES abs 825 cm-1"
                    unit="" />
@@ -653 +656 @@
                    long_name="count_micro after microphysics" 
                    unit="integer" />
-            
+
+            <field id="zdqcloud_ice" 
+                   long_name="cloud ice" 
+                   unit="kg.m-2.s-1" />
+            <field id="zdqcloud_vap"
+                   long_name="cloud vap"
+                   unit="kg.m-2.s-1" />
+            <field id="zdqcloud"
+                   long_name="cloud"
+                   unit="kg.m-2.s-1" />
+
             <!-- CO2 condensation parametrization -->
             <field id="co2condens_pdtc" 

trunk/LMDZ.MARS/deftank/file_def_physics_mars.xml

@@ -12 +12 @@
               time_units="days"
               type="one_file"
-              enabled=".true.">
+              enabled=".false.">
                     
             <!-- VARS 0D -->
@@ -37 +37 @@
                 <field field_ref="u" />
                 <field field_ref="v" />
+                <field field_ref="dustq" />
+                <field field_ref="dustN" />
+                <field field_ref="h2o_vap" />
+                <field field_ref="h2o_ice" />
             </field_group>
 
@@ -83 +87 @@
               type="one_file"
               time_units="days"
-              enabled=".false.">
+              enabled=".true.">
                     
             <!-- VARS 0D -->
@@ -112 +116 @@
                 <field field_ref="h2o_layer1" />
                 <field field_ref="co2_layer1" />
+                <field field_ref="fluxsurf_dn_sw" name="ISR" />
+                <field field_ref="fluxtop_lw" name="OLR" />
             </field_group>
 

trunk/LMDZ.MARS/deftank/pem/modify_startfi_orbit.sh

@@ -4 +4 @@
 name_file="startfi.nc"
 
-new_obl=30.
-new_exc=0.
-new_Lsp=250.
-
-year_day=669.
-halfaxe=227.94
+new_obl=30.0286
+new_exc=0.08212066
+new_Lsp=347.443947402
 
 
 # Calculate modified values
+year_day=669.
+halfaxe=227.94
 pi=$(echo "4*a(1)" | bc -l)
-arg=$(echo "0.5*2.*$pi*$new_Lsp/360." | bc -l)
-zx0=$(echo "-2.*-2.0*a(s($arg)/c($arg)*sqrt((1. - $new_exc)/(1. + $new_exc)))" | bc -l)
-if [ $(echo "$zx0 < 0" | bc -l) -eq 1 ]; then
+degrad=$(echo "360./(2.*$pi)" | bc -l)
+tan=$(echo "s(0.5*$new_Lsp/$degrad)/c(0.5*$new_Lsp/$degrad)" | bc -l)
+zx0=$(echo "-2.*a($tan*sqrt((1.-$new_exc)/(1.+$new_exc)))" | bc -l)
+if [ $(echo "$zx0 <= 0" | bc -l) -eq 1 ]; then
   zx0=$(echo "$zx0 + 2.*$pi" | bc -l)
 fi
@@ -24 +24 @@
 # controle(15) = periheli ! min. Sun-Mars distance (Mkm)  ̃206.66
 # controle(16) = aphelie  ! max. SUn-Mars distance (Mkm)  ̃249.22
-# controle(16) = peri_day ! date of perihelion (sols since N. spring)
+# controle(17) = peri_day ! date of perihelion (sols since N. spring)
 # controle(18) = obliquit ! Obliquity of the planet (deg)  ̃23.98
 ncap2 -O -s "controle(17)=$new_obl" \
@@ -34 +34 @@
 # Rename the temporary file back to the original filename
 mv $name_file.temp $name_file
+
+echo New obliquit = $new_obl
+echo New periheli = $(echo "$halfaxe*(1-$new_exc)" | bc -l)
+echo New aphelie = $(echo "$halfaxe*(1+$new_exc)" | bc -l)
+echo New peri_day = $solp