[MITgcm-devel] tracForcingOutAB
Manfredi Manizza
mmanizza at ocean.mit.edu
Thu Aug 25 17:16:50 EDT 2011
Hi
I just switched form the data.seaice from the latest version of the
code from verification in lab_sea/input_ad
and I have this error even if I set the usedynamics as false in
data.seaice
and if I comment SEAICE_deltaTdyn.
(PID.TID 0000.0001) *** ERROR *** Unsupported combination of
SEAICE_deltaTtherm,
(PID.TID 0000.0001) *** ERROR *** SEAICE_deltaTdyn, and dTtracerLev(1)
Any guess ?
Manfredi
On Thu, 25 Aug 2011, Matthew Mazloff wrote:
> Hi Manfredi,
>
> These are not the exact data.seaice I sent you (e.g. you don't want
> SEAICEdiffKhArea = 10., SEAICEdiffKhHeff = 10., if you have
> SEAICEadvScheme = 2,)
>
> But since you are not using dynamics how about jsut trying:
>
> &SEAICE_PARM01
> SEAICEwriteState = .TRUE.,
> SEAICE_initialHEFF = 0.0,
> SEAICE_deltaTtherm = 900.,
> SEAICEuseFlooding = .TRUE.,
> SEAICE_EPS = 1.E-8,
> SIsalFRAC = 0.3,
> AreaFile = 'SO6input/SIarea_ini_SOSE.bin',
> HeffFile = 'SO6input/SIheff_ini_SOSE.bin',
> SEAICE_area_floor = 1.E-5,
> SEAICE_area_reg = 0.15,
> SEAICE_hice_reg = 0.10,
> IMAX_TICE = 6,
> SEAICE_availHeatFrac = 0.8,
> SEAICEuseDYNAMICS =.FALSE.,
>
>
> Anyway, regardless - do you have:
>
> forcing_In_AB = .FALSE.,
>
> in PARM03 of data?
>
>
>
> -Matt
>
>
>
>
> On Aug 25, 2011, at 11:28 AM, Manfredi Manizza wrote:
>
>>
>> Hi,
>>
>> this is the error I have now :
>>
>> (PID.TID 0001.0001) *** ERROR *** SEAICE_CHECK: Need T,S forcing out of AB
>> (tracForcingOutAB=1)
>>
>> I have hard time to understand what means. Did I say that already ?
>> Aplogies if so.
>>
>> I tried to change tracForcingOutAB to 0 or 1 this in data
>> but the error is the same.
>>
>> At the moment I am using this set of parms from M. Mazloff in data.seaice.
>>
>> I will have a look at the points from Martin and JM.
>>
>> Thanks !
>> Manfredi
>>
>> &SEAICE_PARM01
>> #- old default for useFluxForm and DIFF1:
>> SEAICEwriteState = .TRUE.,
>> SEAICE_initialHEFF = 0.0,
>> SEAICE_deltaTtherm = 900.,
>> SEAICEuseFlooding = .TRUE.,
>> SEAICE_EPS = 1.E-8,
>> SIsalFRAC = 0.3,
>> AreaFile ='/nobackupp1/mmanizza/input/AREA.ini.SOSE.bin',
>> HeffFile ='/nobackupp1/mmanizza/input/HEFF.ini.SOSE.bin',
>> SEAICE_area_floor = 1.E-5,
>> SEAICE_area_reg = 0.15,
>> SEAICE_hice_reg = 0.10,
>> IMAX_TICE = 6,
>> SEAICE_availHeatFrac = 0.8,
>> SEAICEuseDYNAMICS =.FALSE.,
>> # SEAICE_deltaTdyn = 900.,
>> SEAICEadvScheme = 2,
>> SEAICEadvSalt = .TRUE.
>> SEAICEadvAge = .FALSE.,
>> SEAICEadvSnow = .TRUE.,
>> SEAICEdiffKhArea = 10.,
>> SEAICEdiffKhHeff = 10.,
>> DIFF1 = 4.E-3,
>>
>>
>>
>> On Thu, 25 Aug 2011, Jean-Michel Campin wrote:
>>
>>> Hi Manfredi,
>>>
>>> The MITgcm code is tested and this feature (i.e., forcing_In_AB=F)
>>> is well tested in several test-exp.
>>> Furthermore, setting of this parameter (tracForcingOutAB) is not very
>>> subtile, so it's unlikely that recent changes in the MITgcm code has
>>> something to do with this problem.
>>>
>>> For those reasons, I would be more likely to give more details about
>>> this type of issue if this was on MITgcm-support list (general type
>>> of question, general type of answer).
>>>
>>> In short, I would have answered to support list that
>>> the 1rst thing to check, from my point of view (and here I slightly
>>> disagree with Martin, because I would not add print statement until
>>> all the usual check & debug options have been tried), should be
>>> 1) clean built
>>> 2) any source code you added (in one of the "-mods=" dir)
>>>
>>> Cheers,
>>> Jean-Michel
>>>
>>> On Thu, Aug 25, 2011 at 09:26:40AM +0200, Martin Losch wrote:
>>>> Manfredi,
>>>>
>>>> I do not have any idea, the only place where tracForcingOutAB is set is
>>>> in ini_parms.F
>>>> I would check:
>>>> - what does the STDOUT say (from config_summary)
>>>> - what's the value of tracForcingOutAB is before you enter
>>>> seaice_readparms/seaice_check (needs extra print statements) to figure
>>>> out where it is reset to 0.
>>>>
>>>> and I guess you can set this variable explicitly in "data"
>>>>
>>>> Martin
>>>>
>>>> On Aug 24, 2011, at 8:58 PM, Manfredi Manizza wrote:
>>>>
>>>>>
>>>>> Martin,
>>>>>
>>>>> Nothing, I just made a mistake with using that email
>>>>> for the reply to MITgcm developers. Anyway, do you have an
>>>>> idea about the problem stated in the email ?
>>>>>
>>>>> Ciao
>>>>> Manfredi
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> On Wed, 24 Aug 2011, Martin Losch wrote:
>>>>>
>>>>>> Hi Manfredi,
>>>>>>
>>>>>> honestly, what does this question (that I cannot answer), have to do
>>>>>> with the subject "google search"?
>>>>>>
>>>>>> Martin
>>>>>>
>>>>>> On Aug 24, 2011, at 1:37 AM, Manfredi Manizza wrote:
>>>>>>
>>>>>>>
>>>>>>> Hi,
>>>>>>>
>>>>>>> I am running a job starting from obs using Matt's Southern Ocean
>>>>>>> Set-up and it crashes on this :
>>>>>>>
>>>>>>>
>>>>>>> C-- SEAICE needs forcing_In_AB FALSE
>>>>>>> IF (tracForcingOutAB.NE.1) THEN
>>>>>>> WRITE(msgBuf,'(2A)') 'SEAICE_CHECK:',
>>>>>>> & ' Need T,S forcing out of AB (tracForcingOutAB=1)'
>>>>>>> CALL PRINT_ERROR( msgBuf, myThid )
>>>>>>> STOP 'ABNORMAL END: S/R SEAICE_CHECK'
>>>>>>> ENDIF
>>>>>>>
>>>>>>> forcing_In_AB was set to FALSE my data file....
>>>>>>>
>>>>>>> Why then it's complaining ?
>>>>>>>
>>>>>>> Thanks in advance
>>>>>>> Manfredi
>>>>>>
>>>>
>>>>
>>>> _______________________________________________
>>>> MITgcm-devel mailing list
>>>> MITgcm-devel at mitgcm.org
>>>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>>>
>>> _______________________________________________
>>> MITgcm-devel mailing list
>>> MITgcm-devel at mitgcm.org
>>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>>>
>>
>> _______________________________________________
>> MITgcm-devel mailing list
>> MITgcm-devel at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
More information about the MITgcm-devel
mailing list