[Thread Prev][Thread Next][Index]

Re: Heat content in the mixed layer

Hi Satheesh,
For the question of how to do this step

   let temp_vertint = rhow*Cp*temp[z=0:`mix_layer_h`@din]

I think you can use the ideas in this FAQ: "Averaging/integrating over
irregular regions in Ferret", at

You would create a mask using the mixed layer depth and multiply
the variable you are integrating by the mask.  Then everything below
that depth is missing and therefore not included in the integration.
After you have defined mix_layer_h:

   LET mask = IF z LE mix_layer_h THEN 1
   LET temp_vertint = rhow*Cp*mask* temp[z=@din]

Ansley Manke

"Satheesh C. Shenoi" wrote:

> Dear Ferret Guru,
> Perhaps this is a trivial question, but I would like to get some help on
> this.
> I want to compute the heat content in the mixed layer. The mixed layer
> height is computed as follows:
> let mix_layer_dens = dens[k=1]+0.2
> let zero_at_mix_layer_dens = dens - mix_layer_dens
> let mix_layer_h = zero_at_mix_layer_dens[z=@loc:0]
> mix_layer_h has X,Y and T (time) dimensions. Now to compute the heat
> content
> in the water column above the mixed layer. I tried the following
> let temp_vertint = rhow*Cp*temp[z=0:`mix_layer_h`@din]
> it looks, Ferret is unable to substitute for `mix_layer_h`. I am using
> Ferret
> version 5.40. I am getting following error:
> let temp_vertint = rhow*Cp*temp[z=0:`mix_layer_h`@din]
>  **ERROR: invalid command: grave accent doesnt evaluate to scalar
>           mix_layer_h
> DEFINE VARIABLE temp_vertint = rhow*Cp*temp[z=0:`mix_layer_h`@din]
> Command file, command group, or REPEAT execution aborted.
> Please let me know how to compute the heat content in the mixed layer.
> Thanks in advance.
> -satheesh
> Satheesh C. Shenoi                     Ph.:   (91)(832) 226253-56 ext4224/4400
> Physical Oceanography Division,        Fax:   (91)(832) 223340/221360
> National Institute of Oceanography,    email: shenoi@darya.nio.org
>                                               shenoi@csnio.ren.nic.in
> Dona Paula, Goa 403 004, India.        Ph.:   (91)(832) 233453 (home)

Ansley Manke  Pacific Marine Environmental Laboratory  Seattle WA  (206)526-6246

[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement