[MITgcm-devel] netcdf output, gm and visbeck

[Ed Hill]

On Wed, 24 Oct 2007 17:46:57 +0200 Martin Losch <Martin.Losch at awi.de>
wrote:

> A question for the netcdf guru(s).
> I find that the 2D gmredi-variable "Visbeck_K_T" is written to a 3D  
> netcdf field VisbK in in gmredi_output.F. I know that this (tave) is  
> not really the preferred mode of output anymore, but how can I fix  
> this small bug?


Hi Martin,

I think the following diff will fix it.  I think you want to specify
the vertical location to be none ("-") instead of centered ("C") and
tell MNC that the tile indices occupy the 3rd and 4th (not the 4th and
5th) slots.

I don't have the time (or the input files, for that that matter!) to
test it right now.  If it works, that's cool.  And if not, please let
me know what error you get and I'll try again.

Ed

Index: gmredi/gmredi_mnc_init.F
===================================================================
RCS file: /u/gcmpack/MITgcm/pkg/gmredi/gmredi_mnc_init.F,v
retrieving revision 1.1
diff -u -p -u -r1.1 gmredi_mnc_init.F
--- gmredi/gmredi_mnc_init.F    29 Jul 2005 18:50:00 -0000      1.1
+++ gmredi/gmredi_mnc_init.F    24 Oct 2007 22:29:05 -0000
@@ -73,7 +73,7 @@ C    &       'coordinates','XC YC RC ite
      &       'GM Tensor Kwz is K_33 element, Z direction at W point',
      &       myThid)
 
-        CALL MNC_CW_ADD_VNAME('VisbK','Cen_xy_Hn__C__t',4,5,myThid)
+        CALL MNC_CW_ADD_VNAME('VisbK','Cen_xy_Hn__-__t',3,4,myThid)
         CALL MNC_CW_ADD_VATTR_TEXT('VisbK','units','m^2/s', myThid)
         CALL MNC_CW_ADD_VATTR_TEXT('VisbK','description',
      &       'GM mixing/stirring coef after Visbeck et al.',


-- 
Edward H. Hill III, PhD  |  ed at eh3.com  |  http://eh3.com/
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 189 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20071024/dbdcb11c/attachment.sig>