Changeset 4346

Timestamp:

Nov 10, 2022, 4:18:58 PM (23 months ago)

Author:

fhourdin

Message:

New version with -nconly option

Allows to output all the variables of a routine in the phys.nc
file, even with not replaying this particular routine.

Location:

BOL/Replay

Files:

• replay_1D.sh (modified) (1 diff)
• replay_clean.sh (modified) (1 diff)
• replay_equip.sh (modified) (15 diffs)

BOL/Replay/replay_1D.sh

@@ -20 +20 @@
 replay_clean.sh
 replay_equip.sh $routine
+# If wanting to equip other routines with nc outputs only
+# replay_equip.sh -nconly -pre TOTO other_routine
 cd -
 

BOL/Replay/replay_clean.sh

@@ -2 +2 @@
 
 replay_comment="replay automatic include"
-\rm -f dump_param*.h dump_ini.h call_param_replay.F90 dump_ini_module.F90 get_ini_module.F90 input tmp tmp2
+\rm -f dump_replay*.h mcall_param_replay.F90 dump_ini_module.F90 get_ini_module.F90 input tmp tmp2
 for file in `grep "$replay_comment" *F90 | cut -d: -f1` ; do
    sed -i"" -e "/$replay_comment/d" $file

BOL/Replay/replay_equip.sh

@@ -74 +74 @@
 #=============================================================================
 
-param=`basename $1 .F90`
+#-----------------------------------------------------------------------------
+# Reading rguments
+#-----------------------------------------------------------------------------
+nconly=false
+if [ $# = 0 ] ; then $0 -h ; exit ; fi
+while (($# > 0))
+   do
+   case $1 in
+     -h|--help) cat <<........fin
+           $0 [-nconly] [-pre prefix] routine_name
+           The prefix will be put on each variable stored in the nc file
+........fin
+        exit ;;
+     -nconly) nconly=true ; shift ;;
+     -pre) prefix=$2 ; shift ; shift ;;
+     *) if (( $# > 1 )) ; then $0 -h ; exit ; fi ; param=`basename $1 .F90` ; shift
+    esac
+done
+if [ "$param" = "" ] ; then $0 -h ; exit ; fi
+
 case $param in
    thermcell_main|thermcell_plume_6A|thermcell_env|thermcell_height|thermcell_dry|\
@@ -84 +103 @@
 
 replay_comment="replay automatic include"
-paraminc1=dump_param1.h
-paraminc2=dump_param2.h
-iniinc=dump_ini.h
+paraminc1=dump_replay_${param}_head.h
+paraminc2=dump_replay_${param}_nc_${prefix}.h
+iniinc=dump_replay_ini.h
+paramfile=`grep -i "subro.* ${param}[\ (]" *.F90 | sed -e 's/ //g' | grep "${param}(" | cut -d: -f1`
+echo La parametrisation $param est contenue dans le fichier $paramfile
+if [ `echo ${paramfile} | wc -w` != 1 ] ; then echo file $paramfile multiple  ; $0 -h ; exit ; fi
+
 
 #-----------------------------------------------------------------------------
@@ -118 +141 @@
 echo '! dump_param_open >>> '
 }
+
 
 #-----------------------------------------------------------------------------
@@ -133 +157 @@
    \rm -f tmp tmp2
 }
+
 
 #-----------------------------------------------------------------------------
@@ -159 +184 @@
 }
 
-#-----------------------------------------------------------------------------
-function get_param_file(){
-#-----------------------------------------------------------------------------
-  p=$1
-  grep -i "subro.* $p[\ (]" *.F90 | sed -e 's/ //g' | grep "$p(" | cut -d: -f1
-}
+
+
 
 #-----------------------------------------------------------------------------
@@ -181 +202 @@
 
 #-----------------------------------------------------------------------------
