[Ncep.list.emc.fv3gfs_tickets] FV3GFS Ticket #14: Modify global_chgres to output NSST data in NetCDF format
FV3GFS Trac Ticket
ncep.list.emc.fv3gfs_tickets at noaa.gov
Wed Apr 12 19:39:48 UTC 2017
#14: Modify global_chgres to output NSST data in NetCDF format
-----------------------------+----------------------------
Reporter: george.gayno@… | Owner: george.gayno@…
Type: task | Status: accepted
Priority: major | Milestone:
Component: component1 | Version:
Resolution: | Keywords: chgres
-----------------------------+----------------------------
Comment (by george.gayno@…):
Now test changes at r91211 on Luna using the
[https://svnemc.ncep.noaa.gov/trac/fv3gfs/browser/branches/chgres_nsst_netcdf/run_chgres_nemsio_cray.sh?rev=91211
run_chgres_nemsio_cray.sh] test script. The test converts both a surface
and NSST file. The test ran to completion and produced these files:
{{{
/gpfs/hps/stmp/George.Gayno/chgres_save.sfc.nsst\ $ ls -l
total 871680
-rw-r--r-- 1 George.Gayno g01 148674296 Apr 12 19:18 out.sfc.tile1.nc
-rw-r--r-- 1 George.Gayno g01 148674296 Apr 12 19:18 out.sfc.tile2.nc
-rw-r--r-- 1 George.Gayno g01 148674296 Apr 12 19:19 out.sfc.tile3.nc
-rw-r--r-- 1 George.Gayno g01 148674296 Apr 12 19:19 out.sfc.tile4.nc
-rw-r--r-- 1 George.Gayno g01 148674296 Apr 12 19:19 out.sfc.tile5.nc
-rw-r--r-- 1 George.Gayno g01 148674296 Apr 12 19:20 out.sfc.tile6.nc
}}}
Checking the header of one of the files showed both surface and NSST
records:
{{{
/gpfs/hps/stmp/George.Gayno/chgres_save.sfc.nsst\ $ ncdump -h
out.sfc.tile3.nc
netcdf out.sfc.tile3 {
dimensions:
lon = 768 ;
lat = 768 ;
lsoil = 4 ;
variables:
float lon(lon) ;
lon:cartesian_axis = "X" ;
float lat(lat) ;
lat:cartesian_axis = "Y" ;
float lsoil(lsoil) ;
lsoil:cartesian_axis = "Z" ;
float geolon(lat, lon) ;
float geolat(lat, lon) ;
float slmsk(lat, lon) ;
float tsea(lat, lon) ;
float sheleg(lat, lon) ;
float tg3(lat, lon) ;
float zorl(lat, lon) ;
float alvsf(lat, lon) ;
float alvwf(lat, lon) ;
float alnsf(lat, lon) ;
float alnwf(lat, lon) ;
float vfrac(lat, lon) ;
float canopy(lat, lon) ;
float f10m(lat, lon) ;
float t2m(lat, lon) ;
float q2m(lat, lon) ;
float vtype(lat, lon) ;
float stype(lat, lon) ;
float facsf(lat, lon) ;
float facwf(lat, lon) ;
float uustar(lat, lon) ;
float ffmm(lat, lon) ;
float ffhh(lat, lon) ;
float hice(lat, lon) ;
float fice(lat, lon) ;
float tisfc(lat, lon) ;
float tprcp(lat, lon) ;
float srflag(lat, lon) ;
float snwdph(lat, lon) ;
float shdmin(lat, lon) ;
float shdmax(lat, lon) ;
float slope(lat, lon) ;
float snoalb(lat, lon) ;
float stc(lsoil, lat, lon) ;
float smc(lsoil, lat, lon) ;
float slc(lsoil, lat, lon) ;
float tref(lat, lon) ;
float z_c(lat, lon) ;
float c_0(lat, lon) ;
float c_d(lat, lon) ;
float w_0(lat, lon) ;
float w_d(lat, lon) ;
float xt(lat, lon) ;
float xs(lat, lon) ;
float xu(lat, lon) ;
float xv(lat, lon) ;
float xz(lat, lon) ;
float zm(lat, lon) ;
float xtts(lat, lon) ;
float xzts(lat, lon) ;
float d_conv(lat, lon) ;
float ifd(lat, lon) ;
float dt_cool(lat, lon) ;
float qrain(lat, lon) ;
}
}}}
Users may also choose to convert a surface file only. This option was
then checked. The test script was altered as follows:
{{{
/gpfs/hps/emc/global/noscrub/George.Gayno/fv3gfs/branches/chgres_nsst_netcdf\
$ svn diff run_chgres_nemsio_cray.sh
Index: run_chgres_nemsio_cray.sh
===================================================================
--- run_chgres_nemsio_cray.sh (revision 91212)
+++ run_chgres_nemsio_cray.sh (working copy)
@@ -108,7 +108,7 @@
export
CHGRESVARS="use_ufo=.false.,nst_anl=.false.,idvc=2,idvt=21,idsl=1,IDVM=0,nopdpvv=.true."
export SIGINP=NULL
export SFCINP=$COMROT/sfnanl.gfs.$CDATE
-export NSTINP=$COMROT/nsnanl.gfs.$CDATE
+export NSTINP=NULL
}}}
Checking the output files showed surface records only:
{{{
/gpfs/hps/stmp/George.Gayno/chgres_save.sfc\ $ ls -l
total 622848
-rw-r--r-- 1 George.Gayno g01 106197450 Apr 12 19:24 out.sfc.tile1.nc
-rw-r--r-- 1 George.Gayno g01 106197450 Apr 12 19:24 out.sfc.tile2.nc
-rw-r--r-- 1 George.Gayno g01 106197450 Apr 12 19:24 out.sfc.tile3.nc
-rw-r--r-- 1 George.Gayno g01 106197450 Apr 12 19:24 out.sfc.tile4.nc
-rw-r--r-- 1 George.Gayno g01 106197450 Apr 12 19:24 out.sfc.tile5.nc
-rw-r--r-- 1 George.Gayno g01 106197450 Apr 12 19:25 out.sfc.tile6.nc
/gpfs/hps/stmp/George.Gayno/chgres_save.sfc\ $ ncdump -h out.sfc.tile4.nc
netcdf out.sfc.tile4 {
dimensions:
lon = 768 ;
lat = 768 ;
lsoil = 4 ;
variables:
float lon(lon) ;
lon:cartesian_axis = "X" ;
float lat(lat) ;
lat:cartesian_axis = "Y" ;
float lsoil(lsoil) ;
lsoil:cartesian_axis = "Z" ;
float geolon(lat, lon) ;
float geolat(lat, lon) ;
float slmsk(lat, lon) ;
float tsea(lat, lon) ;
float sheleg(lat, lon) ;
float tg3(lat, lon) ;
float zorl(lat, lon) ;
float alvsf(lat, lon) ;
float alvwf(lat, lon) ;
float alnsf(lat, lon) ;
float alnwf(lat, lon) ;
float vfrac(lat, lon) ;
float canopy(lat, lon) ;
float f10m(lat, lon) ;
float t2m(lat, lon) ;
float q2m(lat, lon) ;
float vtype(lat, lon) ;
float stype(lat, lon) ;
float facsf(lat, lon) ;
float facwf(lat, lon) ;
float uustar(lat, lon) ;
float ffmm(lat, lon) ;
float ffhh(lat, lon) ;
float hice(lat, lon) ;
float fice(lat, lon) ;
float tisfc(lat, lon) ;
float tprcp(lat, lon) ;
float srflag(lat, lon) ;
float snwdph(lat, lon) ;
float shdmin(lat, lon) ;
float shdmax(lat, lon) ;
float slope(lat, lon) ;
float snoalb(lat, lon) ;
float stc(lsoil, lat, lon) ;
float smc(lsoil, lat, lon) ;
float slc(lsoil, lat, lon) ;
}
}}}
The changes at r91211 appear to work. Further tests to follow.
--
Ticket URL: <https://svnemc.ncep.noaa.gov/trac/fv3gfs/ticket/14#comment:6>
fv3gfs <https://svnemc.ncep.noaa.gov/trac/fv3gfs>
NGGPS FV3GFS Development
More information about the Ncep.list.emc.fv3gfs_tickets
mailing list