[MITgcm-support] MITgcm-support Digest, Vol 126, Issue 44
katsman
katsman at knmi.nl
Mon Jan 6 09:55:36 EST 2014
Dear Jean-Michel,
The files are attached.
Caroline
On 12/24/2013 06:00 PM, 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. Re: using RBCS to generate flow (Jean-Michel Campin)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Tue, 24 Dec 2013 09:26:52 -0500
> From: Jean-Michel Campin <jmc at ocean.mit.edu>
> To: mitgcm-support at mitgcm.org
> Subject: Re: [MITgcm-support] using RBCS to generate flow
> Message-ID: <20131224142652.GB19034 at ocean.mit.edu>
> Content-Type: text/plain; charset=us-ascii
>
> Hi Caroline,
>
> Could you send few parameter files that you are using when trying to
> work with checkpoint64 code ?
> it would be "data", "data.pkg" and "data.rbcs".
>
> Cheers,
> Jean-Michel
>
> On Tue, Dec 24, 2013 at 01:54:37PM +0100, katsman wrote:
>> Dear MITgcm developers,
>>
>> Over the past years, I used the RBCS package to generate a flow in
>> an (otherwise unforced) basin representing the Labrador Sea, by
>> restoring the flow to a prescibed 3D-temperature and velocity field
>> in a corner of the basin (T and velocity in geostrophic balance).
>> This worked excellent in a (admittedly old) checkpoint58-version in
>> which we manually added a restoring term on U and V analogous to the
>> programmed T, S part.
>>
>> When I got a new workstation I thought it time to upgrade to the
>> most recent MITgcm version, but I cannot get the package to work in
>> the same way (and I tried many things over the past months, so I am
>> getting a bit desperate...).
>>
>> The new checkpoint64 is supposed to be able to restore T,U and V,
>> but I find that - when I prescribe my T,V fields again as before -
>> in the sponge region defined by the mask, T is restored fine, while
>> V is restored fine only during the first few timesteps of the
>> simulation. Then a purely baroclinic U develops as well, while V
>> becomes barotropic (despite RBCS acting on it). Intriguing, but not
>> what I wanted.
>>
>> A test that really makes me think there is a bug in the rbcs package
>> rather than me doing something stupid is that when I only restore T
>> with a gradient in x, I expect to see a purely baroclinic V develop
>> in geostrophic balance with that. In stead, nothing happens (I can
>> run it for a month, no flow develops at all).
>>
>> One of the earlier versions (c62) has yet another version of the
>> RBCS package. When I do a test with T-restore only (U,V is not
>> standard in that one), a baroclinic flow does develop obeying
>> geostrophy.
>>
>> Any ideas what is wrong here? I checked - the model does frequent
>> the appropriate lines of code in rbcs_add_tendency in the c64
>> version to change the gV parameter - I suspect somewhere in the
>> code U and V are mixed up. Notably, the baroclinic flow that
>> develops in the sponge region is roughly the same strength as the V
>> I prescribe
>>
>> Any ideas?? (if the description is unclear I can send pictures)
>>
>> Thanks in advance & happy holidays
>> Caroline
>>
>>
>>
>>
>> _______________________________________________
>> 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 126, Issue 44
> ***********************************************
-------------- next part --------------
#; ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
&PARM01
# tRef= 5.367, 4.102, 3.628, 3.447, 3.341,
# 3.270, 3.197, 3.156, 3.133, 3.115,
# 3.096, 3.050, 2.947, 2.781, 2.655,
# sRef= 15*0.0,
# normal
tRef = 5.36700,4.10200,3.62800,3.44700,3.34100,
3.27710,3.21525,3.17308,3.14642,3.12914,
3.11629,3.10382,3.08450,3.05690,3.00880,
2.95112,2.87507,2.79207,2.72171,2.63647,
# thick upper layer - changed T2+T3
# tRef = 5.36700,3.92425,3.60538,3.44700,3.34100,
# 3.27710,3.21525,3.17308,3.14642,3.12914,
# 3.11629,3.10382,3.08450,3.05690,3.00880,
# 2.95112,2.87507,2.79207,2.72171,2.63647,
sRef= 20*35.0,
viscAz=1.E-5,
# viscAh=40.,
viscA4=1.5e09,
no_slip_sides=.TRUE.,
no_slip_bottom=.TRUE.,
# diffKhT=20.,
diffK4T=0.75e09,
diffKzT=1.0e-05,
diffKhS=0.0,
diffKzS=0.0,
f0=1.16e-4,
beta=1.4e-11,
tAlpha=1.7e-4,
sBeta =0.0,
#gravity=9.81,
rigidLid=.FALSE.,
implicitFreeSurface=.TRUE.,
#useCDscheme=.TRUE.,
#useNHMTerms=.TRUE.,
eosType='LINEAR',
rhoNil=1028.,
hFacMin=0.05,
hFacMinDz=50,
#nonHydrostatic=.TRUE.,
readBinaryPrec=64,
bottomDragLinear=0.E-4,
exactConserv=.TRUE.,
&
# Elliptic solver parameters
&PARM02
cg2dMaxIters=300,
cg2dTargetResidual=1.E-13,
cg3dMaxIters=20,
cg3dTargetResidual=1.E-8,
&
# Time stepping parameters
&PARM03
startTime=0,
# 24 hrs
endTime=86400,
# 3 days
# endTime=259200,
deltaT=600,
deltaTtracer=600,
#deltaTClock =86400.0,
abEps=0.01,
# pChkptFreq=000,
chkptFreq=2592000,
# 1 day
dumpFreq=600,
# 3 day
# dumpFreq=259200,
# dumpFreq=3600,
cadjFreq=600,
monitorFreq=1.,
periodicExternalForcing=.FALSE.,
# periodicExternalForcing=.TRUE.,
# 1 / 12 months
# externForcingPeriod=2592000.,
# externForcingCycle=31104000.,
&
# Gridding parameters
&PARM04
usingCartesianGrid=.TRUE.,
# delZ= 100, 100, 100, 100, 100,
# 150, 150, 150, 150, 200,
# 225, 275, 350, 400, 450,
#
# # normal
delZ= 100,100,100,100,100,
125,125,125,125,125,
125,150,150,150,175,
175,175,200,250,325,
# thick upper layer
# delZ= 150,75,75,100,100,
# 125,125,125,125,125,
# 125,150,150,150,175,
# 175,175,200,250,325,
delX=162*7.5e03,
delY=210*7.5e03,
&
# Input datasets
&PARM05
# bathyFile = '../INPUT_OCT/topogLS_rbcs_narrow_dijk_steep_7500_64.bin',
bathyFile = '../INPUT_OCT/topogLS_rbcs_narrow_dijk_steep_7500_64.bin',
# surfQfile = '../INPUT_OCT/Qfile_100_7.5km.bin',
#
# bathyFile = '../INPUT_IDL/topogLS_rbcs_09_dijk_64.bin',
# surfQfile = '../INPUT_IDL/Qfile_107_7.5km_64.bin',
&
-------------- next part --------------
# RBCS package parameters:
&RBCS_PARM01
#------------------------------------------------------------------------------
# switches
#------------------------------------------------------------------------------
useRBCtemp=.TRUE.,
useRBCvVel=.TRUE.,
useRBCuVel=.TRUE.,
# useRBCsalt=.TRUE.,
#------------------------------------------------------------------------------
#- relaxation times
#------------------------------------------------------------------------------
tauRelaxU=600.,
tauRelaxV=600.,
tauRelaxT=600.,
#------------------------------------------------------------------------------
# masks - #1 and #2 = T,S , #3 = tracers ;
# masks U,V have explicit name - if left out then TMask = used
#------------------------------------------------------------------------------
# -----OCTAVE VERSION----------------------------
relaxMaskFile='../INPUT_OCT/TempMask_7500_north_64_z20.bin','../INPUT_OCT/SaltMask_7500_north_64_z20.bin',
relaxMaskUFile='../INPUT_OCT/UvelMask_7500_north_64_z20.bin',
relaxMaskVFile='../INPUT_OCT/VvelMask_7500_northeast_64_z20.bin',
#
# -----IDL VERSION
# relaxMaskFile='../INPUT_IDL/TempMask_75_east_64_z20.bin','../INPUT_IDL/SaltMask_75_east_64_z20.bin',
# relaxMaskVFile='../INPUT_IDL/VMask_75_east_64_z20.bin',
#------------------------------------------------------------------------------
# files containing relaxation flds have explicit name
#------------------------------------------------------------------------------
# -----OCTAVE VERSION
relaxTFile='../INPUT_OCT/TempFile_7500_drho26_seas055_north_64_z20.bin',
relaxUFile='../INPUT_OCT/UvelFile_7500_drho26_seas055_north_64_z20.bin',
relaxVFile='../INPUT_OCT/VvelFile_7500_drho26_seas055_north_64_z20.bin',
# relaxTFile='../INPUT_OCT/TempFile_7500_drho26_fix_northeast_64_z20.bin',
# relaxVFile='../INPUT_OCT/VvelFile_7500_drho26_fix_northeast_64_z20.bin',
# -----IDL VERSION
# relaxTFile='../INPUT_IDL/TempFile_75_drho26_seas100_east_64_z20.bin',
# relaxVFile='../INPUT_IDL/VFile_75_drho26_seas100_east_64_z20.bin',
#------------------------------------------------------------------------------
rbcsIniter=0,
# 1 / 12 months
# rbcsForcingPeriod=2592000.,
# rbcsForcingCycle=31104000.,
# 5 / 60 days
rbcsForcingPeriod=432000.,
rbcsForcingCycle=5184000.,
&
-------------- next part --------------
A non-text attachment was scrubbed...
Name: data.pkg
Type: application/vnd.apple.installer+xml
Size: 57 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20140106/e47a177d/attachment.pkg>
More information about the MITgcm-support
mailing list