New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7235 for branches/2016/dev_r6999_CONFIGMAN_1/NEMOGCM/TOOLS/SIREN/src/grid_hgr.f90 – NEMO

Ignore:
Timestamp:
2016-11-16T11:40:55+01:00 (7 years ago)
Author:
jpaul
Message:

see ticket #1781

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6999_CONFIGMAN_1/NEMOGCM/TOOLS/SIREN/src/grid_hgr.f90

    r7233 r7235  
    1414!> 
    1515!> ** 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) 
    3831!> 
    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/> 
    4140!> 
    4241!>        This routine is given as an example, it must be modified 
     
    4645!>      second order accuracy schemes. 
    4746!> 
    48 !> N.B. If the domain is periodic, verify that scale factors are also 
     47!> @note If the domain is periodic, verify that scale factors are also 
    4948!>      periodic, and the coriolis term again. 
    5049!> 
    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 
    5856!> 
    5957!> References :   Marti, Madec and Delecluse, 1992, JGR 
Note: See TracChangeset for help on using the changeset viewer.