Changeset 3026 for trunk/LMDZ.MARS/deftank/pem/modify_startfi_orbit.sh

Timestamp:

Aug 9, 2023, 11:56:36 AM (16 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

File:

• trunk/LMDZ.MARS/deftank/pem/modify_startfi_orbit.sh (moved) (moved from trunk/LMDZ.MARS/deftank/pem/modify_startfi.sh) (3 diffs)

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