Changeset 175 for trunk/src/tropflux.pro

Timestamp:

03/22/12 14:40:17 (12 years ago)

Author:

pinsard

Message:

an other bunch of new functions

File:

• trunk/src/tropflux.pro (moved) (moved from trunk/src/TropFlux_19890101_20091231.pro) (34 diffs, 1 prop)

trunk/src/tropflux.pro

@@ -1 +1 @@
 ;+
 ;
-; .. _TropFlux_19890101_20091231.pro:
-;
-; ==============================
-; TropFlux_19890101_20091231.pro
-; ==============================
+; ============
+; tropflux.pro
+; ============
+;
+; .. function:: tropflux(yyyymmddb,yyyymmdde)
+;
+; DESCRITION
+; ==========
 ;
 ; This program computes net heat flux components on the 1° oaflux grid.
@@ -20 +23 @@
 ; :func:`oaflux_mask_30n30s`.
 ;
-; :file:`${PROJECT_OD}/TropFlux_sst_19890101_20091231.nc`
+; :file:`${PROJECT_OD}/TropFlux_sst_{yyyymmdd}_{yyyymmdd}.nc`
 ; containing
 ; sst corrected on OAFLUX grid
 ; has been produced by
-; :ref:`sst_correction_ncdf.pro`.
-;
-; :file:`${PROJECT_OD}/TropFlux_ws_19890101_20091231.nc`
+; :func:`sst_correction_ncdf`.
+;
+; :file:`${PROJECT_OD}/TropFlux_ws_{yyyymmdd}_{yyyymmdd}.nc`
 ; containing
 ; ws corrected on OAFLUX grid
 ; has been produced by
-; :ref:`ws_correction_ncdf.pro`.
-;
-; :file:`${PROJECT_OD}/TropFlux_gustiness_19890101_20091231.nc`
+; :func:`ws_correction_ncdf`.
+;
+; :file:`${PROJECT_OD}/TropFlux_gustiness_{yyyymmdd}_{yyyymmdd}.nc`
 ; containing
 ; ++
@@ -38 +41 @@
 ; :func:`cronin_gustiness_ncdf`.
 ;
-; :file:`${PROJECT_OD}/TropFlux_swr_19890101_20091231_BLND.nc`
+; :file:`${PROJECT_OD}/TropFlux_swr_{yyyymmdd}_{yyyymmdd}_BLND.nc`
 ; containing
 ; ws corrected on OAFLUX grid
@@ -44 +47 @@
 ; :func:`tropflux_swr_blnd`.
 ;
-; :file:`${PROJECT_OD}/TropFlux_lwr_19890101_20091231.nc`
+; :file:`${PROJECT_OD}/TropFlux_lwr_{yyyymmdd}_{yyyymmdd}.nc`
 ; containing
 ; lwr corrected on OAFLUX grid
 ; has been produced by
-; :ref:`lwr_correction_ncdf.pro`.
+; :func:`lwr_correction_ncdf`.
 ;
 ; containing
 ; t2m corrected on OAFLUX grid
 ; has been produced by
-; :ref:`t2m_correction_ncdf.pro`.
-;
-; :file:`${PROJECT_OD}/TropFlux_q2m_19890101_20091231.nc`
+; :func:`t2m_correction_ncdf`.
+;
+; :file:`${PROJECT_OD}/TropFlux_q2m_{yyyymmdd}_{yyyymmdd}.nc`
 ; containing
 ; q2m corrected on OAFLUX grid
@@ -62 +65 @@
 ;
 ; net heat flux components are written
-; in :file:`${PROJECT_OD}/TropFlux_19890101_20091231_coarev3.nc`
+; in :file:`${PROJECT_OD}/TropFlux_{yyyymmdd}_{yyyymmdd}_coarev3.nc`
 ; if this file not already exists.
 ;
 ; This output file
-; :file:`${PROJECT_OD}/TropFlux_19890101_20091231_coarev3.nc`
+; :file:`${PROJECT_OD}/TropFlux_{yyyymmdd}_{yyyymmdd}_coarev3.nc`
 ; will be used by
 ; :ref:`TropFlux_NRT_ncdf.pro`.
