Changeset 2325 in lmdz_wrf for trunk/tools

Timestamp:

Feb 8, 2019, 3:43:19 PM (6 years ago)

Author:

lfita

Message:

Fixing in `compute_slices_stats_areaweighted':

• right coordinates on gridarea variable
• right conputatin of gridarea

Provide right example

Location:

trunk/tools

Files:

• nc_var.py (modified) (1 diff)
• nc_var_tools.py (modified) (3 diffs)

trunk/tools/nc_var.py

@@ -70 +70 @@
 ## e.g. # nc_var.py -o area_weighted -f '/media/lluis/ExtDiskC_ext4/bkup/llamp/estudios/dominios/SA150k/geo_em.d01.nc:west_east;XLONG_M;WRFxbnds;-1;south_north;XLAT_M;WRFybnds;-1,/media/lluis/ExtDiskC_ext4/bkup/llamp/estudios/dominios/SA50k/geo_em.d01.nc:west_east;XLONG_M;WRFxbnds;-1;south_north;XLAT_M;WRFybnds;-1' -S 'no:mean' -v HGT_M
 ## e.g. # nc_var.py -o compute_slices_stats_areaweighted -S 'XLONG,-74.,-36.4,4.;XLAT,-63.,19.,4.;HGT,500.,7000.,500.@Time|WRFtime:west_east|XLONG:south_north|XLAT@Time@west_east|lon_bnds,south_north|lat_bnds@XLONG|lat_bnds;lon_bnds,XLAT|lat_bnds;lon_bnds@Time' -f wrfout_bnds.nc -v T2,U10,V10,Q2
+## e.g. # nc_var.py -o compute_slices_stats_areaweighted -S 'lat,-63.,19.,2.;orog,500.,7000.,500.;rangefaces,fixed,-2.5|-0.5,-0.5|0.5,0.5|2.5@time|time:lon|lon:lat|lat@time@lon|lon_bnds,lat|lat_bnds@lon|lon_bnds,lat|lat_bnds@lat,lon@time' -f /media/lluis/ExtDiskC_ext3/DATA/estudios/Andes/DATA/concatenated/historical/tasmin/tasmin_Amon_ACCESS1-0_historical_r1i1p1_185001-200512_Andes_19600101000000-19900101000000.nc -v tasmin
 
 from optparse import OptionParser

trunk/tools/nc_var_tools.py

@@ -28633 +28633 @@
             basicvardef(anewvar ,dn+'gridarea', "area of the grids cells from " +    \
              dn, vu)
-            anewvar.setncattr('coordinates',' '.join(varfinaldims[::-1]))
+            anewvar.setncattr('coordinates',' '.join(Spgrid[::-1]))
             anewvar[:] = areas[0:Ngridsinmax,...]
 
@@ -28837 +28837 @@
                         ixB = sliceinB[1,inpointsB[iv,jB,iB,jA,iA],jB,iB]
                         iyB = sliceinB[0,inpointsB[iv,jB,iB,jA,iA],jB,iB]
-                        anewvar[iv,jB,iB,jA,iA]= (aA*pA)*(aB*pB)
-                        slicearea = slicearea + (aA*pA)*(aB*pB)
+                        anewvar[iv,jB,iB,jA,iA]= aA*pA*pB
+                        slicearea = slicearea + aA*pA*pB
                     if np.isnan(slicearea):
                         print errormsg
@@ -29002 +29002 @@
                                     ixC = sliceinC[1,ivC,jC,iC]
                                     iyC = sliceinC[0,ivC,jC,iC]
-                                    anewvar[iv,jC,iC,jB,iB,jA,iA]= (aA*pA)*(aB*pB)
+                                    anewvar[iv,jC,iC,jB,iB,jA,iA]= aA*pA*pB
                                     pnewvar[iv,jC,iC,jB,iB,jA,iA]= pA*pB
-                                    slicearea = slicearea + (aA*pA)*(aB*pB)
+                                    slicearea = slicearea + aA*pA*pB
                                 if np.isnan(slicearea):
                                     print errormsg