Changeset 462 for trunk/SRC/Computation/bsf.pro

Timestamp:

04/15/11 19:49:48 (13 years ago)

Author:

smasson

Message:

get time axis even when using /timestep keyword in read_ncdfSRC/ToBeReviewed/LECTURE/read_ncdf.pro

File:

• trunk/SRC/Computation/bsf.pro (modified) (3 diffs)

trunk/SRC/Computation/bsf.pro

@@ -23 +23 @@
 ; @keyword REFVALUE {type=scalar} {default=0.}
 ; the bsf value that we want to speficy at the position defined by refpoint
+;
+; @keyword TRANSPORT {type=scalar: 0 or 1} {default=0.}
+; activate to specify that z3d is not a zonal current but a zonal
+; transport (e2u*e3u*un)
 ;
 ; @returns {type=2D xy array}
@@ -59 +63 @@
 ;-
 ;
-FUNCTION bsf, z3d, NOSTRUCTURE = nostructure, REFPOINT = refpoint, REFVALUE = refvalue
+FUNCTION bsf, z3d, NOSTRUCTURE = nostructure, REFPOINT = refpoint, REFVALUE = refvalue, TRANSPORT = transport
 ;
   compile_opt idl2, strictarrsubs
@@ -81 +85 @@
   un = fitintobox(temporary(un))
 ;
-  e23 = e3u_3d(/e2)
-;
-; mask the array
-  un = temporary(umsk) * temporary(un)
+; current -> transport
+  IF NOT keyword_set(transport) THEN un = temporary(un) * e3u_3d(/e2)
+  utr = temporary(umsk) * temporary(un)
 ; compute the bsf
-  bsf = 1.e-6 * total(total(temporary(un) * temporary(e23), 3), 2, /cumulative)
+  bsf = 1.e-6 * total(total(temporary(utr), 3), 2, /cumulative)
 ; set bsf to 0 in the largest continent... no done...
   IF keyword_set(refpoint) THEN BEGIN