[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