Changeset 13719 for NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/src/ABL

Timestamp:

2020-11-03T12:02:23+01:00 (4 years ago)

Author:

laurent

Message:

Keep up with trunk r13718 + figure generation for STATION_ASF (ocean & ice).

Location:

NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/src/ABL

Files:

• ablmod.F90 (modified) (4 diffs)
• sbcabl.F90 (modified) (3 diffs)

NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/src/ABL/ablmod.F90

@@ -19 +19 @@
    USE sbc_oce, ONLY  : ght_abl, ghw_abl, e3t_abl, e3w_abl, jpka, jpkam1, rhoa
    USE sbcblk         ! use rn_efac, cdn_oce
-   USE sbc_phy        ! use some physical constants for flux computation
+   USE sbc_phy        ! Catalog of functions for physical/meteorological parameters in the marine boundary layer
    !
    USE prtctl         ! Print control                    (prt_ctl routine)
@@ -252 +252 @@
                   &                     + rDt_abl * fft_abl(ji, jj) * v_abl( ji, jj, jk, nt_n ) )  &
                   &                               / (1._wp + gamma_Cor*gamma_Cor*zcff)
-   
+
                v_abl( ji, jj, jk, nt_a ) =  e3t_abl(jk) *(                                         &
                   &        (1._wp-gamma_Cor*(1._wp-gamma_Cor)*zcff) * v_abl( ji, jj, jk, nt_n )    &
@@ -590 +590 @@
 
       DO_2D( 0, 1, 0, 1 )
-         zwnd_i(ji,jj) = u_abl(ji  ,jj,2,nt_a) - 0.5_wp * ( pssu(ji  ,jj) + pssu(ji-1,jj) )  
-         zwnd_j(ji,jj) = v_abl(ji,jj  ,2,nt_a) - 0.5_wp * ( pssv(ji,jj  ) + pssv(ji,jj-1) ) 
+         zwnd_i(ji,jj) = u_abl(ji  ,jj,2,nt_a) - 0.5_wp * ( pssu(ji  ,jj) + pssu(ji-1,jj) )
+         zwnd_j(ji,jj) = v_abl(ji,jj  ,2,nt_a) - 0.5_wp * ( pssv(ji,jj  ) + pssv(ji,jj-1) )
       END_2D
       !
@@ -633 +633 @@
       !    Wind stress relative to the moving ice ( U10m - U_ice )   !
       ! ------------------------------------------------------------ !
-      DO_2D( 0, 0, 0, 0 )           
+      DO_2D( 0, 0, 0, 0 )
          ptaui_ice(ji,jj) = 0.5_wp * ( rhoa(ji+1,jj) * pCd_du_ice(ji+1,jj) + rhoa(ji,jj) * pCd_du_ice(ji,jj)      )   &
             &                      * ( 0.5_wp * ( u_abl(ji+1,jj,2,nt_a) + u_abl(ji,jj,2,nt_a) ) - pssu_ice(ji,jj) )

NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/src/ABL/sbcabl.F90

@@ -22 +22 @@
    USE sbc_oce        ! Surface boundary condition: ocean fields
    USE sbcblk         ! Surface boundary condition: bulk formulae
-   USE sbc_phy        ! Surface boundary condition: bulk formulae
+   USE sbc_phy        ! Catalog of functions for physical/meteorological parameters in the marine boundary layer
    USE dom_oce, ONLY  : tmask
    !
@@ -175 +175 @@
          IF(nn_amxl==0) WRITE(numout,*) 'Deardorff 80 length-scale '
          IF(nn_amxl==1) WRITE(numout,*) 'Modified Deardorff 80 length-scale '
-         IF(nn_amxl==2) WRITE(numout,*) 'Bougeault and Lacarrere length-scale '      
-         IF(nn_amxl==3) WRITE(numout,*) 'Rodier et al. length-scale '   
+         IF(nn_amxl==2) WRITE(numout,*) 'Bougeault and Lacarrere length-scale '
+         IF(nn_amxl==3) WRITE(numout,*) 'Rodier et al. length-scale '
          WRITE(numout,*) ' Minimum value of atmospheric TKE           = ',tke_min,' m^2 s^-2'
          WRITE(numout,*) ' Minimum value of atmospheric mixing length = ',mxl_min,' m'
@@ -356 +356 @@
          !! 3 - Advance ABL variables from now (n) to after (n+1)
          !!-------------------------------------------------------------------------------------------
-   
+
          CALL abl_stp( kt, tsk_m, ssu_m, ssv_m, zssq,                          &   !   <<= in
             &              sf(jp_wndi)%fnow(:,:,:), sf(jp_wndj)%fnow(:,:,:),   &   !   <<= in