Ignore:
Timestamp:
2019-11-04T11:11:58+01:00 (13 months ago)
Author:
laurent
Message:

Improving syntax consistency

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_coare3p0.F90

    r11817 r11845  
    178178      LOGICAL :: l_zt_equal_zu = .FALSE.      ! if q and t are given at same height as U 
    179179      ! 
    180       REAL(wp), DIMENSION(jpi,jpj) ::  & 
    181          &  u_star, t_star, q_star, & 
    182          &  dt_zu, dq_zu,    & 
    183          &  znu_a,           & !: Nu_air, Viscosity of air 
    184          &  z0, z0t 
    185       REAL(wp), DIMENSION(jpi,jpj) ::   zeta_u        ! stability parameter at height zu 
    186       REAL(wp), DIMENSION(jpi,jpj) ::   ztmp0, ztmp1, ztmp2 
    187       REAL(wp), DIMENSION(:,:), ALLOCATABLE ::   zeta_t        ! stability parameter at height zt 
     180      REAL(wp), DIMENSION(jpi,jpj) :: u_star, t_star, q_star 
     181      REAL(wp), DIMENSION(jpi,jpj) :: dt_zu, dq_zu 
     182      REAL(wp), DIMENSION(jpi,jpj) :: znu_a,        !: Nu_air, Viscosity of air 
     183      REAL(wp), DIMENSION(jpi,jpj) :: z0, z0t 
     184      REAL(wp), DIMENSION(jpi,jpj) :: zeta_u        ! stability parameter at height zu 
     185      REAL(wp), DIMENSION(jpi,jpj) :: ztmp0, ztmp1, ztmp2 
    188186      ! 
    189187      REAL(wp), DIMENSION(:,:), ALLOCATABLE :: & 
     188         &              zeta_t,   &  ! stability parameter at height zt 
    190189         &                zsst,   &  ! to back up the initial bulk SST 
    191190         &                pdTc,   &  ! SST increment "dT" for cool-skin correction           [K] 
     
    299298         !! Roughness lengthes z0, z0t (z0q = z0t) : 
    300299         ztmp2 = u_star/vkarmn*LOG(10./z0)                                 ! Neutral wind speed at 10m 
    301          z0    = alfa_charn_3p0(ztmp2)*ztmp1/grav + 0.11_wp*znu_a/u_star   ! Roughness length (eq.6) 
     300         z0    = alfa_charn_3p0(ztmp2)*ztmp1/grav + 0.11_wp*znu_a/u_star   ! Roughness length (eq.6) [ ztmp1==u*^2 ] 
    302301         z0     = MIN( MAX(ABS(z0), 1.E-9) , 1._wp )                      ! (prevents FPE from stupid values from masked region later on) 
    303302 
Note: See TracChangeset for help on using the changeset viewer.