[MITgcm-support] EmPmRfile, EmPmR

Dmitri Leonov dleonov at u.washington.edu
Tue Apr 3 15:58:44 EDT 2007


A recent version of CONFIG_CHECK told me that nonlinFreeSurf cannot be 
used without exactConserv.

Matthew Mazloff wrote:
> Hi Patrick,
>
> I use real fresh water flux and do not use exact conserv and 
> everything is fine.  The mean sea level does change and this is OK 
> too.  The best way to do it would be to use non-linear free surface 
> and turn off all salt flux....this is done by setting  convertFW2Salt 
> = 0.,
>
> I think you can also use  rhoConstFresh=1000. to set the density of 
> the rain water...not sure..check in external_forcing_surf.F
>
> -matt
>
>
>
>
> On Mar 22, 2007, at 11:39 AM, Patrick Rosendahl wrote:
>
>> Hi Martin,
>>
>> thanks! I use the non-hydrostatic formulation. I think you are 
>> referring to this
>> >> http://mitgcm.org/pelican/code_reference/vdb/code/167.htm#192_L
>> section of the code, where the Temp & Salt is actually added in a way 
>> that is correct for the rainy surface height increase. (Please scroll 
>> up a few lines until you see the "#ifdef EXACT_CONSERV").
>>
>> Do you think I can change the code to execute it anyway (it is not 
>> executed because of
>> >> nonlinFreeSurf.GT.0 .OR. usingPCoords
>> is false)?
>> Additionally, I probably have to change the eta-code, so that the 
>> rain is increasing the surface level. But I cannot find a good place 
>> in the code to do it.
>>
>> best,
>> Patrick
>>
>> Martin Losch wrote:
>>> Patrick,
>>> there is a flag called "useRealFreshWaterFlux", that makes the 
>>> model  add volume to the topmost grid cell. I recommend to use it 
>>> along with  exactConserv = .true. (you need to have the 
>>> corresponding CPP flag  compiled with it). That should prevent your 
>>> surface cells from  getting negative salinities.
>>> If you feel really adventurous turn on the non-linear free surface  
>>> with r*-coordinates (global_ocean.cs32x15 has an example of this),  
>>> but that's not really required for the model use actual freshwater  
>>> flux (just makes is even more "physical" and stable for the kind of  
>>> forcing that you are trying to apply).
>>> Martin
>>> PS. the EmPmR variable is used always (with or without seaice,  
>>> virtual salt flux or not).




More information about the MITgcm-support mailing list