;;  bsub -Is -q geyser -W 4:00 -n 1 -P  P35171148 ncl< JRA25.ncl

;;  On yellowstone:
;   OLD method from DSS
;;  copy from glade  ...
;;  /glade/p/rda/data/ds625.0.fcst_phy2m/1979/197901.fcst_phy2m.subc60.tar
;;     note that JPL's tars are named 197901.fcst_phy2m.tar, and that
;;        their lat,lon values are slightly off from NCAR's version

;; hsi get /DSS/DS625.0/fcst_phy2m/1979/'*subc60.tar' .

;;tar -xvf 197901.fcst_phy2m.subc60.tar
;;tar -xvf 197902.fcst_phy2m.subc60.tar
;;tar -xvf 197903.fcst_phy2m.subc60.tar
;;tar -xvf 197904.fcst_phy2m.subc60.tar
;;tar -xvf 197905.fcst_phy2m.subc60.tar
;;tar -xvf 197906.fcst_phy2m.subc60.tar
;;tar -xvf 197907.fcst_phy2m.subc60.tar
;;tar -xvf 197908.fcst_phy2m.subc60.tar
;;tar -xvf 197909.fcst_phy2m.subc60.tar
;;tar -xvf 197910.fcst_phy2m.subc60.tar
;;tar -xvf 197911.fcst_phy2m.subc60.tar
;;tar -xvf 197912.fcst_phy2m.subc60.tar
;;mv 1979*/* .
;;rm *tar
;;rmdir 1979*

;; ncl<JRA25.ncl
;;  bsub -Is -q geyser -W 4:00 -n 1 -P  P35171148 ncl< JRA25.ncl 

;;;                 ncl_convert2nc    ncdump -h 
;;;............................................
yr=2014

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
setfileoption("bin","WriteByteOrder","BigEndian") 

;  Uroot1="mss:/DSS/DS627.0/ei.oper.an.sfc/"
;  Uroot2="/ei.oper.an.sfc.regn128sc."
Oroot="fcst_phy2m."
OOroot="JRA25g_"

fils=systemfunc("ls " + Oroot + yr + "*")
print (fils)
nfils=dimsizes(fils)
print (nfils)

vars3d = (/ "LHTFL_GDS4_SFC_ave6h" ,  \
   "SHTFL_GDS4_SFC_ave6h" , "U_FLX_GDS4_SFC_ave6h" , \
   "V_FLX_GDS4_SFC_ave6h" , "USWRF_GDS4_SFC_ave6h" , \
   "DSWRF_GDS4_SFC_ave6h" , "ULWRF_GDS4_SFC_ave6h" , \
   "DLWRF_GDS4_SFC_ave6h" , "WTMP_GDS4_SFC_ave6h" , \
   "ACPCP_GDS4_SFC_ave6h" , "NCPCP_GDS4_SFC_ave6h" ,  \
   "T_CDC_GDS4_SFC_ave6h" , \
   "U_GRD_GDS4_SFC" , "V_GRD_GDS4_SFC" , \
   "TMP_GDS4_SFC" , "SPF_H_GDS4_SFC" , "PRES_GDS4_SFC", \
   "T_CDC_GDS4_SFC_ave6h"  /)

names3d = (/ "xlhf_" , "xshf_" , "xustr_" , \
   "xvstr_" , "xusw_" , "xdsw_" , \
   "xulw_" , "xdlw_" , "xtmp_degC_" , \
   "xacpcp_" , "xncpcp_"  , "xcldcv_" , "u10m_" , "v10m_" , \
   "tmp2m_degC_" , "spfh2m_" , \
   "pres_", \
   "xcldcv_" /)
offset = (/0, 0, 0, 0, 0, 0, 0, 0, -273.15, 0, 0, 0, 0, 0, -273.15, 0, 0, 0 /)
scalef = (/1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, .01, 1, 1, 1, 1, 1, .01 /)
nvars=dimsizes (vars3d)
do nv=0,nvars-1
	filbin= OOroot + names3d(nv) + yr  ; name of binary output file
	system("/bin/rm -f "+filbin)            ; remove if it already exist
end do



     do nf=0,nfils-1
        f = addfile (fils(nf)+".grb", "r")  ; open
        print(""+fils(nf))
     
        do nv=0,nvars-1                   ; loop over selected variables
          ;print("variable: "+vars3d(nv))
	   filbin= OOroot + names3d(nv) + yr  ; name of binary output file
;          print(filbin)
           x = f->$vars3d(nv)$
          ;printVarSummary( x ) 
	   x =  x*scalef(nv) + offset(nv)
	   X = x(::-1,:)  ; leftmost dimension is slowest varying
           dimx = dimsizes( X )
           fbindirwrite(filbin, X )       ; write to binary file    

        end do
        
        delete( x )                   ; may change size on the next
        delete( X )                   ; iterations due to different
     end do

;-----------------------------------------------------------------
; Remove the current year of files  ... then return and read next 12 files
;-----------------------------------------------------------------
;    SYS_RM = "/bin/rm Uroot+"*"
;    print ( SYS_RM )
;; system( SYS_RM )    ; not until debugged