Changeset 661 in lmdz_wrf for trunk

Timestamp:

Dec 9, 2015, 12:43:18 PM (9 years ago)

Author:

lfita

Message:

Adding las improvements in the script

File:

• trunk/tools/getting_fields-trajectories.bash (modified) (6 diffs)

trunk/tools/getting_fields-trajectories.bash

@@ -16 +16 @@
 # tevol_boxtraj.nc: file with the values following the trajectory
 
-models='wrf'
+models='wrf:lwrf'
+
+ifold=/media/ExtDiskD/bkup/ciclad/etudes/WL_HyMeX/superstorm
+ofold=/home/lluis/etudes/WL_HyMeX/superstorm
 
 #WRFexps='control:mp1:pbl1:pbl2:rad1'
-WRFexps='control:pbl2:rad1'
-LMDZexps='lmdza:lmdzb:lmdzora:lmdzorb'
+WRFexps='control:cu1:pblsfc1:rad1:shallow'
+LMDZexps='lmdza:lmdzb'
 LWRFexps='lmdza:lmdzb'
 
-
 # Domain to analyse
-dom='2'
-
-ofold=/homedata/lfita/etudes/WRFdynamicoSENS
-HOMEpy=/home/lfita/etudes/WRF_LMDZ/svn/LMDZ_WRF/tools
-
-use_concatenated=false
+dom='1'
+
+#ofold=/homedata/lfita/etudes/WRFdynamicoSENS
+#HOMEpy=/home/lfita/etudes/WRF_LMDZ/svn/LMDZ_WRF/tools
+HOMEpy=/home/lluis/etudes/WRF_LMDZ/svn/LMDZ_WRF/tools
+
+use_concatenated=true
 
 ####### ###### ##### #### ### ## #
@@ -37 +40 @@
 
 #WRFvars='all'
-WRFvars='Times:XLONG:XLAT:T2,U10,V10,P,PB,PSFC,HGT,T,QVAPOR,RAINC,RAINNC'
+WRFvars='Times,XLONG,XLAT,T2,U10,V10,P,PB,PSFC,HGT,T,QVAPOR,RAINC,RAINNC'
 
 #LMDZvars='all'
-LMDZars='Time_counter:longitude:latitude:u10m:v10m:t2m,wind10m,slp,precip'
-
-LWRFvars='Times:XLONG:XLAT:T2,U10,V10,P,PB,PSFC,HGT,T,QVAPOR,RAINC,RAINNC'
+LMDZars='Time_counter,longitude,latitude,u10m,v10m,t2m,wind10m,slp,precip'
+
+LWRFvars='Times,XLONG,XLAT,T2,U10,V10,P,PB,PSFC,HGT,T,QVAPOR,RAINC,RAINNC'
 #LWRFvars='LU_INDEX,ZNU,ZNW,ZS,DZS,U,RU_TEND,V,RV_TEND,W,PH,PHB,T,HFX_FORCE,LH_FORCE,'
 #LWRFvars=${LWRFvars}'TSK_FORCE,HFX_FORCE_TEND,LH_FORCE_TEND,TSK_FORCE_TEND,MU,MUB,'
@@ -117 +120 @@
 LMDZtrjvars='time_counter,u10m,v10m,precip,slp,psol'
 #LMDZ lon/lat file
-if test ! -f lmdz_lonlat.nc; then
-  python ${HOMEpy}/nc_var.py -o selvar -v t2m -f \
-    lmdz/run/lmdza/lmdzout/histins.nc -S 'lon@lon,lat@lat'
-  mv selvar_new.nc lmdz_lonlat.nc
-fi
+#if test ! -f lmdz_lonlat.nc; then
+#  python ${HOMEpy}/nc_var.py -o selvar -v t2m -f \
+#    ${ifold}/lmdz/run/lmdza/lmdzout/histins.nc -S 'lon@lon,lat@lat'
+#  mv selvar_new.nc lmdz_lonlat.nc
+#fi
 
 WRFvdims='XLONG,XLAT,ZNU,time'
@@ -167 +170 @@
     osim=${ofold}/${exp}
     if test ${mod} = 'wrf' || test ${mod} = 'lwrf'; then
-      isim=${osim}/wrfout
+      isim=${ifold}/${mod}/run/${exp}/wrfout
 # Trajectory file
 ##
@@ -235 +238 @@
           for file in ${isim}/wrfout_d0${dom}*; do
             fdate=`echo ${file} | tr '_' ' ' | awk '{print $3"_"$4}'`
-            echo "  "${file}
-            python ${HOMEpy}/diagnostics.py -f ${osim}/full_concatenated.nc              \
-              -d 'Time@Times,bottom_top@ZNU,south_north@XLAT,west_east@XLONG' -v ${diags}\
-              -f ${file}
-            #exit
-            if test $? -ne 0; then
-              echo ${errormsg}
-              echo "  "${main}": python does not work !!"
-              echo "  python "${HOMEpy}"/diagnostics.py -f "${osim}"/full_concatenated.nc \
-                -d 'Time@Times,bottom_top@ZNU,south_north@XLAT,west_east@XLONG' -v "${diags} \
-                "-f "${file}
-              exit
-            else
-              mv diagnostics.nc ${osim}/diagnostics_${fdate}.nc
+            if test ! ${osim}/diagnostics_${fdate}.nc; then
+              echo "  "${file}
+              python ${HOMEpy}/diagnostics.py -f ${osim}/full_concatenated.nc              \
+                -d 'Time@Times,bottom_top@ZNU,south_north@XLAT,west_east@XLONG' -v ${diags}\
+                -f ${file}
+              #exit
+              if test $? -ne 0; then
+                echo ${errormsg}
+                echo "  "${main}": python does not work !!"
+                echo "  python "${HOMEpy}"/diagnostics.py -f "${osim}"/full_concatenated.nc \
+                  -d 'Time@Times,bottom_top@ZNU,south_north@XLAT,west_east@XLONG' -v "${diags} \
+                  "-f "${file}
+                exit
+              else
+                mv diagnostics.nc ${osim}/diagnostics_${fdate}.nc
+              fi
             fi
           done
@@ -292 +297 @@
 # Following trajectory values
 ##
-    if test ! -f ${osim}/tevol_boxtraj.nc && -f ${osim}/trajectory.dat; then
+    if test ! -f ${osim}/tevol_boxtraj.nc && test -f ${osim}/trajectory.dat; then
       tfirst=`cat ${osim}/trajectory.dat | grep -v '#' | head -n 1 | awk '{print $1-1}'`
       echo "tfirst: "${tfirst}