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

[Martin Losch]

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