Changeset 189 in lmdz_wrf for trunk

Timestamp:

Oct 16, 2014, 3:34:25 PM (10 years ago)

Author:

lfita

Message:

Adding zeroing of SNOW and LANDUSEF(==16) == 1

File:

• trunk/tools/make_met-aquaplanet.bash (modified) (7 diffs)

trunk/tools/make_met-aquaplanet.bash

@@ -19 +19 @@
   rootsh=`pwd`
   export PATH=/u/lflmd/bin/gcc_Python-2.7.5/bin:${PATH}
+  ncvarhome=${HOME}/UNSW-CCRC-WRF/tools/python
   errormsg='ERROR -- error -- ERROR -- error'
 
@@ -34 +35 @@
   rm ${ofile} > /dev/null
 
-  zerovars='LANDMASK:HGT_M:HGT_U:HGT_V:SLPX:SLPY:CON:VAR:OA1:OA2:OA3:OA4:OL1:OL2:OL3:OL4:SKINTEMP:SEAICE'
+  zerovars='LANDMASK:HGT_M:HGT_U:HGT_V:SLPX:SLPY:CON:VAR:OA1:OA2:OA3:OA4:OL1:OL2:OL3:OL4:SKINTEMP:SEAICE:SNOW'
 # Fixing p_top at 0.0001 hPa)
 ##  zTopvars='PRES:GHT:RH:VV:UU:TT'
@@ -41 +42 @@
 ##  keepvars='XLAT_M:XLONG_M:XLAT_U:XLONG_U:XLAT_V:XLONG_V:CLAT:CLONG:MAPFAC_M:MAPFAC_V:MAPFAC_U:MAPFAC_MX:MAPFAC_VX:MAPFAC_UX:MAPFAC_MY:MAPFAC_VY:MAPFAC_UY:E:F:SINALPHA:COSALPHA'
   keepvars=''
-  specificvars='LANDUSEF|1:LU_INDEX|1:SOILTEMP|0:SOILCTOP|0:SCT_DOM|14:SOILCBOT|0:SCB_DOM|14:ALBEDO12M|8:GREENFRAC|0:SNOALB|0:SLOPECAT|0'
+  specificvars='LANDUSEF|0:LU_INDEX|1:SOILTEMP|0:SOILCTOP|0:SCT_DOM|14:SOILCBOT|0:SCB_DOM|14:ALBEDO12M|8:GREENFRAC|0:SNOALB|0:SLOPECAT|0'
 #  specificvars=''
   thresvars='SST|273.15|below|275.:'
+  dimvars='LANDUSEF|z-dimension0024,15,sumc,1:'
 
   cp ${met_emfile} ${ofile}
@@ -54 +56 @@
   for var in ${zerov}; do
     echo "  "${var}"..."
-    python ${HOME}/climatescripts/python/nc_var.py -f ${ofile} -v ${var} -o valmod -S mulc,0. 
+    python ${ncvarhome}/nc_var.py -f ${ofile} -v ${var} -o valmod -S mulc,0. 
     ivar=`expr ${ivar} + 1`
 ##    exit
@@ -67 +69 @@
     varv=`echo ${var} | tr '|' ' ' | awk '{print $2}'`
     echo "  "${varn}"..."
-    python ${HOME}/climatescripts/python/nc_var.py -f ${ofile} -v ${varn} -o valmod -S mulc,0. 
-    python ${HOME}/climatescripts/python/nc_var.py -f ${ofile} -v ${varn} -o valmod -S sumc,${varv} 
+    python ${ncvarhome}/nc_var.py -f ${ofile} -v ${varn} -o valmod -S mulc,0. 
+    python ${ncvarhome}/nc_var.py -f ${ofile} -v ${varn} -o valmod -S sumc,${varv} 
     ivar=`expr ${ivar} + 1`
 ##    exit
@@ -84 +86 @@
     if test ${vark} = 'below'; then
       echo "  "${varn}" values below "${vart}" to "${varv}"..."
-      python ${HOME}/climatescripts/python/nc_var.py -f ${ofile} -v ${varn} -o valmod -S lowthres,${vart},${varv} 
+      python ${ncvarhome}/nc_var.py -f ${ofile} -v ${varn} -o valmod -S lowthres,${vart},${varv} 
     else
       echo "  "${varn}" values above "${vart}" to "${varv}"..."
-      python ${HOME}/climatescripts/python/nc_var.py -f ${ofile} -v ${varn} -o valmod -S upthres,${vart},${varv}   
+      python ${ncvarhome}/nc_var.py -f ${ofile} -v ${varn} -o valmod -S upthres,${vart},${varv}   
     fi
+    ivar=`expr ${ivar} + 1`
+##    exit
+  done # end of thresvars
+
+  # dimension values
+  ##
+  ivar=1
+  echo "dimension values..."
+  dimv=`echo ${dimvars} | tr ':' ' '`
+  for var in ${dimv}; do
+    varn=`echo ${var} | tr '|' ' ' | awk '{print $1}'`
+    varv=`echo ${var} | tr '|' ' ' | awk '{print $2}'`
+    python ${ncvarhome}/nc_var.py -f ${ofile} -v ${varn} -o valmod_dim -S ${varv}
     ivar=`expr ${ivar} + 1`
 ##    exit
@@ -117 +132 @@
     zvarv=`echo ${zvar} | tr '@' ' ' | awk '{print $2}'`
     echo "    "${zvar}
-    python ${HOME}/PY/nc_var.py -f ${ofile} -v ${zvarn} -o addvals -S '1:last:1:'${zvarv}
+    python ${ncvarhome}/nc_var.py -f ${ofile} -v ${zvarn} -o addvals -S '1:last:1:'${zvarv}
   done
-  python ${HOME}/PY/nc_var.py -f ${ofile} -o gaddattrk -S 'BOTTOM-TOP_GRID_DIMENSION|39|Inp32'
+  python ${ncvarhome}/nc_var.py -f ${ofile} -o gaddattrk -S 'BOTTOM-TOP_GRID_DIMENSION|39|Inp32'
 
 # Giving the geopotential and temperature at the top (assuming p_top 0.0001 hPa)