| 1 |
gforget |
1.1 |
;; bsub -Is -q geyser -W 4:00 -n 1 -P P35171148 ncl< JRA25.ncl |
| 2 |
|
|
|
| 3 |
|
|
;; On yellowstone: |
| 4 |
|
|
; OLD method from DSS |
| 5 |
|
|
;; copy from glade ... |
| 6 |
|
|
;; /glade/p/rda/data/ds625.0.fcst_phy2m/1979/197901.fcst_phy2m.subc60.tar |
| 7 |
|
|
;; note that JPL's tars are named 197901.fcst_phy2m.tar, and that |
| 8 |
|
|
;; their lat,lon values are slightly off from NCAR's version |
| 9 |
|
|
|
| 10 |
|
|
;; hsi get /DSS/DS625.0/fcst_phy2m/1979/'*subc60.tar' . |
| 11 |
|
|
|
| 12 |
|
|
;;tar -xvf 197901.fcst_phy2m.subc60.tar |
| 13 |
|
|
;;tar -xvf 197902.fcst_phy2m.subc60.tar |
| 14 |
|
|
;;tar -xvf 197903.fcst_phy2m.subc60.tar |
| 15 |
|
|
;;tar -xvf 197904.fcst_phy2m.subc60.tar |
| 16 |
|
|
;;tar -xvf 197905.fcst_phy2m.subc60.tar |
| 17 |
|
|
;;tar -xvf 197906.fcst_phy2m.subc60.tar |
| 18 |
|
|
;;tar -xvf 197907.fcst_phy2m.subc60.tar |
| 19 |
|
|
;;tar -xvf 197908.fcst_phy2m.subc60.tar |
| 20 |
|
|
;;tar -xvf 197909.fcst_phy2m.subc60.tar |
| 21 |
|
|
;;tar -xvf 197910.fcst_phy2m.subc60.tar |
| 22 |
|
|
;;tar -xvf 197911.fcst_phy2m.subc60.tar |
| 23 |
|
|
;;tar -xvf 197912.fcst_phy2m.subc60.tar |
| 24 |
|
|
;;mv 1979*/* . |
| 25 |
|
|
;;rm *tar |
| 26 |
|
|
;;rmdir 1979* |
| 27 |
|
|
|
| 28 |
|
|
;; ncl<JRA25.ncl |
| 29 |
|
|
;; bsub -Is -q geyser -W 4:00 -n 1 -P P35171148 ncl< JRA25.ncl |
| 30 |
|
|
|
| 31 |
|
|
;;; ncl_convert2nc ncdump -h |
| 32 |
|
|
;;;............................................ |
| 33 |
|
|
yr=2014 |
| 34 |
|
|
|
| 35 |
|
|
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" |
| 36 |
|
|
setfileoption("bin","WriteByteOrder","BigEndian") |
| 37 |
|
|
|
| 38 |
|
|
; Uroot1="mss:/DSS/DS627.0/ei.oper.an.sfc/" |
| 39 |
|
|
; Uroot2="/ei.oper.an.sfc.regn128sc." |
| 40 |
|
|
Oroot="fcst_phy2m." |
| 41 |
|
|
OOroot="JRA25g_" |
| 42 |
|
|
|
| 43 |
|
|
fils=systemfunc("ls " + Oroot + yr + "*") |
| 44 |
|
|
print (fils) |
| 45 |
|
|
nfils=dimsizes(fils) |
| 46 |
|
|
print (nfils) |
| 47 |
|
|
|
| 48 |
|
|
vars3d = (/ "LHTFL_GDS4_SFC_ave6h" , \ |
| 49 |
|
|
"SHTFL_GDS4_SFC_ave6h" , "U_FLX_GDS4_SFC_ave6h" , \ |
| 50 |
|
|
"V_FLX_GDS4_SFC_ave6h" , "USWRF_GDS4_SFC_ave6h" , \ |
| 51 |
|
|
"DSWRF_GDS4_SFC_ave6h" , "ULWRF_GDS4_SFC_ave6h" , \ |
| 52 |
|
|
"DLWRF_GDS4_SFC_ave6h" , "WTMP_GDS4_SFC_ave6h" , \ |
| 53 |
|
|
"ACPCP_GDS4_SFC_ave6h" , "NCPCP_GDS4_SFC_ave6h" , \ |
| 54 |
|
|
"T_CDC_GDS4_SFC_ave6h" , \ |
| 55 |
|
|
"U_GRD_GDS4_SFC" , "V_GRD_GDS4_SFC" , \ |
| 56 |
|
|
"TMP_GDS4_SFC" , "SPF_H_GDS4_SFC" , "PRES_GDS4_SFC", \ |
| 57 |
|
|
"T_CDC_GDS4_SFC_ave6h" /) |
| 58 |
|
|
|
| 59 |
|
|
names3d = (/ "xlhf_" , "xshf_" , "xustr_" , \ |
| 60 |
|
|
"xvstr_" , "xusw_" , "xdsw_" , \ |
| 61 |
|
|
"xulw_" , "xdlw_" , "xtmp_degC_" , \ |
| 62 |
|
|
"xacpcp_" , "xncpcp_" , "xcldcv_" , "u10m_" , "v10m_" , \ |
| 63 |
|
|
"tmp2m_degC_" , "spfh2m_" , \ |
| 64 |
|
|
"pres_", \ |
| 65 |
|
|
"xcldcv_" /) |
| 66 |
|
|
offset = (/0, 0, 0, 0, 0, 0, 0, 0, -273.15, 0, 0, 0, 0, 0, -273.15, 0, 0, 0 /) |
| 67 |
|
|
scalef = (/1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, .01, 1, 1, 1, 1, 1, .01 /) |
| 68 |
|
|
nvars=dimsizes (vars3d) |
| 69 |
|
|
do nv=0,nvars-1 |
| 70 |
|
|
filbin= OOroot + names3d(nv) + yr ; name of binary output file |
| 71 |
|
|
system("/bin/rm -f "+filbin) ; remove if it already exist |
| 72 |
|
|
end do |
| 73 |
|
|
|
| 74 |
|
|
|
| 75 |
|
|
|
| 76 |
|
|
do nf=0,nfils-1 |
| 77 |
|
|
f = addfile (fils(nf)+".grb", "r") ; open |
| 78 |
|
|
print(""+fils(nf)) |
| 79 |
|
|
|
| 80 |
|
|
do nv=0,nvars-1 ; loop over selected variables |
| 81 |
|
|
;print("variable: "+vars3d(nv)) |
| 82 |
|
|
filbin= OOroot + names3d(nv) + yr ; name of binary output file |
| 83 |
|
|
; print(filbin) |
| 84 |
|
|
x = f->$vars3d(nv)$ |
| 85 |
|
|
;printVarSummary( x ) |
| 86 |
|
|
x = x*scalef(nv) + offset(nv) |
| 87 |
|
|
X = x(::-1,:) ; leftmost dimension is slowest varying |
| 88 |
|
|
dimx = dimsizes( X ) |
| 89 |
|
|
fbindirwrite(filbin, X ) ; write to binary file |
| 90 |
|
|
|
| 91 |
|
|
end do |
| 92 |
|
|
|
| 93 |
|
|
delete( x ) ; may change size on the next |
| 94 |
|
|
delete( X ) ; iterations due to different |
| 95 |
|
|
end do |
| 96 |
|
|
|
| 97 |
|
|
;----------------------------------------------------------------- |
| 98 |
|
|
; Remove the current year of files ... then return and read next 12 files |
| 99 |
|
|
;----------------------------------------------------------------- |
| 100 |
|
|
; SYS_RM = "/bin/rm Uroot+"*" |
| 101 |
|
|
; print ( SYS_RM ) |
| 102 |
|
|
;; system( SYS_RM ) ; not until debugged |
Parent Directory
| 
Revision Log
| 
Revision Graph