Lancement des simulations 1D
mkdir 1D
cd 1D
cp ~/LMDZ/phyparam/def_files/* .
cp ~/LMDZ/phyparam/utils/run_eqrc.sh .
cp ~/LMDZ/phyparam/utils/*py .
cp ../scm.e
run_eqrc.sh


Figure1D="fig-vol3-chp-gcm-zprofphys0theta.pdf fig-vol3-chp-gcm-zprofphys4theta.pdf fig-vol3-chp-gcm-zprofphys0temp.pdf fig-vol3-chp-gcm-zprofphys4temp.pdf fig-vol3-chp-gcm-evol0.pdf fig-vol3-chp-gcm-wth.pdf"
Eqrad="fig-vol3-chp-gcm-EqRad.pdf"
Evol="fig-vol3-chp-gcm-evol.pdf"
Flux="fig-vol3-chp-gcm-flux0.pdf fig-vol3-chp-gcm-flux1.pdf"

\cp -f $Figure1D $Eqrad $Evol $Flux ~/TEX/LIVRES/ModelPlanet/tex/fig

-------------------------------------------------------------------------------
Lancement des simulations 1D pour le livre
-------------------------------------------------------------------------------
cd PARAM20230218/modipsl/modeles/LMDZ
1D
cp ~/LMDZ/phyparam/def_files/* 1D/
mv scm.e 1D
cd 1D
cp ~/LMDZ/phyparam/utils/run_eqrc.sh .

-------------------------------------------------------------------------------
Inversion des latitudes entre phys et dynzon ...

ncrcat -v ts,swtopdn,swtopup,lwtop SIM[1-5]/phys.nc -O srf.nc



simu=mars
case $simu in
   mars) dir=ghio:/data/hourdin/phyparam/PARAM20230114.trunk/modipsl/modeles/LMDZ/MARS0/SIM4 ;;
   terre)  dir=ghio:/data/hourdin/phyparam/PARAM20230114.trunk/modipsl/modeles/LMDZ/TERRE0/SIM4 ;;
   *) exit 1
esac

if [ ! -f phys_${simu}.nc ] ; then scp $dir/phys.nc phys_${simu}.nc ; fi
if [ ! -f zon_${simu}.nc ; then
   scp $dir/dynzon.nc zon_${simu}.nc ; fi
   ncatted -a positive,presnivs,o,c,down zon_${simu}.nc 
   ncrename -v T,temp -O zon_${simu}.nc
fi



Calendrier partant de l'equinoxe de printemps, 22 mars, jour 80 de l'année
juin - aout : 151->240  -80   70 -> 160 
dec  - fev  :                250 -> 340

ncks -d time,250,340 phys.nc -O ete.nc
cdo timavg -zonavg ete.nc ../zete.nc
ncks -d time,310 phys.nc -O ../aout1.nc


ncatted -a positive,presnivs,o,c,down dynzon.nc 
ncrename -v t,temp -O dynzon.nc
ncks -d time,25,34 phys.nc -O ete.nc
cdo timavg -zonavg ete.nc ../zete.nc
ncks -d time,31 phys.nc -O ../aout1.nc

-------------------------------------------------------
Gestion git
mkdir git_param
cd git_param/
git init
git add dynphy_lonlat/phyparam phyparam/
git commit -m "physique a 20 parametre"
-------------------------------------------------

Nouvelle version de phyparam

sous git.
Le repertoire de ref est sur 
phyparam

On cree aussi des archives par securité avec le meme nom que la svn avec laquelle le code à tourné.

==============================================================
Tendances et thermique
======================
On vise l'hiver de l'hémispère nord
Les simulations commençant le 1er mars, il s'agit des pas de temps
26 à 34
Refaire tourner un an de simulation
Puis
ncks -d time,25:33 phys.nc -O tmp.nc
cdo zonavg -timavg tmp.nc zon_terre.nc
Puis refaire tourner 290 jours avant de lancer 10 jours avec de la haute fréquence.
ncks -d lon,0 phys.nc phys0.nc
==============================================================

Il faut faire une simulation où on stoke les tendances tous les 10 jours
puis une moyenne sur les pas de temps

On refait tourner les deux planètes sur 10 jours avec une sortie par heure.

MARS0_4 et TERRE0C_16


