[MITgcm-support] Bulk Force and SSS restoring

gael forget gforget at mit.edu
Tue Feb 23 20:07:09 EST 2010


Hi Amato, Martin, et al.

I concur with Holly's statement.

I would add that it may be good practice to put all of the forcing
specifications in data.exf if you use exf -- or in data if you do not
use exf -- but not half and half. In fact, I made the recent MITgcm
modification for SST/SSS relaxation because I had tripped several
times over their specification by means of data+data.exf.

Cheers,
Gael

On Feb 22, 2010, at 10:22 AM, Holly Dail wrote:

> Martin -
>
> In my experience, when using EXF it is not necessary to set  
> externForcingPeriod or externForcingCycle in the data file.   The  
> data.exf repeatPeriod and field specific periods (e.g. atempperiod)  
> seem sufficient:
>
>  &EXF_NML_01
>  exf_albedo        = 0.15,
>  exf_iprec         = 32,
>  exf_yftype        = 'RL',
>  useExfYearlyFields= .FALSE.,
>  twoDigitYear      = .FALSE.,
>  useExfCheckRange  = .FALSE.,
>  repeatPeriod      = 31104000.,
>  &
> #
>  &EXF_NML_02
>  atempfile         = 'RG_tmp2m06',
> <snip>
> #
>  atempstartdate1   = 20010101,
>  atempstartdate2   = 060000,
>  atempperiod       = 21600.0,
> #
> <snip>
>
> Holly
>
> On Feb 18, 2010, at Feb 18 , 1:57 AM, Martin Losch wrote:
>
>> Hi Amato,
>> On Feb 17, 2010, at 8:57 PM, Amato Evan wrote:
>>
>>>>
>>>> In exf, you can specify wind (or stress), lwdown, swdown,  
>>>> humidity, surface air temperature, precipitation (and with an  
>>>> additional flag, also evaporation), or net lwflux and net swflux  
>>>> (depends on the CPP-flags, see EXF_OPTIONS.h for more details).
>>> This is what I am doing (u & v stress, lw down, sw down, air q &  
>>> temp). I'd like to leave off precip and just restore to SSS since  
>>> this is how I spun up the model to begin with (using "external  
>>> fields load").
>> when you define EXF_READ_EVAP in EXF_OPTIONS.h, you can read both  
>> precipitation and evaporation with exf. If you then do not specify  
>> input files both fields (precip and evap) will be zero, so no FW  
>> flux. with
>> #define ALLOW_CLIMSSS_RELAXATION
>> you can read SSS for restoring (the restoring time scale set in  
>> data, and not data.exf unfortunately. However, if I interpret the  
>> most recent check-ins by Gael correctly, then this has just been  
>> changed last night and you can now, after a cvs-update set  
>> climsssTauRelax in data.exf, instead of tauSaltClimRelax in data.  
>> Thank you, Gael, that was long overdue.
>> What I don't know: Do I still need to set
>> externForcingPeriod=2592000.,
>> externForcingCycle=31104000.,
>> or to some non-zero value in data to get periodic forcing, even  
>> with exf, or has this been fixed, too?)
>>
>> M
>>
>> _______________________________________________
>> 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




More information about the MITgcm-support mailing list