Changeset 9154


Ignore:
Timestamp:
2017-12-21T14:25:59+01:00 (3 years ago)
Author:
clem
Message:

SAS_BIPER test case: change the position of the grid to match Jerome's changes in VORTEX test case

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  
    6363      ! 
    6464      INTEGER  ::   ji, jj   ! dummy loop indices 
    65       REAL(wp) ::   zphi0, zlam0 
     65      REAL(wp) ::   zphi0, zlam0, zbeta, zf0 
    6666      REAL(wp) ::   zti, zui, zvi, zfi   ! local scalars 
    6767      REAL(wp) ::   ztj, zuj, zvj, zfj   !   -      - 
     
    7575 
    7676      !                          ========== 
    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 
    8082#if defined key_agrif  
    8183      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 
    8490         rn_dx = Agrif_Parent(rn_dx)/Agrif_Rhox() 
    8591         rn_dy = Agrif_Parent(rn_dy)/Agrif_Rhoy()           
     
    146152      kff = 1                       !  indicate not to compute Coriolis parameter afterward 
    147153      ! 
     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 
    148158      pff_f(:,:) = 0._wp            ! here No earth rotation: f=0 
    149159      pff_t(:,:) = 0._wp 
Note: See TracChangeset for help on using the changeset viewer.