Changeset 204 for trunk/src/interp_erai_ws.pro

Timestamp:

03/26/14 11:58:56 (10 years ago)

Author:

pinsard

Message:

fix thanks to coding rules; typo

File:

• trunk/src/interp_erai_ws.pro (modified) (10 diffs)

trunk/src/interp_erai_ws.pro

@@ -198 +198 @@
 ; - fplod 20120320
 ;
-;   * taking project_overwite into account
+;   * taking project_overwrite into account
 ;   * try to add compile_opt seems to be incompatible with ncdf_quickwrite
 ;   * pro -> function
@@ -252 +252 @@
 nparam = N_PARAMS()
 IF (nparam NE 3) THEN BEGIN
-ras = report(['Incorrect number of arguments.' $
-      + '!C' $
-      + 'Usage : ' + usage])
-      return, result
+    ras = report(['Incorrect number of arguments.' $
+    + '!C' $
+    + 'Usage : ' + usage])
+    return, result
 ENDIF
 ;
@@ -268 +268 @@
 ; test if ${PROJECT_ID} defined
 CASE project_id_env OF
-    ''  :  BEGIN
-     msg = 'eee : ${PROJECT_ID} is not defined'
-     ras = report(msg)
-     return, result
-           END
- ELSE: BEGIN
-     msg = 'iii : ${PROJECT_ID} is ' + project_id_env
-     ras = report(msg)
-       END
+    '' : BEGIN
+        msg = 'eee : ${PROJECT_ID} is not defined'
+        ras = report(msg)
+        return, result
+    END
+    ELSE : BEGIN
+        msg = 'iii : ${PROJECT_ID} is ' + project_id_env
+        ras = report(msg)
+    END
 ENDCASE
 ;
@@ -283 +283 @@
 ; existence and protection of ${PROJECT_ID}
 IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN
-   msg = 'eee : the directory' + iodirin  + ' is not accessible.'
-   ras = report(msg)
-   return, result
+    msg = 'eee : the directory' + iodirin  + ' is not accessible.'
+    ras = report(msg)
+    return, result
 ENDIF
 ;
@@ -296 +296 @@
 fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST)
 IF fullfilename_msk[0] EQ '' THEN BEGIN
-   msg = 'eee : the file ' + fullfilename_msk + ' was not found.'
-   ras = report(msg)
-   return, result
+    msg = 'eee : the file ' + fullfilename_msk + ' was not found.'
+    ras = report(msg)
+    return, result
 ENDIF
 ;
@@ -322 +322 @@
 fullfilename = isafile(iodirin + filename, NEW=0, /MUST_EXIST)
 IF fullfilename[0] EQ '' THEN BEGIN
-   msg = 'eee : the file ' + fullfilename + ' was not found.'
-   ras = report(msg)
-   return, result
+    msg = 'eee : the file ' + fullfilename + ' was not found.'
+    ras = report(msg)
+    return, result
 ENDIF
 ;
 ; test if ${PROJECT_OD} defined
 CASE project_od_env OF
-  '' : BEGIN
-         msg = 'eee : ${PROJECT_OD} is not defined'
-         ras = report(msg)
-       return, result
-       END
-  ELSE: BEGIN
-          msg = 'iii : ${PROJECT_OD} is ' + project_od_env
-          ras = report(msg)
-        END
- ENDCASE
+    '' : BEGIN
+        msg = 'eee : ${PROJECT_OD} is not defined'
+        ras = report(msg)
+        return, result
+    END
+    ELSE : BEGIN
+        msg = 'iii : ${PROJECT_OD} is ' + project_od_env
+        ras = report(msg)
+    END
+ENDCASE
 ;
 ; check if output data will be possible
@@ -355 +355 @@
 ; in order to avoid unexpected overwritten
 IF ((FILE_TEST(fullfilename_out) EQ 1)  AND (project_overwrite EQ 0)) THEN BEGIN
-   msg = 'eee : the file ' + fullfilename_out  + ' already exists.'
-   ras = report(msg)
-   return, result
+    msg = 'eee : the file ' + fullfilename_out  + ' already exists.'
+    ras = report(msg)
+    return, result
 ENDIF
 ;
@@ -389 +389 @@
 fullfilename = isafile(iodirin + filename, NEW=0, /MUST_EXIST)
 IF fullfilename[0] EQ '' THEN BEGIN
-   msg = 'eee : the file ' + fullfilename + ' was not found.'
-   ras = report(msg)
-   return, result
+    msg = 'eee : the file ' + fullfilename + ' was not found.'
+    ras = report(msg)
+    return, result
 ENDIF
 initncdf, fullfilename
@@ -430 +430 @@
 ;
 for jt=0,jptin-1 do begin
-  ; ++print, 'Interpolation jt=',jt,' / ',jptin-1
-  tab=reform(u10in[*,*,jt])
-  u10out[*,*,jt]=call_interp2d(tab,lonin,latin,mskin $
-      , lonout,latout,method='bilinear' $
-      , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout)
-  u10out[*,*,jt]=u10out[*,*,jt]*mskout+(1.-mskout)*1.e20
-;
-  tab=reform(v10in[*,*,jt])
-  v10out[*,*,jt]=call_interp2d(tab,lonin,latin,mskin $
-      , lonout,latout,method='bilinear' $
-      , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout)
-  v10out[*,*,jt]=v10out[*,*,jt]*mskout+(1.-mskout)*1.e20
-;
+    ; ++print, 'Interpolation jt=',jt,' / ',jptin-1
+    tab=reform(u10in[*,*,jt])
+    u10out[*,*,jt]=call_interp2d(tab,lonin,latin,mskin $
+    , lonout,latout,method='bilinear' $
+    , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout)
+    u10out[*,*,jt]=u10out[*,*,jt]*mskout+(1.-mskout)*1.e20
+    ;
+    tab=reform(v10in[*,*,jt])
+    v10out[*,*,jt]=call_interp2d(tab,lonin,latin,mskin $
+    , lonout,latout,method='bilinear' $
+    , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout)
+    v10out[*,*,jt]=v10out[*,*,jt]*mskout+(1.-mskout)*1.e20
+    ;
 endfor
 ;
@@ -466 +466 @@
 ;
 ncfields = 'u10[longitude,latitude,*time]=u10out:u10_attr; ' $
-          +'v10[longitude,latitude,*time]=v10out:v10_attr; ' $
-                      + 'longitude[]=xlon:lon_attr; ' $
-                      + 'latitude[]=ylat:lat_attr; ' $
-                      + 'time[]=timein:time_attr ' $
-                      + ' @ globattr'
++ 'v10[longitude,latitude,*time]=v10out:v10_attr; ' $
++ 'longitude[]=xlon:lon_attr; ' $
++ 'latitude[]=ylat:lat_attr; ' $
++ 'time[]=timein:time_attr ' $
++ ' @ globattr'
 ;
 @ncdf_quickwrite