@@ -72 +75 @@
 ;     .. graphviz::
 ;
-;        digraph tropflux_19890101_20091231 {
+;        digraph tropflux {
 ;
 ;           mask [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/mask_oaflux_30N30S.nc"];
-;           file_sst [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_sst_19890101_20091231.nc"];
-;           file_ws [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_ws_19890101_20091231.nc"];
-;           file_wg [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_gustiness_19890101_20091231.nc"];
-;           file_swr [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_swr_19890101_20091231_BLND.nc"];
-;           file_lwr [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_lwr_19890101_20091231.nc"];
-;           file_t2m [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_t2m_19890101_20091231.nc"];
-;           file_q2m [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_q2m_19890101_20091231.nc"];
-;
-;           file_out[shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_19890101_20091231_coarev3.nc"];
-;
-;           tropflux_19890101_20091231 [shape=box,
+;           file_sst [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_sst_{yyyymmdd}_{yyyymmdd}.nc"];
+;           file_ws [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_ws_{yyyymmdd}_{yyyymmdd}.nc"];
+;           file_wg [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_gustiness_{yyyymmdd}_{yyyymmdd}.nc"];
+;           file_swr [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_swr_{yyyymmdd}_{yyyymmdd}_BLND.nc"];
+;           file_lwr [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_lwr_{yyyymmdd}_{yyyymmdd}.nc"];
+;           file_t2m [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_t2m_{yyyymmdd}_{yyyymmdd}.nc"];
+;           file_q2m [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_q2m_{yyyymmdd}_{yyyymmdd}.nc"];
+;
+;           file_out[shape=ellipse,fontname=Courier,label="${PROJECT_OD}/TropFlux_{yyyymmdd}_{yyyymmdd}_coarev3.nc"];
+;
+;           tropflux [shape=box,
 ;           fontname=Courier,
 ;           color=blue,
-;           URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/TropFlux_19890101_20091231.pro",
-;           label="${PROJECT}/src/TropFlux_19890101_20091231.pro"];
-;
-;           {mask file_sst file_ws file_wg file_swr file_lwr file_t2m file_q2m} -> {tropflux_19890101_20091231} -> {file_out}
+;           URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/tropflux.pro",
+;           label="${PROJECT}/src/tropflux.pro"];
+;
+;           {mask file_sst file_ws file_wg file_swr file_lwr file_t2m file_q2m} -> {tropflux} -> {file_out}
 ;
 ;        }
@@ -99 +102 @@
 ;
 ; :ref:`project_profile.sh`
+;
+; Used by :ref:`tropflux.sh`
 ;
 ; :func:`report <saxo:report>`
@@ -124 +129 @@
 ; ::
 ;
-;  IDL> .compile TropFlux_19890101_20091231
-;  IDL> tropflux_19890101_20091231
-;
+;  IDL> yyyymmddb = 20000101L
+;  IDL> yyyymmdde = 20001231L
+;  IDL> result = tropflux(yyyymmddb, yyyymmdde)
+;  IDL> print, result
+
 ; TODO
 ; ====
 ;
 ; describe usage of tau
+;
+; compare time between all inputs after reading
 ;
 ; make it work ::
@@ -180 +189 @@
 ;   wg=read_ncdf('wg',da1-1,da2,file=fullfilename_wg,/nostr)
 ;
-; avoid mix lower/uppercase in pro name to avoid compile
-;
 ; coding rules
 ;
@@ -196 +203 @@
 ;
 ; $URL$
+;
+; - fplod 20120322
+;
+;   * pro -> func
+;   * rename with lower case TropFlux_19890101_20091231 become tropflux
+;   * taking project_overwite into account
+;   * try to add compile_opt
+;   * hard coded da1 and da2 replaced by yyyymmddb and yyyymmdde parameters
+;   * get rid of timegen
 ;
 ; - fplod 20110830T135832Z cratos (Linux)
@@ -241 +257 @@
 ;-
 ;
-pro TropFlux_19890101_20091231
+function tropflux $
+         , yyyymmddb $
+         , yyyymmdde
+;
+compile_opt idl2, strictarrsubs, logical_predicate
 ;
 @cm_4cal
@@ -249 +269 @@
 @cm_project
 ;
+; Return to caller if errors
+ON_ERROR, 2
+;
+result = -1
+;
+usage = 'result = tropflux_swr_nrt(yyyymmddb, yyyymmdde)'
+nparam = N_PARAMS()
+IF (nparam NE 2) THEN BEGIN
+    ras = report(['Incorrect number of arguments.' $
+          + '!C' $
+          + 'Usage : ' + usage])
+    return, result
+ENDIF
+;
 ; check for input directory
 ;
@@ -256 +290 @@
      msg = 'eee : ${PROJECT_ID} is not defined'
      ras = report(msg)
-     STOP
+     return, result
            END
  ELSE: BEGIN
@@ -270 +304 @@
    msg = 'eee : the directory' + iodirin  + ' is not accessible.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
@@ -283 +317 @@
    msg = 'eee : the file ' + fullfilename_msk + ' was not found.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
@@ -291 +325 @@
          msg = 'eee : ${PROJECT_OD} is not defined'
          ras = report(msg)
-       STOP
+       return, result
        END
   ELSE: BEGIN
@@ -306 +340 @@
    msg = 'eee : the directory' + iodirout  + ' is not accessible.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
@@ -313 +347 @@
     msg = 'eee : the directory' + iodirout  + ' was not found.'
     ras = report(msg)
-    STOP
+    return, result
 ENDIF
 ;
 ; build sst data filename
-filename_sst='TropFlux_sst_19890101_20091231.nc'
+filename_sst='TropFlux_sst_' +  string(yyyymmddb,format='(I8.8)') + '_' + string(yyyymmdde,format='(I8.8)') + '.nc'
 ;
 ; check if this file exists
@@ -326 +360 @@
    msg = 'eee : the file ' + fullfilename_sst + ' was not found.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
 ; build ws data filename
-filename_ws='TropFlux_ws_19890101_20091231.nc'
+filename_ws='TropFlux_ws_' +  string(yyyymmddb,format='(I8.8)') + '_' + string(yyyymmdde,format='(I8.8)') + '.nc'
 ;
 ; check if this file exists
@@ -339 +373 @@
    msg = 'eee : the file ' + fullfilename_ws + ' was not found.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
 ; build swr data filename
-filename_swr='TropFlux_swr_19890101_20091231_BLND.nc'
+filename_swr='TropFlux_swr_' +  string(yyyymmddb,format='(I8.8)') + '_' + string(yyyymmdde,format='(I8.8)') + '_BLND.nc'
 ;
 ; check if this file exists
@@ -352 +386 @@
    msg = 'eee : the file ' + fullfilename_swr + ' was not found.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
 ; build lwr data filename
-filename_lwr='TropFlux_lwr_19890101_20091231.nc'
+filename_lwr='TropFlux_lwr_' +  string(yyyymmddb,format='(I8.8)') + '_' + string(yyyymmdde,format='(I8.8)') + '.nc'
 ;
 ; check if this file exists
@@ -365 +399 @@
    msg = 'eee : the file ' + fullfilename_lwr + ' was not found.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
 ; build t2m data filename
-filename_t2m='TropFlux_t2m_19890101_20091231.nc'
+filename_t2m='TropFlux_t2m_' +  string(yyyymmddb,format='(I8.8)') + '_' + string(yyyymmdde,format='(I8.8)') + '.nc'
 ;
 ; check if this file exists
@@ -378 +412 @@
    msg = 'eee : the file ' + fullfilename_t2m + ' was not found.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
 ; build q2m data filename
-filename_q2m='TropFlux_q2m_19890101_20091231.nc'
+filename_q2m='TropFlux_q2m_' +  string(yyyymmddb,format='(I8.8)') + '_' + string(yyyymmdde,format='(I8.8)') + '.nc'
 ;
 ; check if this file exists
@@ -391 +425 @@
    msg = 'eee : the file ' + fullfilename_q2m + ' was not found.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
 ; build wg data filename
-filename_wg='TropFlux_gustiness_19890101_20091231.nc'
+filename_wg='TropFlux_gustiness_' +  string(yyyymmddb,format='(I8.8)') + '_' + string(yyyymmdde,format='(I8.8)') + '.nc'
 ;
 ; check if this file exists
@@ -404 +438 @@
    msg = 'eee : the file ' + fullfilename_wg + ' was not found.'
    ras = report(msg)
-   STOP
+   return, result
 ENDIF
 ;
 ; build output filename
-filename_out = 'TropFlux_19890101_20091231_coarev3.nc'
+filename_out = 'TropFlux_' +  string(yyyymmddb,format='(I8.8)') + '_' + string(yyyymmdde,format='(I8.8)') + '_coarev3.nc'
 fullfilename_out = iodirout + filename_out
 ; in order to avoid unexpected overwritten
@@ -414 +448 @@
    msg = 'eee : the file ' + fullfilename_out  + ' already exists.'
    ras = report(msg)
-   STOP
-ENDIF
-;
-da1=19880101
-da2=20101231
+   return, result
+ENDIF
 ;
 initncdf, fullfilename_msk
@@ -425 +456 @@
 initncdf, fullfilename_sst
 ;
-ws=read_ncdf('ws',da1,da2,file=fullfilename_ws,/nostr)
-wg=read_ncdf('wg',da1-1,da2,file=fullfilename_wg,/nostr)
+ws=read_ncdf('ws',yyyymmddb-.5d,yyyymmdde,file=fullfilename_ws,/nostr)
+timein=24.d*(time-julday(1,1,1957,0,0,0))
+jpt=n_elements(timein)
+da=jul2date(time[0])
+cda0=string(da,format='(i8.8)')
+da=jul2date(time[jpt-1])
+cda1=string(da,format='(i8.8)')
+print, 'ws in tropflux first date ', cda0
+print, 'ws in ws_correction_ncdf last date ' , cda1
+;
+wg=read_ncdf('wg',yyyymmddb-.5d-1,yyyymmdde,file=fullfilename_wg,/nostr)
+timein=24.d*(time-julday(1,1,1957,0,0,0))
+jpt=n_elements(timein)
+da=jul2date(time[0])
+cda0=string(da,format='(i8.8)')
+da=jul2date(time[jpt-1])
+cda1=string(da,format='(i8.8)')
+print, 'wg in tropflux first date ', cda0
+print, 'wg in ws_correction_ncdf last date ' , cda1
 tt=time
 jpt=n_elements(time)
-sst=read_ncdf('sst',da1,da2,file=fullfilename_sst,/nostr)
+sst=read_ncdf('sst',yyyymmddb-.5d,yyyymmdde,file=fullfilename_sst,/nostr)
+timein=24.d*(time-julday(1,1,1957,0,0,0))
+jpt=n_elements(timein)
+da=jul2date(time[0])
+cda0=string(da,format='(i8.8)')
+da=jul2date(time[jpt-1])
+cda1=string(da,format='(i8.8)')
+print, 'sst in tropflux first date ', cda0
+print, 'sst in ws_correction_ncdf last date ' , cda1
 sst=reform(sst-273.15)
-swd=read_ncdf('swr',da1,da2,file=fullfilename_swr,/nostr)
-lw=read_ncdf('lwr',da1,da2,file=fullfilename_lwr,/nostr)
+swd=read_ncdf('swr',yyyymmddb-.5d,yyyymmdde,file=fullfilename_swr,/nostr)
+timein=24.d*(time-julday(1,1,1957,0,0,0))
+jpt=n_elements(timein)
+da=jul2date(time[0])
+cda0=string(da,format='(i8.8)')
+da=jul2date(time[jpt-1])
+cda1=string(da,format='(i8.8)')
+print, 'swr in tropflux first date ', cda0
+print, 'swr in ws_correction_ncdf last date ' , cda1
+lw=read_ncdf('lwr',yyyymmddb-.5d,yyyymmdde,file=fullfilename_lwr,/nostr)
+timein=24.d*(time-julday(1,1,1957,0,0,0))
+jpt=n_elements(timein)
+da=jul2date(time[0])
+cda0=string(da,format='(i8.8)')
+da=jul2date(time[jpt-1])
+cda1=string(da,format='(i8.8)')
+print, 'lwr in tropflux first date ', cda0
+print, 'lwr in ws_correction_ncdf last date ' , cda1
 swd=swd/0.94  ; converting from net swr to downward swr
 ;
-t2m=read_ncdf('t2m',da1,da2,file=fullfilename_t2m,/nostr)-273.15   ; in C
-q2m=read_ncdf('q2m',da1,da2,file=fullfilename_q2m,/nostr)    ; in g/kg
+t2m=read_ncdf('t2m',yyyymmddb-.5d,yyyymmdde,file=fullfilename_t2m,/nostr)-273.15   ; in C
+timein=24.d*(time-julday(1,1,1957,0,0,0))
+jpt=n_elements(timein)
+da=jul2date(time[0])
+cda0=string(da,format='(i8.8)')
+da=jul2date(time[jpt-1])
+cda1=string(da,format='(i8.8)')
+print, 't2m in tropflux first date ', cda0
+print, 't2m in ws_correction_ncdf last date ' , cda1
+q2m=read_ncdf('q2m',yyyymmddb-.5d,yyyymmdde,file=fullfilename_q2m,/nostr)    ; in g/kg
+timein=24.d*(time-julday(1,1,1957,0,0,0))
+jpt=n_elements(timein)
+da=jul2date(time[0])
+cda0=string(da,format='(i8.8)')
+da=jul2date(time[jpt-1])
+cda1=string(da,format='(i8.8)')
+print, 'q2m in tropflux first date ', cda0
+print, 'q2m in ws_correction_ncdf last date ' , cda1
 ;
 w=sqrt(ws*ws+wg*wg)  ; wind corrected for gustiness
@@ -471 +559 @@
 ;
 for jt=0,jpt-1 do begin
-  jday=time(jt)-julday(1,1,yea(jt))
+  jday=time[jt]-julday(1,1,yea[jt])
   print, 'Computing Fluxes ',jt,' / ',jpt-1
 ;
@@ -478 +566 @@
 ;
   ; wind speed (m/s)
-  wn=w(*,*,jt)
-  wn=wn(ocean)
+  wn=w[*,*,jt]
+  wn=wn[ocean]
   ; Bulk sst (°C)
-  ts=sst(*,*,jt)
-  ts=ts(ocean)
+  ts=sst[*,*,jt]
+  ts=ts[ocean]
   ; 2m Air T (°C)
-  t=t2m(*,*,jt)
-  t=t(ocean)
+  t=t2m[*,*,jt]
+  t=t[ocean]
   ; Sea surface sat. spec. humidity (g/kg)
   qs=qsee(ts,P)
   ; 2m AIr specific humidity  (g/kg)
-  q=q2m(*,*,jt)
-  q=q(ocean)
+  q=q2m[*,*,jt]
+  q=q[ocean]
   ; Downward solar flux (W/m2)
-  Rs=swd(*,*,jt)
-  Rs=Rs(ocean)
-  ylat=gphit(ocean)
+  Rs=swd[*,*,jt]
+  Rs=Rs[ocean]
+  ylat=gphit[ocean]
   ;  cld=calc_cloud(jday,Rs,ylat)
   cld=calc_cloud_vlat(jday,Rs,ylat)
   ; Downward IR flux (W/m2)
   ;  Rl=lwdown_clark(ts,q,cld,t,P)
-  Rl=lw(*,*,jt)
-  Rl=Rl(ocean)
+  Rl=lw[*,*,jt]
+  Rl=Rl[ocean]
   rain=0.
   lw_clrk=-lwnet_clark(ts,q,cld,t,P)
-;  junk(*,*,jt)=lw_clrk
+;  junk[*,*,jt]=lw_clrk
 ;
 ;stop
@@ -509 +597 @@
 ; A few punctual missing values (coare does not converge): filled by spatial extrapolation
   tab=fltarr(jpi,jpj)+!values.f_nan
-  x=reform(y(*,0))
-  tab(ocean)=x
+  x=reform(y[*,0])
+  tab[ocean]=x
   m=finite(tab)
   help, tab, m, swr
   tab=extrapolate(tab,m)
-  swr(*,*,jt)=tab*msk+valmask*(1-msk)
+  swr[*,*,jt]=tab*msk+valmask*(1-msk)
 ;
   tab=fltarr(jpi,jpj)+!values.f_nan
-  x=reform(y(*,1))
-  tab(ocean)=x
+  x=reform(y[*,1])
+  tab[ocean]=x
   m=finite(tab)
   help, tab, m, lwr
@@ -524 +612 @@
   print, 'www : extrapolation for lwr is not done'
   tab=tab
-  lwr(*,*,jt)=tab*msk+valmask*(1-msk)
+  lwr[*,*,jt]=tab*msk+valmask*(1-msk)
 ;
   tab=fltarr(jpi,jpj)+!values.f_nan
-  x=reform(y(*,2))
-  tab(ocean)=x
+  x=reform(y[*,2])
+  tab[ocean]=x
   m=finite(tab)
   help, tab, m, lat
@@ -534 +622 @@
   print, 'www : extrapolation for lat is not done'
   tab=tab
-  lat(*,*,jt)=tab*msk+valmask*(1-msk)
+  lat[*,*,jt]=tab*msk+valmask*(1-msk)
   tab=fltarr(jpi,jpj)+!values.f_nan
-  x=reform(y(*,3))
-  tab(ocean)=x
+  x=reform(y[*,3])
+  tab[ocean]=x
   m=finite(tab)
   help, tab, m, sen
@@ -543 +631 @@
   print, 'www : extrapolation for sen is not done'
   tab=tab
-  sen(*,*,jt)=tab*msk+valmask*(1-msk)
-  ;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(lw_clrk) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & lwnet_clrk(*,*,jt)=tab*msk+valmask*(1-msk)
-  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,5)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & tau(*,*,jt)=tab*msk+valmask*(1-msk)
-  ;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,6)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & Ch(*,*,jt)=tab*msk+valmask*(1-msk)
-  ;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,7)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & Ce(*,*,jt)=tab*msk+valmask*(1-msk)
-  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,8)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & wg(*,*,jt)=tab*msk+valmask*(1-msk)
+  sen[*,*,jt]=tab*msk+valmask*(1-msk)
+  ;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(lw_clrk) & tab[ocean]=x & m=finite(tab) & tab=extrapolate(tab,m) & lwnet_clrk[*,*,jt]=tab*msk+valmask*(1-msk)
+  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y[*,5]) & tab[ocean]=x & m=finite(tab) & tab=extrapolate(tab,m) & tau[*,*,jt]=tab*msk+valmask*(1-msk)
+  ;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y[*,6]) & tab[ocean]=x & m=finite(tab) & tab=extrapolate(tab,m) & Ch[*,*,jt]=tab*msk+valmask*(1-msk)
+  ;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y[*,7]) & tab[ocean]=x & m=finite(tab) & tab=extrapolate(tab,m) & Ce[*,*,jt]=tab*msk+valmask*(1-msk)
+  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y[*,8]) & tab[ocean]=x & m=finite(tab) & tab=extrapolate(tab,m) & wg[*,*,jt]=tab*msk+valmask*(1-msk)
 endfor
-time=timegen(7670, start=julday(1,1,1989,0), units='days')
 jpt=n_elements(time)
 ;
-cda0=string(jul2date(time(0)),format='(i8.8)')
-cda1=string(jul2date(time(jpt-1)),format='(i8.8)')
+cda0=string(jul2date(time[0]),format='(i8.8)')
+cda1=string(jul2date(time[jpt-1]),format='(i8.8)')
 time=time-julday(1,1,1950,00,00,00)
-xlon=reform(glamt(*,0) )
-ylat=reform(gphit(0,*))
+xlon=reform(glamt[*,0])
+ylat=reform(gphit[0,*])
 ;
 ncfile='!' + fullfilename_out
@@ -590 +677 @@
 @ncdf_quickwrite
 ;
+result = 0
+return, result
+;
 end