Changeset 205 for trunk/procs/macros/make_energetics.pro

Timestamp:

01/26/10 10:46:13 (14 years ago)

Author:

pinsard

Message:

homegenize THEN BEGIN ... ENDIF

File:

• trunk/procs/macros/make_energetics.pro (modified) (9 diffs)

trunk/procs/macros/make_energetics.pro

@@ -190 +190 @@
          valmask = 1.e20
          idxw = where (tmaskr EQ 0)
-         IF idxw[0] NE -1 THEN w[idxw] = valmask
+         IF idxw[0] NE -1 THEN BEGIN
+          w[idxw] = valmask
+         ENDIF
          idxt = where (tmaskr EQ 0)
-         IF idxt[0] NE -1 THEN t [idxt] = 0.
-         IF idxt[0] NE -1 THEN s [idxt] = 0.
+         IF idxt[0] NE -1 THEN BEGIN
+          t[idxt] = 0.
+         ENDIF
+         IF idxt[0] NE -1 THEN BEGIN
+          s[idxt] = 0.
+         ENDIF
          tmaskr = 0 ; free memory
       ENDIF ELSE BEGIN
          idxt = where (t EQ valmask)
-         IF idxt[0] NE -1 THEN t [idxt] = 0.
-         IF idxt[0] NE -1 THEN s [idxt] = 0.
+         IF idxt[0] NE -1 THEN BEGIN
+          t[idxt] = 0.
+         ENDIF
+         IF idxt[0] NE -1 THEN BEGIN
+          s[idxt] = 0.
+         ENDIF
       ENDELSE
 
@@ -223 +233 @@
 
       idx_2d = where (u[*, *, 0, 0] GT valmask/10.)
-      IF idx_2d[0] NE -1 THEN tx[idx_2d] = valmask
+      IF idx_2d[0] NE -1 THEN BEGIN
+       tx[idx_2d] = valmask
+      ENDIF
       idx_2d = where (v[*, *, 0, 0] GT valmask/10.)
-      IF idx_2d[0] NE -1 THEN ty[idx_2d] = valmask
+      IF idx_2d[0] NE -1 THEN BEGIN
+       ty[idx_2d] = valmask
+      ENDIF
       idx = where (t LT valmask/10.)
       IF t_unit NE "C" THEN BEGIN
-         IF idx[0] NE -1 THEN t[idx] = t[idx]-273.15
+         IF idx[0] NE -1 THEN BEGIN
+          t[idx] = t[idx]-273.15
+         ENDIF
       ENDIF
       idx1 = n_elements(idx_t)
@@ -241 +257 @@
       idxt=where(t GT valmask/10.)
       idxs=where(s GT valmask/10.)
-      IF source_model EQ 'ipcc' THEN print, '       check these are 0 (T,S): ', min(idxt -idxs), max(idxt-idxs)
-      IF idxt[0] NE -1 THEN t[idxt]=0.
-      IF idxs[0] NE -1 THEN s[idxs]=0.
+      IF source_model EQ 'ipcc' THEN BEGIN
+       print, '       check these are 0 (T,S): ', min(idxt -idxs), max(idxt-idxs)
+      ENDIF
+      IF idxt[0] NE -1 THEN BEGIN
+       t[idxt]=0.
+      ENDIF
+      IF idxs[0] NE -1 THEN BEGIN
+       s[idxs]=0.
+      ENDIF
       idxs = 0                  ; free memory
 
@@ -262 +284 @@
    rhop = ( ( 4.8314E-4*s + r3*sr +r2)*s +r1)
    print, '      rhop min/max', min(rhop), max(rhop)
-   IF idxt[0] NE -1 THEN rhop[idxt] = valmask
+   IF idxt[0] NE -1 THEN BEGIN
+    rhop[idxt] = valmask
+   ENDIF
    sr = 0 & r1 = 0 &  r2 = 0 &  r3 = 0 & t = 0 & s = 0 ; free memory
 
@@ -305 +329 @@
 
    int_val2 = ((rhop-rho_s4d)^2)*stab_inv
-   IF idxt[0] NE -1 THEN int_val2[idxt] = 0.
+   IF idxt[0] NE -1 THEN BEGIN
+    int_val2[idxt] = 0.
+   ENDIF
 
    print, '     compute APE...'
@@ -321 +347 @@
 
    int_val = (rhop-rho_s4d)*(w_T)
-   IF idxt[0] NE -1 THEN int_val[idxt] = 0.
+   IF idxt[0] NE -1 THEN BEGIN
+    int_val[idxt] = 0.
+   ENDIF
 
 ; remove first 2 levels (MXL too unstable)
@@ -348 +376 @@
    idyv = where(vmean GT valmask/10.)
 
-   IF idx[0] NE -1 THEN tx[idx] = 0.
-   IF idy[0] NE -1 THEN ty[idy] = 0.
-   IF idxu[0] NE -1 THEN umean[idxu] = 0.
-   IF idyv[0] NE -1 THEN vmean[idyv] = 0.
+   IF idx[0] NE -1 THEN BEGIN
+    tx[idx] = 0.
+   ENDIF
+   IF idy[0] NE -1 THEN BEGIN
+    ty[idy] = 0.
+   ENDIF
+   IF idxu[0] NE -1 THEN BEGIN
+    umean[idxu] = 0.
+   ENDIF
+   IF idyv[0] NE -1 THEN BEGIN
+    vmean[idyv] = 0.
+   ENDIF
 
    dot_prodx = tx*umean
@@ -387 +423 @@
    tvlct, red, green, blue
 
-   IF cmd_wrk.out EQ 'ps' THEN ps = 1
-
-   IF ps EQ 1 THEN openps
+   IF cmd_wrk.out EQ 'ps' THEN BEGIN
+    ps = 1
+   ENDIF
+
+   IF ps EQ 1 THEN BEGIN
+    openps
+   ENDIF
 
    pltt, ape, 't', petit = [2, 4, 1], landscape = 1, /rempli, /BASICMARGES, title = 'APE (full) '+cmd_wrk.exp, window=3
@@ -422 +462 @@
 
    maxdepth = 10
-   IF gdept[0] GT maxdepth THEN maxdepth = gdept[0]
+   IF gdept[0] GT maxdepth THEN BEGIN
+    maxdepth = gdept[0]
+   ENDIF
 
    ;domdef, [210., 280., -5., 5., 0., maxdepth]