[MITgcm-support] MITgcm-support Digest, Vol 104, Issue 1

liujinliang.627 liujinliang.627 at gmail.com
Wed Feb 1 19:24:01 EST 2012


Hi  Martin,
    Thanks for the answers ! It does help a lot. I read the documentations,  but there are many problems that I can not understand, so I sent an email to you. And I think your answers are very understandable.
  
    Thanks for your help !

Liu Jinliang




At 2012-02-02 01:00:01,mitgcm-support-request at mitgcm.org wrote:
>Send MITgcm-support mailing list submissions to
>	mitgcm-support at mitgcm.org
>
>To subscribe or unsubscribe via the World Wide Web, visit
>	http://mitgcm.org/mailman/listinfo/mitgcm-support
>or, via email, send a message with subject or body 'help' to
>	mitgcm-support-request at mitgcm.org
>
>You can reach the person managing the list at
>	mitgcm-support-owner at mitgcm.org
>
>When replying, please edit your Subject line so it is more specific
>than "Re: Contents of MITgcm-support digest..."
>
>
>Today's Topics:
>
>   1. Format of initialization data and relaxation data
>      (liujinliang.627)
>   2. Help: Format of initialization data and relaxation	data
>      (liujinliang.627)
>   3. Re: Format of initialization data and relaxation	data
>      (Martin Losch)
>
>
>----------------------------------------------------------------------
>
>Message: 1
>Date: Wed, 1 Feb 2012 15:47:29 +0800 (CST)
>From: "liujinliang.627" <liujinliang.627 at gmail.com>
>To: mitgcm-support <mitgcm-support at mitgcm.org>
>Subject: [MITgcm-support] Format of initialization data and relaxation
>	data
>Message-ID:
>	<16638319.7c57.13537e1f62e.Coremail.liujinliang.627 at gmail.com>
>Content-Type: text/plain; charset="gbk"
>
> Hello,
>
>    I'm trying to have a 10 years run for climatology with the example "\verification\tutorial_global_oce_latlon".
>
>1.  But I don't know how to set initial tracer data "hydrogThetaFile" and "hydrogSaltFile".
>  
>     If I want to run for 10 years(Global Ocean Simulation, 1 degree * 1 degree, Nx=360, Ny=160,Nr=15 ), which format should I use for "lev_t.bin" and "lev_s.bin"?  Can it be in an (x,y) or (x,y,r) fashion or  others ?  which software should be use for input, matlab or fortran?
>
>2.  Can the relaxation term be turned off ?  If could, which part should I change?
>
>3.  If I want to using wind force for 10 years continuously (externForcingPeriod = 1month, externForcingCycle = 10years ), what format should "zonalWindFile" and  "meridWindFile" be ?  Should it be in  (Nx,Ny) from first month to the end( i.e. (Nx*12*10, Ny) ), or in (Nx, Ny, t) (t is the number of months) fashion?
>
>Thank you !
>
>Liu Jinliang
>
> 
>-------------- next part --------------
>An HTML attachment was scrubbed...
>URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120201/3e89c84a/attachment.html>
>
>------------------------------
>
>Message: 2
>Date: Wed, 1 Feb 2012 16:01:44 +0800 (CST)
>From: "liujinliang.627" <liujinliang.627 at gmail.com>
>To: mitgcm-support at mitgcm.org, 	"mitgcm-support-request at mitgcm.org"
>	<mitgcm-support-request at mitgcm.org>
>Subject: [MITgcm-support] Help: Format of initialization data and
>	relaxation	data
>Message-ID:
>	<7048cda3.813a.13537ef0192.Coremail.liujinliang.627 at gmail.com>
>Content-Type: text/plain; charset="gbk"
>
> Hello,
>
>    I'm trying to have a 10 years run for climatology with the example "\verification\tutorial_global_oce_latlon".
>
>1.  But I don't know how to set initial tracer data "hydrogThetaFile" and "hydrogSaltFile".
>  
>     If I want to run for 10 years(Global Ocean Simulation, 1 degree * 1 degree, Nx=360, Ny=160,Nr=15 ), which format should I use for "lev_t.bin" and "lev_s.bin"?  Can it be in an (x,y) or (x,y,r) fashion or  others ?  which software should be use for input, matlab or fortran?
>
>2.  Can the relaxation term be turned off ?  If could, which part should I change?
>
>3.  If I want to using wind force for 10 years continuously (externForcingPeriod = 1month, externForcingCycle = 10years ), what format should "zonalWindFile" and  "meridWindFile" be ?  Should it be in  (Nx,Ny) from first month to the end( i.e. (Nx*12*10, Ny) ), or in (Nx, Ny, t) (t is the number of months) fashion?
>
>Thank you !
>
>Liu Jinliang
>
>
>-------------- next part --------------
>An HTML attachment was scrubbed...
>URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120201/5bcb73bd/attachment-0001.htm>
>
>------------------------------
>
>Message: 3
>Date: Wed, 01 Feb 2012 12:26:50 +0100
>From: Martin Losch <Martin.Losch at awi.de>
>To: mitgcm-support at mitgcm.org
>Subject: Re: [MITgcm-support] Format of initialization data and
>	relaxation	data
>Message-ID: <5F62DE45-B6A7-439E-9BD2-D2F856C90A09 at awi.de>
>Content-Type: text/plain; CHARSET=US-ASCII
>
>Hi Liu Jinliang,
>
>see answers below. (And I really recommend reading the documenation.)
>
>Martin
>On Feb 1, 2012, at 8:47 AM, liujinliang.627 wrote:
>
>>  Hello,
>> 
>>     I'm trying to have a 10 years run for climatology with the example "\verification\tutorial_global_oce_latlon".
>> 
>> 1.  But I don't know how to set initial tracer data "hydrogThetaFile" and "hydrogSaltFile". 
>>    
>>      If I want to run for 10 years(Global Ocean Simulation, 1 degree * 1 degree, Nx=360, Ny=160,Nr=15 ), which format should I use for "lev_t.bin" and "lev_s.bin"?  Can it be in an (x,y) or (x,y,r) fashion or  others ?  which software should be use for input, matlab or fortran? 
>hydrogThetaFile/hydrogSaltFile require 3D fields (x,y,r). 
>The format the all input files (including hydrogThetaFile/hydrogSaltFile) is ieee-be. You can use any tool to generate such files. In matlab it is very simple:
>t = zeros(nx,ny,nr); % here you need to use realistic values of course
>fid = fopen('theta.init','w','b');
>fwrite(fid,t,'real*4'); % readBinaryPrec=32, for readBinaryPrec=64 use 'real*8'
>fclose(fid).
>In Fortran you have to open your file with ACCESS='DIRECT' (and make sure that you write the correct precision and endian type)
>> 
>> 2.  Can the relaxation term be turned off ?  If could, which part should I change?
>yes. In your example like this (no relaxation for either temperature nor salinity):
> tauThetaClimRelax=  0.,
> tauSaltClimRelax = 0.,
>> 
>> 3.  If I want to using wind force for 10 years continuously (externForcingPeriod = 1month, externForcingCycle = 10years ), what format should "zonalWindFile" and  "meridWindFile" be ?  Should it be in  (Nx,Ny) from first month to the end( i.e. (Nx*12*10, Ny) ), or in (Nx, Ny, t) (t is the number of months) fashion?
>If you want to precribe monthly data for 10 years (and all years are different), then your data needs to have  the format (Nx,Ny,10*12) (and ieee-be). But setting these flags requires that all forcing data is in this format, i.e., also for these
> thetaClimFile=  'lev_sst.bin',
> saltClimFile=   'lev_sss.bin',
> surfQFile=      'ncep_qnet.bin',
> EmPmRFile=      'ncep_emp.bin',
>
>To have more freedom (and more problems) use the exf-package. 
>> 
>> Thank you !
>> 
>> Liu Jinliang
>>  
>> 
>> 
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>
>
>
>------------------------------
>
>_______________________________________________
>MITgcm-support mailing list
>MITgcm-support at mitgcm.org
>http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>
>End of MITgcm-support Digest, Vol 104, Issue 1
>**********************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120202/232f17e2/attachment.htm>


More information about the MITgcm-support mailing list