[MITgcm-support] [SPAM] Re: [SPAM] Re: Transport through a section

michael schaferkotter schaferk at bellsouth.net
Tue Mar 8 17:57:30 EST 2016


you probably want to look at the values of hFacMin, and  hFacMinDr

it seems a bit strange that hFac[WCS] ==0.


have a look at http://mitgcm.org/public/r2_manual/latest/online_documents/node49.html


On Mar 7, 2016, at 8:21 PM, lianzhan at fio.org.cn wrote:

> In my case, I have turned off the partial cells and I have check the hFacW, hFacC and hFacS in the results. All of these variables are zero. 
> 
> Should I consider the changing of thickness in the surface layer caused by eta (free-surface_r-anomaly)?
> 
>>>> With my best regards
>>>> 
>>>> Lian Zhan
> 
> First Institute of Oceanography (FIO),State Oceanic Administration (SOA),P.R.China
> Add: No.6 Xianxialing Road,Qingdao,266061,P.R.China
> Tel: +86-532-88967320
> E-mail: lianzhan at fio.org.cn
>  
> From: michael schaferkotter
> Date: 2016-03-08 07:08
> To: mitgcm-support
> Subject: [SPAM] Re: [MITgcm-support] [SPAM] Re: Transport through a section
> excerpt of fortran90 code
> 
>       allocate(OBWmask(NNY,NR))
> 
>       OBWmask             = hFacW(  2,  :,:) 
>       OBWmask(  1,:)    = 0 
>       OBWmask(NNY,:) = 0 
> 
> !***********************************************************************
> !  compute area of each face
> !***********************************************************************
> 
>       do k=1,NR
>         OBWmask(:,k) = thk(k)*OBWmask(:,k)*DYG(2  ,:  )
>        enddo
> 
> !thk(k) is thickness of cells, found in RF.data
> 
>       allocate(uobw(NNY,NR))
> 
>       ! read OBWuFile into uobw for timestep with a routine
> 
>         OBW = sum(uobw*OBWmask)
> 
> 
> 
> On Mar 7, 2016, at 8:18 AM, Ryan Abernathey wrote:
> 
>> You need to multiply by hFacW in order to account for partial cells
>> 
>> Sent from my iPhone
>> 
>> On Mar 7, 2016, at 8:54 AM, "lianzhan at fio.org.cn" <lianzhan at fio.org.cn> wrote:
>> 
>>> Thank you for your reminder. I have made a mistake in last mail, actually I did caculate transport by v(k)*dy(k)*dz(k). And furthermore, the dy is uniform in this case, so I think it is not the key of problem. 
>>> 
>>>> With my best regards
>>>> 
>>>> Lian Zhan
>>> 
>>> First Institute of Oceanography (FIO),State Oceanic Administration (SOA),P.R.China
>>> Add: No.6 Xianxialing Road,Qingdao,266061,P.R.China
>>> Tel: +86-532-88967320
>>> E-mail: lianzhan at fio.org.cn
>>>  
>>> From: Michael Schaferkotter
>>> Date: 2016-03-07 19:37
>>> To: mitgcm-support
>>> Subject: [SPAM] Re: [MITgcm-support] Transport through a section
>>> transport will be velocity times area of face, presumably in this case \sum_{k=1}^km v(k)*dy(k)*dz(k).
>>> 
>>> Sent from Here3.
>>> 
>>> On Mar 7, 2016, at 04:03, "lianzhan at fio.org.cn" <lianzhan at fio.org.cn> wrote:
>>> 
>>>>>>>>> Dear MITgcm users,
>>>> 
>>>> 
>>>> I am simulating the  barotropic gyre in an ideal rectangular basin, there is an island in this basin and there is a sill on the left side of the island (the sketch map is shown in attachment). I want to study the transports of water through section A and B (in the north and south of the strait with sill). There two sections , the western boundary of island and the weastern boundary of basin form the boundary of an closed region, therefore the transports through these two sections are suppose to be strictly equal. But they are not, as seen from the result of numerical model. Is it possibly and how could this be? Every suggestion would be greatly appreciated!
>>>> 
>>>> I calculate the transport <Catch5AB7.jpg>by: <Catch.jpg>, where the <Catch3834.jpg> is meridinal velocity and dz is the thickness of layers. The number of layers in this model is "km".  And here is my "data" file:
>>>> 
>>>> 
>>>>  &PARM01
>>>>  tRef=24*20.,
>>>>  sRef=24*10.,
>>>>  viscAz=1.E-4,
>>>>  viscAh=2.E4,
>>>>  diffKhT=1.E3,
>>>>  diffKzT=1.E-4,
>>>>  f0=0.2E-4,
>>>>  beta=2.E-11,
>>>>  tAlpha=2.E-4,
>>>>  sBeta =0.,
>>>>  gravity=9.81,
>>>>  gBaro=9.81,
>>>>  rigidLid=.FALSE.,
>>>>  implicitFreeSurface=.TRUE.,
>>>>  rhonil=1035.,
>>>>  eosType='LINEAR',
>>>>  readBinaryPrec=32,
>>>>  no_slip_sides=.FALSE.,
>>>>  no_slip_bottom=.FALSE.,
>>>>  &
>>>> 
>>>> # Elliptic solver parameters
>>>>  &PARM02
>>>>  cg2dMaxIters=1000,
>>>>  cg2dTargetResidual=1.E-7,
>>>>  &
>>>> 
>>>> # Time stepping parameters
>>>>  &PARM03
>>>>  startTime=0,
>>>>  endTime=311040000,
>>>>  deltaTmom=200.0,
>>>> # tauCD=172800.,
>>>>  abEps=0.1,
>>>>  pChkptFreq=311040000.0,
>>>>  chkptFreq=311040000.0,
>>>>  dumpFreq=31104000.0,
>>>>  monitorFreq=311040000.,
>>>>  &
>>>> 
>>>> # Gridding parameters
>>>>  &PARM04
>>>>  usingCartesianGrid=.TRUE.,
>>>>  delX=500*20E3,
>>>>  delY=250*20E3, 
>>>>  delZ=5.,5.,5.,5.,10.,10.,10.,10.,20.,20.,20.,20.,50.,50.,50.,50.,100.,100.,100.,100.,100.,100.,100.,100.,
>>>>  &
>>>> 
>>>> # Input datasets
>>>>  &PARM05
>>>>  bathyFile='dep_1i_sill_50.bin',
>>>>  zonalWindFile='taux01.bin',
>>>> 
>>>> With my best regards
>>>> 
>>>> Lian Zhan
>>>> 
>>>> First Institute of Oceanography (FIO),State Oceanic Administration (SOA),P.R.China
>>>> Add: No.6 Xianxialing Road,Qingdao,266061,P.R.China
>>>> Tel: +86-532-88967320
>>>> E-mail: lianzhan at fio.org.cn
>>>> <dep.gif>
>>>> _______________________________________________
>>>> 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
>> _______________________________________________
>> 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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160308/816a34aa/attachment-0001.htm>


More information about the MITgcm-support mailing list