ana_grid.h from upwelling test case

[angelolemos]

People,

Anyone can tell me what means the values 84.5, 66.526 (the complete equation too)and strings not declared JstrR, JendR, IstrR and IendR below from ana_grid.h?

#elif defined UPWELLING
DO j=JstrR,JendR
IF (j.le.Mm(ng)/2) THEN
val1=REAL(j,r8)
ELSE
val1=REAL(Mm(ng)+1-j,r8)
END IF
val2=MIN(depth,84.5_r8+66.526_r8*TANH((val1-10.0_r8)/7.0_r8))
DO i=IstrR,IendR
h(i,j)=val2
END DO
END DO

Thanks everybody.

[kate]

IstrR and friends are shown in a diagram here.

The shape of the bathymetry in UPWELLING is a flat interior, with hyperbolic tangent ramps up on either side. The parameters are chosen to make the tanh (almost) meet the flat interior, to have a certain width and minimum depth, etc. Note that the range of tanh is -1 to 1.