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 12062 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/MY_SRC/isfcavgam.F90 – NEMO

Ignore:
Timestamp:
2019-12-05T10:57:46+01:00 (4 years ago)
Author:
mathiot
Message:

changes required by N.J. review

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/MY_SRC/isfcavgam.F90

    r11987 r12062  
    6565         ! gamma is constant (specified in namelist) 
    6666         ! nothing to do 
    67       CASE ('ad15', 'hj99') 
     67      CASE ('vel', 'vel_stab') 
    6868         ! compute velocity in tbl 
    6969         CALL isf_tbl(un(:,:,:) ,zutbl(:,:),'U', miku, rhisf_tbl_cav) 
     
    8989         pgt(:,:) = rn_gammat0 
    9090         pgs(:,:) = rn_gammas0 
    91       CASE ( 'ad15' ) ! gamma is proportional to u* 
    92          CALL gammats_AD15 (              zutbl, zvtbl, rCd0_top, rn_vtide**2,               pgt, pgs ) 
    93       CASE ( 'hj99' ) ! gamma depends of stability of boundary layer and u* 
    94          CALL gammats_HJ99 (pttbl, pstbl, zutbl, zvtbl, rCd0_top, r_ke0_top, pqoce, pqfwf, pgt, pgs ) 
     91      CASE ( 'vel' ) ! gamma is proportional to u* 
     92         CALL gammats_vel      (              zutbl, zvtbl, rCd0_top, rn_vtide**2,             pgt, pgs ) 
     93      CASE ( 'vel_stab' ) ! gamma depends of stability of boundary layer and u* 
     94         CALL gammats_vel_stab (pttbl, pstbl, zutbl, zvtbl, rCd0_top, r_ke0_top, pqoce, pqfwf, pgt, pgs ) 
    9595      CASE DEFAULT 
    9696         CALL ctl_stop('STOP','method to compute gamma (cn_gammablk) is unknown (should not see this)') 
     
    115115   !!----------------------------------------------------------------------------------------------------- 
    116116   ! 
    117    SUBROUTINE gammats_AD15(putbl, pvtbl, pCd, pke2, &   ! <<== in 
     117   SUBROUTINE gammats_vel( putbl, pvtbl, pCd, pke2, &   ! <<== in 
    118118      &                                  pgt, pgs   )   ! ==>> out gammats [m/s] 
    119119      !!---------------------------------------------------------------------- 
     
    145145      CALL iom_put('isfustar',zustar(:,:)) 
    146146      ! 
    147    END SUBROUTINE gammats_AD15 
    148  
    149    SUBROUTINE gammats_HJ99( pttbl, pstbl, putbl, pvtbl, pCd, pke2, pqoce, pqfwf, &  ! <<== in 
    150       &                                                            pgt  , pgs    )  ! ==>> out gammats [m/s] 
     147   END SUBROUTINE gammats_vel 
     148 
     149   SUBROUTINE gammats_vel_stab( pttbl, pstbl, putbl, pvtbl, pCd, pke2, pqoce, pqfwf, &  ! <<== in 
     150      &                                                                pgt  , pgs    )  ! ==>> out gammats [m/s] 
    151151      !!---------------------------------------------------------------------- 
    152152      !! ** Purpose    : compute the coefficient echange coefficient  
     
    252252      END DO 
    253253 
    254    END SUBROUTINE gammats_HJ99 
     254   END SUBROUTINE gammats_vel_stab 
    255255 
    256256END MODULE isfcavgam 
Note: See TracChangeset for help on using the changeset viewer.