Update experiments

tests/Pantanal_n48/BuildHTUs_Pantanal.pbs → tests/AmSud/BuildHTUs_SAM.pbs

@@ -3,10 +3,9 @@
 #PBS -N BuildHTUs_P
 #
 #PBS -j oe
-#PBS -l nodes=1:ppn=48
-#PBS -l walltime=128:00:00
-#PBS -l mem=180gb
-#PBS -l vmem=180gb
+#PBS -l nodes=1:ppn=64
+#PBS -l walltime=6:00:00
+#PBS -l vmem=240g,mem=200g
 #
 cd ${PBS_O_WORKDIR}
 export NSLOTS=$(($PBS_NUM_NODES*$PBS_NUM_PPN))
@@ -17,7 +16,7 @@ source ../../Environment
 #
 # Clean-up. Weights are kept for future runs.
 #
-/bin/rm -f DocumentationInterface *.nc *.txt
+/bin/rm -f DocumentationInterface *graph.nc *.txt check.out
 #
 # Run the Python code to generate the HTUs and write them into a netCDF file.
 #

tests/Iberia_n48_regular/BuildHTUs_IP.pbs → tests/AmSud/Weight.pbs

 #!/bin/bash
 #
-#PBS -N BuildHTUs_IP
+#PBS -N Weights
 #
 #PBS -j oe
 #PBS -l nodes=1:ppn=64
-#PBS -l walltime=12:00:00
+#PBS -l walltime=128:00:00
 #PBS -l mem=240gb
 #PBS -l vmem=240gb
 #
@@ -20,9 +20,8 @@ source ../../Environment
 /bin/rm -f DocumentationInterface *graph.nc *.txt
 #
 # Run the Python code to generate the HTUs and write them into a netCDF file.
-# It will also generate the weights if needed.
 #
-mpirun -n ${NSLOTS} python ../../RoutingPreProc.py
+mpirun -n ${NSLOTS} python ../../WeightsOnly.py
 if [ $? -gt 0 ] ; then
     echo "XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX"
     echo "X    Run on Iberian Peninsula failed    X"

tests/Pantanal_n48/run.def → tests/AmSud/run.def

@@ -5,9 +5,9 @@
 EarthRadius = 6370000.
 #
 ModelGridFile = /bdd/MEDI/workspaces/polcher/WRF_Forcing/AmSud-A/geo_em.d01.nc
-WEST_EAST = -66.75, -52.25
-SOUTH_NORTH = -23.5, -13.5
-HydroFile = /homedata/aschrapffer/RoutingNc_AS/SouthAmerica_FP/routing.nc
+# WEST_EAST = -90.75, -30.25
+# SOUTH_NORTH = -55.25, 25.5
+HydroFile = /homedata/aschrapffer/routingnc/AmSud/routing.nc
 #
 # FORTRAN interface parameters
 #
@@ -23,7 +23,7 @@ numop = 100
 #
 # Output
 #
-GraphFile = AmSud_A_Pantanal_graph.nc
+GraphFile = AmSud_A_graph.nc
 #
 # Diagnostics
 # You need to provide an interval in longitude and Latitude.

tests/EuroSW/BuildHTUs_EuroSW.pbs

