[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