[MITgcm-support] small rStarFac on land....

[Andreas Klocker]

Hi Martin,
Yes I updated the code in between....I wouldn't have thought that should 
break the configuration though...
I definitely looked at  (i,j)= (26,34) - does the netcdf output 
interpret the coordinates differently?
I'll try going back to my old DIAG_SIZE even though I can't see how that 
should change anything.
A>

On 16/03/15 20:25, Martin Losch wrote:
> Hi Andreas,
>
> did you only recompile with a different DIAG_SIZE.h? Or did you update the code in between? (Not that that would explain anything, but the statement "without changing anything else I think” cannot be true (o:) Did your try going back to your original DIAG_SIZE.h?
>
> (i,j)=  (26,34) can’t be a land point either, if rStarFacC,H,eta =  0.010000  1.286500E+02 -1.273635E+02 (ie. your unperturbed water colum thickness is 128.65m and your elevation is -127.36m. I guess the maybe confusing format of the error message is my fault because I always insist on vectrizable code, and these error messages (print/write/stop statements) prohibit vectoriziation within loops, so that the code (calc_r_star.F in this case) checks for cases when a rStarFac condition is violated, counts them (in icntc1, icnts, or icntw), and then, if there are any errors (icnt* >0, l204 in calc_r_star.F), checks again (without worrying about vectorization), so that we can in the second line of output your only get the “bi,bj” count (1,1 in your case). Did you accidentally interpret that a (i,j)?
>
> M.
>
>> On 16 Mar 2015, at 01:43, Andreas Klocker <Andreas.Klocker at utas.edu.au> wrote:
>>
>> Hi everyone,
>>
>> I'm running a channel model with realistic topography at 1/20 degree resolution and have run it happily for approx. 30 years...and then I made the mistake to recompile the code to allow for more variables in diagnostics (without changing anything else I think)....and now I get:
>>
>> fail at i,j=  26  34 ; rStarFacC,H,eta =  0.010000  1.286500E+02 -1.273635E+02
>> WARNING: r*FacC < hFacInf at       1 pts : bi,bj,Thid,Iter=   1 1   1  11676434
>> STOP in CALC_R_STAR : too SMALL rStarFac[C,W,S] !
>>
>> The weird thing is that the grid point above is land, i.e. Depth = 0 and HFacC = 0, according to the grid.nc file of the last output.
>>
>> Any help would be appreciated!
>>
>> ANdreas
>>
>>
>>
>> -- 
>> ===============================================================
>> Dr. Andreas Klocker
>> ARC Fellow (DECRA)
>>
>> Institute for Marine and Antarctic Studies
>> University of Tasmania
>> 20 Castray Esplanade
>> Battery Point, TAS
>> 7004 Australia
>>
>> M:     +61 437 870 182
>> W:https://sites.google.com/site/andreasoceanographywebsite
>> skype: andiklocker
>> ===============================================================
>>
>>
>>
>> University of Tasmania Electronic Communications Policy (December, 2014).
>> This email is confidential, and is for the intended recipient only. Access, disclosure, copying, distribution, or reliance on any of it by anyone outside the intended recipient organisation is prohibited and may be a criminal offence. Please delete if obtained in error and email confirmation to the sender. The views expressed in this email are not necessarily the views of the University of Tasmania, unless clearly intended otherwise.
>>
>>
>> _______________________________________________
>> 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

-- 
===============================================================
Dr. Andreas Klocker
ARC Fellow (DECRA)

Institute for Marine and Antarctic Studies
University of Tasmania
20 Castray Esplanade
Battery Point, TAS
7004 Australia

M:     +61 437 870 182
W:     https://sites.google.com/site/andreasoceanographywebsite
skype: andiklocker
===============================================================



University of Tasmania Electronic Communications Policy (December, 2014).
This email is confidential, and is for the intended recipient only. Access, disclosure, copying, distribution, or reliance on any of it by anyone outside the intended recipient organisation is prohibited and may be a criminal offence. Please delete if obtained in error and email confirmation to the sender. The views expressed in this email are not necessarily the views of the University of Tasmania, unless clearly intended otherwise.