+#!/bin/bash
+#
+#PBS -N BuildHTUs_EuroSW
+#
+#PBS -j oe
+#PBS -l nodes=2:ppn=32
+#PBS -l walltime=12:00:00
+#PBS -l vmem=500g,mem=240g
+#
+cd ${PBS_O_WORKDIR}
+export NSLOTS=$(($PBS_NUM_NODES*$PBS_NUM_PPN))
+#
+# Set the right Python 3 Anaconda environment
+#
+source ../../Environment
+#
+# Clean-up. Weights are kept for future runs.
+#
+/bin/rm -f DocumentationInterface *.txt
+#
+# If weights file does not exist, then compute it quickly
+#
+if [ ! -e Weights.nc ] ; then
+    mpirun -n ${NSLOTS} python ../../WeightsOnly.py
+    if [ $? -gt 0 ] ; then
+	exit
+    else
+	echo "======================================="
+	echo "= Weight computed                     ="
+	ls -l Weights.nc
+	echo "======================================="	
+    fi
+fi
+#
+# Run the Python code to generate the HTUs and write them into a netCDF file.
+#
+mpirun -n ${NSLOTS} python ../../RoutingPreProc.py
+if [ $? -gt 0 ] ; then
+    echo "XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX"
+    echo "X    Run on EuroSW failed               X"
+    echo "XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX"
+else
+    echo "======================================="
+    echo "= Run on EuroSW successful            ="
+    echo "======================================="
+fi
+#
+ls -l

tests/Pantanal_n48_regular/run.def → tests/EuroSW/run.def

@@ -4,10 +4,8 @@
 #
 EarthRadius = 6370000.
 #
-ModelGridFile = /bdd/ORCHIDEE_Forcing/BC/OOL/OL2/WFDEI_CRU/WFDEI_CRU_1994.nc 
-WEST_EAST = -68.75, -53.25
-SOUTH_NORTH = -24.5, -14.25
-HydroFile = /homedata/aschrapffer/RoutingNc_AS/SouthAmerica_FP/routing.nc
+ModelGridFile = /bdd/MEDI/workspaces/polcher/NewRouting/geo_EuroSW.nc
+HydroFile = /bdd/MEDI/workspaces/polcher/NewRouting/routing_MED.nc
 #
 # FORTRAN interface parameters
 #
@@ -17,13 +15,9 @@ Documentation = true
 #
 nbasmax = 35
 #
-# Number of operation of simplification performed together
-#
-numop = 50
-#
 # Output
 #
-GraphFile = WFDEI_CRU_Pantanal_test_graph.nc
+GraphFile = LC_Spain_test_graph.nc
 #
 # Diagnostics
 # You need to provide an interval in longitude and Latitude.

tests/Iberia_n48/BuildHTUs_IP.pbs → tests/Iberia/BuildHTUs_IP.pbs

@@ -17,7 +17,7 @@ source ../../Environment
 #
 # Clean-up. Weights are kept for future runs.
 #
-/bin/rm -f DocumentationInterface *graph.nc *.txt
+/bin/rm -f DocumentationInterface *graph.nc *.txt check.out
 #
 # If weights file does not exist, then compute it quickly
 #

tests/Pantanal_n48_regular/BuildHTUs_Pantanal.pbs → tests/Iberia_regular/BuildHTUs_IP.pbs

 #!/bin/bash
 #
-#PBS -N BuildHTUs_P
+#PBS -N BuildHTUs_IP
 #
 #PBS -j oe
-#PBS -l nodes=1:ppn=48
-#PBS -l walltime=128:00:00
+#PBS -l nodes=1:ppn=64
+#PBS -l walltime=12:00:00
 #PBS -l mem=240gb
 #PBS -l vmem=240gb
 #
@@ -17,9 +17,22 @@ source ../../Environment
 #
 # Clean-up. Weights are kept for future runs.
 #
-/bin/rm -f DocumentationInterface *.nc *.txt
+/bin/rm -f DocumentationInterface *graph.nc *.txt check.out
+#
+if [ ! -e Weights.nc ] ; then
+    mpirun -n ${NSLOTS} python ../../WeightsOnly.py
+    if [ $? -gt 0 ] ; then
+        exit
+    else
+        echo "======================================="
+        echo "= Weight computed                     ="
+        ls -l Weights.nc
+        echo "======================================="  
+    fi
+fi
 #
 # Run the Python code to generate the HTUs and write them into a netCDF file.
+# It will also generate the weights if needed.
 #
 mpirun -n ${NSLOTS} python ../../RoutingPreProc.py
 if [ $? -gt 0 ] ; then