+function iotd_calls(){
+#-----------------------------------------------------------------------------
+    klev=$1 ; shift
+    pre=$1 ; shift
+    if (( $# >> 0 )) ; then
+       vars=$*
+       for var in $vars ; do
+          echo "call iotd_ecrit_seq('"$pre$var"',1,'"$pre$var in $param"',' ',"$var")"
+       done
+    fi
+}
+         
+
+#-----------------------------------------------------------------------------
 # On nettoye les inclusions précédente dans les fichiers .F90
 #-----------------------------------------------------------------------------
 
-for file in `grep "$replay_comment" *.F90 2> /dev/null | cut -d: -f1` ; do
-   sed -i"" -e "/$replay_comment/d" $file
-done
-line=`sed -n -e "/CALL wake_ini/=" physiq_mod.F90 | head -1`
-sed -i"" -e "${line}s/$/\n   CALL iophys_ini(pdtphys) ! $replay_comment  ! $replay_comment/" physiq_mod.F90
-
-#-----------------------------------------------------------------------------
-# Modifying the parameterization routine
-#-----------------------------------------------------------------------------
-
-paramfile=`get_param_file $param`
-echo La parametrisation $param est contenue dans le fichier $paramfile
-if [ `echo ${paramfile} | wc -w` != 1 ] ; then echo file $paramfile multiple  exit ; fi
-extract_subroutine $param $paramfile
+if [ $nconly = false ] ; then
+   for file in `grep "$replay_comment" *.F90 2> /dev/null | cut -d: -f1` ; do
+      sed -i"" -e "/$replay_comment/d" $file
+   done
+   line=`sed -n -e "/CALL wake_ini/=" physiq_mod.F90 | head -1`
+   sed -i"" -e "${line}s/$/\n   CALL iophys_ini(pdtphys) ! $replay_comment  ! $replay_comment/" physiq_mod.F90
+fi
 
 #-----------------------------------------------------------------------------
@@ -203 +231 @@
 #-----------------------------------------------------------------------------
 
+extract_subroutine $param $paramfile # -> input file
 varin0=`grep inten.*.in input | sed -e 's/\!.*$//' | cut -d: -f3 | sed -e 's/,/ /g'`
 varinout0=`grep inten.*.inout input | sed -e '/\!.*$/d' | cut -d: -f3 | sed -e 's/,/ /g'`
@@ -213 +242 @@
    full) search_str='real'
 esac
+var_1D_inout=`grep -i "$search_str" input | grep intent.*inout | grep $klon | sed -e 's/!.*$//' -e /$klev/d | cut -d: -f3 | sed -e 's/,/ /g'`
+var_2D_inout=`grep -i "$search_str" input | grep intent.*inout | grep $klon | grep $klev | cut -d: -f3 | sed -e 's/!.*$//' -e 's/,/ /g'`
 var_1D_noarg=`grep -i "$search_str" input | sed -e /intent/d | grep $klon | sed -e 's/!.*$//' -e /$klev/d | cut -d: -f3 | sed -e 's/,/ /g'`
 var_2D_noarg=`grep -i "$search_str" input | sed -e /intent/d | grep $klon | grep $klev | cut -d: -f3 | sed -e 's/!.*$//' -e 's/,/ /g'`
@@ -218 +249 @@
 var_2D=`grep -i "$search_str" input | grep $klon | grep $klev | cut -d: -f3 | sed -e 's/!.*$//' -e 's/,/ /g'`
 echo varin  : $varin
-echo var_1D_noarg : $var_1D_noarg
-echo var_2D_noarg : $var_2D_noarg
-echo var_1D : $var_1D
-echo var_2D : $var_2D
+echo var_1D_inout : $var_1D_inout ; echo
+echo var_2D_inout : $var_2D_inout ; echo
+echo var_1D_noarg : $var_1D_noarg ; echo
+echo var_2D_noarg : $var_2D_noarg ; echo
+echo var_1D : $var_1D ; echo
+echo var_2D : $var_2D ; echo
 
 #-----------------------------------------------------------------------------
 # Ecriture de la routine d'appel a la parametrisation en mode replay
 #-----------------------------------------------------------------------------
+
+if [ $nconly = false ] ; then
 
 cat > call_param_replay.F90 <<eod
@@ -268 +303 @@
 eod
 
-
-#-----------------------------------------------------------------------------
-# Creation d'un .h d'ecriture de toutes les variables intent in & out
+fi # nconly = false
+
+#-----------------------------------------------------------------------------
+# Creating file dump_${param}_head.h
 #-----------------------------------------------------------------------------
 
 \rm $paraminc1
+if [ $nconly = false ] ; then
 cat> $paraminc1 <<eod
 logical, save :: first=.true.
@@ -286 +323 @@
 write(81,rec=irec) $varin
 eod
+iotd_calls 1     in_${prefix} $var_1D_inout >> $paraminc1
+iotd_calls $klev in_${prefix} $var_2D_inout >> $paraminc1
+fi # nconly = false
+
 for var in $var_1D_noarg $var_2D_noarg ; do echo $var=0. >> $paraminc1 ; done
 include_line "include \"$paraminc1\"" $param $paramfile after_declarations
 
+#-----------------------------------------------------------------------------
+# Creating file dump_${param}_nc_${prefix}.h
+#-----------------------------------------------------------------------------
+
 \rm $paraminc2
-for var in $var_1D ; do
-echo "call iotd_ecrit_seq('"$var"',1,'"$var"',' ',"$var")" >> $paraminc2
-done
-for var in $var_2D ; do
-echo "call iotd_ecrit_seq('"$var"',"$klev",'"$var"',' ',"$var")" >> $paraminc2
-done
+iotd_calls 1     "${prefix}" $var_1D >> $paraminc2
+iotd_calls $klev "${prefix}" $var_2D >> $paraminc2
 include_line "include \"$paraminc2\"" $param $paramfile before_return
   
@@ -303 +344 @@
 #-----------------------------------------------------------------------------
 
+if [ $nconly = false ] ; then
 varinmod=`grep -i 'intent.in' $inimod.F90 | sed -e 's/\!.*$//' | cut -d: -f3 | sed -e 's/,/ /g'`
 varinmodv=`echo $varinmod | sed -e 's/ /,/g'`
@@ -311 +353 @@
 eod
 include_line "include \"$iniinc\"" $paramini $inimod.F90  before_return
+fi # nconly = false
 
 
@@ -316 +359 @@
 # get_ini_module gere l'initialisation du module en mode replay
 #-----------------------------------------------------------------------------
-
+if [ $nconly = false ] ; then
 cat > get_ini_module.F90 <<eod
 subroutine get_ini_module
@@ -334 +377 @@
 end
 eod
+fi # nconly = false
 
 #-----------------------------------------------------------------------------