[MITgcm-support] urgent question

Martin Losch Martin.Losch at awi.de
Wed Dec 27 03:43:49 EST 2006


Riema,

I frequently make this mistake myself: The topography needs to be  
negative, that is, it does not contain "depth" (which would be a  
positive number, right?), but vertical position of the bottom z=-H.
Your temp.new has only two values, the model expect 700x2.

MDSIO etc. You should not have to modify anything in these routines.  
But this is how I look for routines. Change into your "build"  
directory and use grep:

csysm3::build> grep -i 'ine mdswritefield' *.F | grep -v -i call
mdsio_gl.F:      SUBROUTINE MDSWRITEFIELD_3D_GL(
mdsio_gl.F:      SUBROUTINE MDSWRITEFIELD_2D_GL(
mdsio_gl_slice.F:      SUBROUTINE MDSWRITEFIELD_XZ_GL(
mdsio_gl_slice.F:      SUBROUTINE MDSWRITEFIELD_YZ_GL(
mdsio_slice.F:      SUBROUTINE MDSWRITEFIELDXZ(
mdsio_slice.F:      SUBROUTINE MDSWRITEFIELDYZ(
mdsio_slice_loc.F:      SUBROUTINE MDSWRITEFIELDXZ_LOC(
mdsio_slice_loc.F:      SUBROUTINE MDSWRITEFIELDYZ_LOC(
mdsio_writefield.F:      SUBROUTINE MDSWRITEFIELD(            <<<<<<  
this is what you want
mdsio_writefield_loc.F:      SUBROUTINE MDSWRITEFIELD_LOC(
mdsio_writefield_new.F:      SUBROUTINE MDSWRITEFIELD_NEW(

I you want to write (not read or prescribe?) the barotropic tidal  
velocity field along OBW, then you should use MDSWRITEFIELDYZ in  
obcs_calc.F
But I think you want to prescribe this field at every time step. For  
that, you can either use the code that's there (with  
ALLOW_OBCS_PRESCRIBE and useOBCSprescribe, see exp4 for an example  
how to use this stuff), or use MDSREADFIELDYZ.

Martin

On 25 Dec 2006, at 06:49, Riema Rachmayani wrote:

> dear MITgcm,
>
> as Martin told me before that all input files of for mitgcm are  
> binary ieee big endian files, such as bathyfile and hydrogThetaFile  
> in internal wave verification. after i generate to binary ieee big  
> endian files, i still have zero result for the output (i thought  
> because uncorrect while read the input) and sometimes i have error  
> message :
>
> unit 9 named temp.ok lately reading direct unformatted external IO
>
> i was looking where is the statement of input/output files format  
> for example like : OPEN(dUnit,file=fn 
> (1:iLen),status='unknown',form='formatted'), specially for  
> bathyfile and hydrothetafile or others input. is there a routine  
> that show the statement?(formatt)
>
> for example in my new case, i have bathyfile in ascii file  
> (SILL.BTM) and i generate to binary ieee big endian (SILL.OK), and   
> hydrogThetaFile in ascii file (temp.new) ==> binary ieee big endian  
> (temp.ok) in matlab :
>
> load SILL.BTM;
> % write data into big endian format
> fid = fopen(SILL.OK','w','ieee-be'); % open new file in big endian  
> mode
> wrt = fwrite(fid,temp,'float32');       % write data into opened  
> file, using stream and float32 press
> fclose(fid);
>
> load temp.new;
> % write data into big endian format
> fid = fopen('temp.OK','w','ieee-be'); % open new file in big endian  
> mode
> wrt = fwrite(fid,temp,'float32');       % write data into opened  
> file, using stream and float32 press
> fclose(fid);
>
> my bathyfile in ascii file is in row data for Nx=700, Ny=1, and my  
> hidrothetafile in 2 row (2 layers) Nr=2 (you can see in bathy and  
> hydroTheta attachment file). give me some correction if i was wrong  
> to generate the file.
>
> i change nothing in internal wave verification except bathy and  
> hidrothetafile,
> anyone can help me for my case??
>
> thank you MITgcm....
>
> -rima-
>
>
>
> Send instant messages to your online friends http:// 
> uk.messenger.yahoo.com
> <SILL.BTM>
> <temp.new>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list