[MITgcm-support] calc_r_star

[Elja Huibregtse]

Hi all,

I am using the example of the cube-sphere integration with sea-ice on the cs32 grid (mitgcm_contrib/high_res_cube/README.cs32). 

Now, I would like to add a tidal potential. Therefore, I added a couple of lines in 'external_forcing.f' to calculate the gradient of this tidal potential field, both in zonal and meridional direction, covering the entire body. However, I immediately get the next error message:
ABNORMAL END: S/R CALC_R_STAR

The message-file gives:
WARNING: r*FacW < hFacInf at     499 pts : bi,bj,Thid,Iter=   1   1   1         1
 e.g. at i,j=   8  32 ; rStarFac,H,eta = -0.318636  9.750000E+04 -1.805742E+05 -8.111568E+04
STOP in CALC_R_STAR : too SMALL rStarFacW !
WARNING: r*FacS < hFacInf at     474 pts : bi,bj,Thid,Iter=   1   1   1         1
 e.g. at i,j=  11  32 ; rStarFac,H,eta = -0.049461  9.750000E+04 -4.993763E+05  2.969881E+05
STOP in CALC_R_STAR : too SMALL rStarFacS !
WARNING: r*FacC < hFacInf at     469 pts : bi,bj,Thid,Iter=   1   1   1         1
 e.g. at i,j=   7  32 ; rStarFac,H,eta = -0.852043  9.750000E+04 -1.805742E+05
STOP in CALC_R_STAR : too SMALL rStarFacC !

Do you maybe have some suggestions what is going on and how to fix it?
Thanks a lot!

Elja



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20080427/ac0afbbe/attachment.htm>