- Timestamp:
- 2017-12-21T14:25:59+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/MY_SRC/usrdef_hgr.F90
r9124 r9154 63 63 ! 64 64 INTEGER :: ji, jj ! dummy loop indices 65 REAL(wp) :: zphi0, zlam0 65 REAL(wp) :: zphi0, zlam0, zbeta, zf0 66 66 REAL(wp) :: zti, zui, zvi, zfi ! local scalars 67 67 REAL(wp) :: ztj, zuj, zvj, zfj ! - - … … 75 75 76 76 ! ========== 77 zlam0 = 0._wp 78 zphi0 = 0._wp 79 77 !clem zlam0 = 0._wp 78 !clem zphi0 = 0._wp 79 zlam0 = -(jpiglo-1)/2 * 1.e-3 * rn_dx 80 zphi0 = -(jpjglo-1)/2 * 1.e-3 * rn_dy 81 80 82 #if defined key_agrif 81 83 IF( .NOT. Agrif_Root() ) THEN 82 zlam0 = Agrif_Parent(zlam0) + (Agrif_ix())*Agrif_Parent(rn_dx) * 1.e-5 83 zphi0 = Agrif_Parent(zphi0) + (Agrif_iy())*Agrif_Parent(rn_dy) * 1.e-5 84 !clem zlam0 = Agrif_Parent(zlam0) + (Agrif_ix())*Agrif_Parent(rn_dx) * 1.e-5 85 !clem zphi0 = Agrif_Parent(zphi0) + (Agrif_iy())*Agrif_Parent(rn_dy) * 1.e-5 86 zlam0 = ( 0.5_wp - ( Agrif_parent(jpiglo) - 1 ) / 2 ) * 1.e-3 * Agrif_irhox() * rn_dx & 87 & + ( Agrif_Ix() + nbghostcells - 1 ) * Agrif_irhox() * rn_dx * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dx * 1.e-3 88 zphi0 = ( 0.5_wp - ( Agrif_parent(jpjglo) - 1 ) / 2 ) * 1.e-3 * Agrif_irhoy() * rn_dy & 89 & + ( Agrif_Iy() + nbghostcells - 1 ) * Agrif_irhoy() * rn_dy * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dy * 1.e-3 84 90 rn_dx = Agrif_Parent(rn_dx)/Agrif_Rhox() 85 91 rn_dy = Agrif_Parent(rn_dy)/Agrif_Rhoy() … … 146 152 kff = 1 ! indicate not to compute Coriolis parameter afterward 147 153 ! 154 !!jerome zbeta = 2._wp * omega * COS( rad * rn_ppgphi0 ) / ra 155 !! zf0 = 2._wp * omega * SIN( rad * rn_ppgphi0 ) 156 !! pff_f(:,:) = zf0 + zbeta * pphif(:,:) * 1.e+3 157 !!jerome pff_t(:,:) = zf0 + zbeta * pphit(:,:) * 1.e+3 148 158 pff_f(:,:) = 0._wp ! here No earth rotation: f=0 149 159 pff_t(:,:) = 0._wp
Note: See TracChangeset
for help on using the changeset viewer.