- Timestamp:
- 2016-11-16T11:40:55+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6999_CONFIGMAN_1/NEMOGCM/TOOLS/SIREN/src/grid_hgr.f90
r7233 r7235 14 14 !> 15 15 !> ** Method : The geographical position of the model grid-points is 16 !> defined from analytical functions, fslam and fsphi, the deriva- 17 !> tives of which gives the horizontal scale factors e1,e2. 18 !> Defining two function fslam and fsphi and their derivatives in 19 !> the two horizontal directions (fse1 and fse2), the model grid- 20 !> point position and scale factors are given by: 21 !> t-point:<br/> 22 !> glamt(i,j) = fslam(i ,j ) e1t(i,j) = fse1(i ,j )<br/> 23 !> gphit(i,j) = fsphi(i ,j ) e2t(i,j) = fse2(i ,j )<br/> 24 !> u-point:<br/> 25 !> glamu(i,j) = fslam(i+1/2,j ) e1u(i,j) = fse1(i+1/2,j )<br/> 26 !> gphiu(i,j) = fsphi(i+1/2,j ) e2u(i,j) = fse2(i+1/2,j )<br/> 27 !> v-point:<br/> 28 !> glamv(i,j) = fslam(i ,j+1/2) e1v(i,j) = fse1(i ,j+1/2)<br/> 29 !> gphiv(i,j) = fsphi(i ,j+1/2) e2v(i,j) = fse2(i ,j+1/2)<br/> 30 !> f-point:<br/> 31 !> glamf(i,j) = fslam(i+1/2,j+1/2) e1f(i,j) = fse1(i+1/2,j+1/2)<br/> 32 !> gphif(i,j) = fsphi(i+1/2,j+1/2) e2f(i,j) = fse2(i+1/2,j+1/2)<br/> 33 !> Where fse1 and fse2 are defined by:<br/> 34 !> fse1(i,j) = ra * rad * SQRT( (cos(phi) di(fslam))**2 35 !> + di(fsphi) **2 )(i,j)<br/> 36 !> fse2(i,j) = ra * rad * SQRT( (cos(phi) dj(fslam))**2 37 !> + dj(fsphi) **2 )(i,j)<br/> 16 !> defined from analytical functions, fslam and fsphi, the derivatives of which gives the horizontal scale factors e1,e2. 17 !> Defining two function fslam and fsphi and their derivatives in the two horizontal directions (fse1 and fse2), 18 !> the model grid-point position and scale factors are given by: 19 !> - t-point: 20 !> - glamt(i,j) = fslam(i ,j ) e1t(i,j) = fse1(i ,j ) 21 !> - gphit(i,j) = fsphi(i ,j ) e2t(i,j) = fse2(i ,j ) 22 !> - u-point: 23 !> - glamu(i,j) = fslam(i+1/2,j ) e1u(i,j) = fse1(i+1/2,j ) 24 !> - gphiu(i,j) = fsphi(i+1/2,j ) e2u(i,j) = fse2(i+1/2,j ) 25 !> - v-point: 26 !> - glamv(i,j) = fslam(i ,j+1/2) e1v(i,j) = fse1(i ,j+1/2) 27 !> - gphiv(i,j) = fsphi(i ,j+1/2) e2v(i,j) = fse2(i ,j+1/2) 28 !> - f-point: 29 !> - glamf(i,j) = fslam(i+1/2,j+1/2) e1f(i,j) = fse1(i+1/2,j+1/2) 30 !> - gphif(i,j) = fsphi(i+1/2,j+1/2) e2f(i,j) = fse2(i+1/2,j+1/2) 38 31 !> 39 !> The coriolis factor is given at z-point by:<br/> 40 !> ff = 2.*omega*sin(gphif) (in s-1)<br/> 32 !> Where fse1 and fse2 are defined by: 33 !> - fse1(i,j) = ra * rad * SQRT( (cos(phi) di(fslam))**2 34 !> + di(fsphi) **2 )(i,j) 35 !> - fse2(i,j) = ra * rad * SQRT( (cos(phi) dj(fslam))**2 36 !> + dj(fsphi) **2 )(i,j) 37 !> 38 !> The coriolis factor is given at z-point by: 39 !> - ff = 2.*omega*sin(gphif) (in s-1)<br/> 41 40 !> 42 41 !> This routine is given as an example, it must be modified … … 46 45 !> second order accuracy schemes. 47 46 !> 48 !> N.B.If the domain is periodic, verify that scale factors are also47 !> @note If the domain is periodic, verify that scale factors are also 49 48 !> periodic, and the coriolis term again. 50 49 !> 51 !> ** Action : - define glamt, glamu, glamv, glamf: longitude of t-, 52 !> u-, v- and f-points (in degre) 53 !> - define gphit, gphiu, gphiv, gphit: latitude of t-, 54 !> u-, v- and f-points (in degre) 55 !> define e1t, e2t, e1u, e2u, e1v, e2v, e1f, e2f: horizontal 56 !> scale factors (in meters) at t-, u-, v-, and f-points. 57 !> define ff: coriolis factor at f-point 50 !> ** Action : 51 !> - define glamt, glamu, glamv, glamf: longitude of t-, u-, v- and f-points (in degre) 52 !> - define gphit, gphiu, gphiv, gphit: latitude of t-, u-, v- and f-points (in degre) 53 !> - define e1t, e2t, e1u, e2u, e1v, e2v, e1f, e2f: horizontal 54 !> - scale factors (in meters) at t-, u-, v-, and f-points. 55 !> - define ff: coriolis factor at f-point 58 56 !> 59 57 !> References : Marti, Madec and Delecluse, 1992, JGR
Note: See TracChangeset
for help on using the changeset viewer.