[MITgcm-support] Global_ocean query

Martin Losch mlosch at awi-bremerhaven.de
Mon Mar 21 02:50:58 EST 2005


Amitabh,

I don't know what went wrong exactly, but clearly you got something 
wrong with the grid. To me it looks like you have somehow scrambled 
your input fields: from the second plot (Qdiag and (E-P)diag, it is 
pretty clear that you have shifted the input sst and sss fields by 1 or 
2 per horizontal row (latitude), so your field have maybe 181 instead 
of 180 entries per row? (Qdiag is computed from the different of T at 
the surface and whatever you have specified at surface restoring field, 
thetaClimField, in data). To check this, try this in matlab:

fname='your2dinputfield';
fid=fopen(fname,'r','b');
t=fread(fid,[nx ny nt],'real*4'); % or 'real*8', whatever you use, 
nx,ny,nt should be the dimension of your field (you should know them, 
as you have created these files)
fclose(fid);
% then plot the fields
t(find(t==0))=NaN;
pcolor(t(:,:,1)'); shading flat
% if that doesn't look like the earth, but like the warped map that you 
see on you second figure, then you have to redo your input fields. Make 
sure that they have exactly the same dimensions as in SIZE.h, also save 
them like this:
t=your field(nx,ny,nt);
fid=fopen('your2dinputfield','w','b');fwrite(fid,t,'real*4/8'); 
fclose(fid);
% that should work. Do this for all of your input fields:
surfThetaClim = SST,
surfSaltClim = SSS,
surfQnetFile = QNET,
surfQswFile = shortwave heat,
EmPmRFile = EMP
zonalWindFile = TAUX,
meridWindFile =  TAUY,
ploadFile = SURFACE Pressure,
see PARAMS.h for all possible files and FFIELDS.h for typical ranges 
(and sign convention!!!) of the fields.

Martin

On Mar 18, 2005, at 6:47 PM, Ed Hill wrote:

>
>> -------- Forwarded Message --------
>>> From: Amitabh  Mitra <amitabhmitra1 at rediffmail.com>
>>> Reply-To: Amitabh  Mitra <amitabhmitra1 at rediffmail.com>
>>> To: mitgcm-support at mitgcm.org
>>> Subject: Re: Re: [MITgcm-support] Global_ocean query
>>> Date: 18 Mar 2005 13:19:13 -0000
>
>>> Dear Sir,
>>>
>>> I plotted the model output
>>>
>>> using matlab. I made certain
>>>
>>> changes in mit_oceanmasks.m. The ps file is
>>>
>>> attached along with the mail.
>>>
>>> The map is not correct. I
>>>
>>> suppose something is wrong
>>>
>>> with the input.
>>>
>>>   Please can you state where have I
>>>
>>> gone wrong.
>>>
>>> Thanking You for your kind support,
>>>
>>> Amitabh Mitra




More information about the MITgcm-support mailing list