[MITgcm-support] Top boundary conditions
Martin Losch
Martin.Losch at awi.de
Sat Feb 6 12:25:39 EST 2010
Hi Georgy,
what do you need this for? In "realistic" cases, oceanic surface friction (at the air-sea interface) is realized by imposing wind stress. If you want an ocean with a solid lid on it, you could try (and "abuse") the ice-shelf package (pkg/shelfice, see verificationi/isomip for an example application) and put thin ice-shelp on top of your ocean. Then you can add friction a the ice-ocean interface (have a look at the possible parameters in shelfice_readparms.F)
Martin
On Feb 5, 2010, at 11:37 PM, Georgy Manucharyan wrote:
> Hello,
> For the ocean gcm, is there a way to impose friction at the surface (no slip b.c.)?
> Thanks for your help,
> Georgy.
>
>
> 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: bottom pressure diagnostics (Jean-Michel Campin)
>> 2. Re: bottom pressure diagnostics (Oliver Jahn)
>> 3. Re: bottom pressure diagnostics (Oliver Jahn)
>>
>>
>> ----------------------------------------------------------------------
>>
>> Message: 1
>> Date: Fri, 5 Feb 2010 05:22:12 -0500
>> From: Jean-Michel Campin <jmc at ocean.mit.edu>
>> To: mitgcm-support at mitgcm.org
>> Subject: Re: [MITgcm-support] bottom pressure diagnostics
>> Message-ID: <20100205102212.GA28909 at ocean.mit.edu>
>> Content-Type: text/plain; charset=us-ascii
>>
>> Hi Dimitris,
>>
>> I would assume the same thing as you.
>> I will check when I am back from vacation.
>> Cheers,
>> Jean-Michel
>>
>> On Thu, Feb 04, 2010 at 08:33:07AM -0800, Dimitris Menemenlis wrote:
>>
>>> Jean-Michel, I will be very interested to hear your answer too.
>>>
>>> Martin, up to now, I have always assumed that PHIBOT is the bottom pressure,
>>> specifically that PHIBOT in m^2/s^2 is the anomaly relative to Depth * rhonil * g
>>> so that the absolute bottom pressure in Pa is:
>>> Depth * rhonil * g + PHIBOT * rhonil
>>>
>>> 6-hourly output of ETAN and PHIBOT is available for cube78 and cub81
>>> is available online and on skylla.
>>>
>>> Dimitris
>>>
>>> On Feb 4, 2010, at 8:00 AM, Martin Losch wrote:
>>>
>>>
>>>> Hi Jean-Michel (and Dimitris),
>>>>
>>>> we would like to cube78 or cube81 bottom pressure results. Can you confirm that for a time series of anomalies (in time) we need both PHIBOT and ETAN (for the linear free surface contribution):
>>>>
>>>> bottom pressure = gravity*rhoConst*ETAN + PHIBOT/rhoConst (+terms that are constant in time)
>>>>
>>>> From looking at calc_phi_hyd I gather, that the non-linear free surface contribution etaH is only included, when selectRstar=0, which is not the case in cube78/81, correct? But in diags_phi_rlow, it says in the head that phi_hyd_rlow is actually bottom pressure. I am confused.
>>>>
>>>> Martin
>>>>
>>>>
>>>> _______________________________________________
>>>> 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
>>>
>>
>>
>>
>> ------------------------------
>>
>> Message: 2
>> Date: Fri, 05 Feb 2010 09:43:23 -0500
>> From: Oliver Jahn <jahn at mit.edu>
>> To: mitgcm-support at mitgcm.org
>> Subject: Re: [MITgcm-support] bottom pressure diagnostics
>> Message-ID: <4B6C2E8B.3080000 at mit.edu>
>> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>>
>> Hi,
>>
>> as far as I can tell, the surface contribution is added at the end of diags_phi_rlow.F. HOWEVER, prior to cvs 2010-01-15, the vertical integral did not take the stretching of the column with nonlinFreeSurf=4 into account, i.e., the diagnostic computed
>>
>> PHIBOT = gravity*ETAN + Depth/(Depth+ETAN) * gravity * \int (rho/rhoConst) dr
>>
>> rather than
>>
>> PHIBOT = gravity*ETAN + gravity * \int (rho/rhoConst) dr
>>
>> You may want to correct for this.
>>
>> Oliver
>>
>>
>> On 2010-02-05 05:22, Jean-Michel Campin wrote:
>>
>>> Hi Dimitris,
>>>
>>> I would assume the same thing as you.
>>> I will check when I am back from vacation.
>>> Cheers,
>>> Jean-Michel
>>>
>>> On Thu, Feb 04, 2010 at 08:33:07AM -0800, Dimitris Menemenlis wrote:
>>>
>>>> Jean-Michel, I will be very interested to hear your answer too.
>>>>
>>>> Martin, up to now, I have always assumed that PHIBOT is the bottom pressure,
>>>> specifically that PHIBOT in m^2/s^2 is the anomaly relative to Depth * rhonil * g
>>>> so that the absolute bottom pressure in Pa is:
>>>> Depth * rhonil * g + PHIBOT * rhonil
>>>>
>>>> 6-hourly output of ETAN and PHIBOT is available for cube78 and cub81
>>>> is available online and on skylla.
>>>>
>>>> Dimitris
>>>>
>>>> On Feb 4, 2010, at 8:00 AM, Martin Losch wrote:
>>>>
>>>>
>>>>> Hi Jean-Michel (and Dimitris),
>>>>>
>>>>> we would like to cube78 or cube81 bottom pressure results. Can you confirm that for a time series of anomalies (in time) we need both PHIBOT and ETAN (for the linear free surface contribution):
>>>>>
>>>>> bottom pressure = gravity*rhoConst*ETAN + PHIBOT/rhoConst (+terms that are constant in time)
>>>>>
>>>>> From looking at calc_phi_hyd I gather, that the non-linear free surface contribution etaH is only included, when selectRstar=0, which is not the case in cube78/81, correct? But in diags_phi_rlow, it says in the head that phi_hyd_rlow is actually bottom pressure. I am confused.
>>>>>
>>>>> Martin
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> 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
>>>
>>
>>
>>
>> ------------------------------
>>
>> Message: 3
>> Date: Fri, 05 Feb 2010 10:03:23 -0500
>> From: Oliver Jahn <jahn at mit.edu>
>> To: mitgcm-support at mitgcm.org
>> Subject: Re: [MITgcm-support] bottom pressure diagnostics
>> Message-ID: <4B6C333B.4060409 at mit.edu>
>> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>>
>> Sorry, it is the anomaly relative to Depth*rhoConst, of course,
>>
>> cvs prior to 2010-1-15:
>>
>> PHIBOT = gravity*ETAN + Depth/(Depth+ETAN) * gravity * \int (rho/rhoConst - 1) dr
>>
>> should be:
>>
>> PHIBOT = gravity*ETAN + gravity * \int (rho/rhoConst - 1) dr
>>
>> Oliver
>>
>>
>> On 2010-02-05 09:43, Oliver Jahn wrote:
>>
>>> Hi,
>>>
>>> as far as I can tell, the surface contribution is added at the end of
>>> diags_phi_rlow.F. HOWEVER, prior to cvs 2010-01-15, the vertical
>>> integral did not take the stretching of the column with nonlinFreeSurf=4
>>> into account, i.e., the diagnostic computed
>>>
>>> PHIBOT = gravity*ETAN + Depth/(Depth+ETAN) * gravity * \int
>>> (rho/rhoConst) dr
>>>
>>> rather than
>>>
>>> PHIBOT = gravity*ETAN + gravity * \int (rho/rhoConst) dr
>>>
>>> You may want to correct for this.
>>>
>>> Oliver
>>>
>>>
>>> On 2010-02-05 05:22, Jean-Michel Campin wrote:
>>>
>>>> Hi Dimitris,
>>>>
>>>> I would assume the same thing as you.
>>>> I will check when I am back from vacation.
>>>> Cheers,
>>>> Jean-Michel
>>>>
>>>> On Thu, Feb 04, 2010 at 08:33:07AM -0800, Dimitris Menemenlis wrote:
>>>>
>>>>> Jean-Michel, I will be very interested to hear your answer too.
>>>>>
>>>>> Martin, up to now, I have always assumed that PHIBOT is the bottom
>>>>> pressure,
>>>>> specifically that PHIBOT in m^2/s^2 is the anomaly relative to Depth
>>>>> * rhonil * g
>>>>> so that the absolute bottom pressure in Pa is:
>>>>> Depth * rhonil * g + PHIBOT * rhonil
>>>>>
>>>>> 6-hourly output of ETAN and PHIBOT is available for cube78 and cub81
>>>>> is available online and on skylla.
>>>>>
>>>>> Dimitris
>>>>>
>>>>> On Feb 4, 2010, at 8:00 AM, Martin Losch wrote:
>>>>>
>>>>>
>>>>>> Hi Jean-Michel (and Dimitris),
>>>>>>
>>>>>> we would like to cube78 or cube81 bottom pressure results. Can you
>>>>>> confirm that for a time series of anomalies (in time) we need both
>>>>>> PHIBOT and ETAN (for the linear free surface contribution):
>>>>>>
>>>>>> bottom pressure = gravity*rhoConst*ETAN + PHIBOT/rhoConst (+terms
>>>>>> that are constant in time)
>>>>>>
>>>>>> From looking at calc_phi_hyd I gather, that the non-linear free
>>>>>> surface contribution etaH is only included, when selectRstar=0,
>>>>>> which is not the case in cube78/81, correct? But in diags_phi_rlow,
>>>>>> it says in the head that phi_hyd_rlow is actually bottom pressure. I
>>>>>> am confused.
>>>>>>
>>>>>> Martin
>>>>>>
>>>>>>
>>>>>> _______________________________________________
>>>>>> 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
>>>
>>
>>
>>
>> ------------------------------
>>
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>
>>
>> End of MITgcm-support Digest, Vol 80, Issue 3
>> *********************************************
>>
>>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list