Changeset 11334 for NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/src/OCE/SBC/sbcblk.F90

Timestamp:

2019-07-24T11:20:42+02:00 (5 years ago)

Author:

flemarie

Message:

Follow-on to the ABL implementation (see ticket #2131)

• Modification of diawri.F90 to output ABL variables with IOIPSL (i.e. without key_iomput)
• Modification of ice_var_agg in icevar.F90 to allow the use of tm_su in ABL and bulk
• Error handling in case ln_abl = .true. and ABL sources were not compiled
• ABL now works with SI3 considering an average over ice categories
• Update reference namelist to avoid runtime warnings due to nam_tide

• Sanity checks with ORCA2 for the following configurations :

1 - ABL src + ln_blk = .true.
2 - ABL src + ln_abl = .true.
3 - no ABL src + ln_blk = .true.

All configurations are MPP reproducible and configurations 1 and 3 results are identical

File:

• NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/src/OCE/SBC/sbcblk.F90 (modified) (4 diffs)

NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/src/OCE/SBC/sbcblk.F90

@@ -816 +816 @@
       !
       INTEGER  ::   ji, jj    ! dummy loop indices
-      REAL(wp) ::   zwndi_t , zwndj_t             ! relative wind components at T-point
+      REAL(wp) ::   zwndi_t , zwndj_t, zootm_su, zztmp1, zztmp2 ! relative wind components at T-point
       REAL(wp), DIMENSION(jpi,jpj) ::   zrhoa     ! transfer coefficient for momentum      (tau)
       REAL(wp), DIMENSION(jpi,jpj) ::   zcd_dui     ! transfer coefficient for momentum      (tau)
@@ -855 +855 @@
       zcd_dui(:,:) = wndm_ice(:,:) * Cd_atm(:,:)
       
- !     IF( ln_blk ) THEN 
+      IF( ln_blk ) THEN 
          ! ------------------------------------------------------------ !
          !    Wind stress relative to the moving ice ( U10m - U_ice )   !
@@ -874 +874 @@
          IF(ln_ctl)   CALL prt_ctl( tab2d_1=putaui  , clinfo1=' blk_ice: putaui : '   &
             &                     , tab2d_2=pvtaui  , clinfo2='          pvtaui : ' )
- !     ELSE
- !        DO jj = 1, jpj
- !           DO ji = 1, jpi
- !              pcd_dui(ji,jj) = zcd_dui(ji,jj)
- !              pseni  (ji,jj) = wndm_ice(ji,jj) * Ch_atm(ji,jj)
- !              pevpi  (ji,jj) = wndm_ice(ji,jj) * Ce_atm(ji,jj)
- !              pssqi  (ji,jj) = 11637800.0_wp * EXP( -5897.8_wp / (ptsui(ji,jj)+rt0) ) / zrhoa(ji,jj)
- !           END DO
- !        END DO
- !     ENDIF
+      ELSE
+         zztmp1 = 11637800.0_wp
+       zztmp2 =    -5897.8_wp
+       DO jj = 1, jpj
+            DO ji = 1, jpi
+               pcd_dui(ji,jj) = zcd_dui (ji,jj)
+               pseni  (ji,jj) = wndm_ice(ji,jj) * Ch_atm(ji,jj)
+               pevpi  (ji,jj) = wndm_ice(ji,jj) * Ce_atm(ji,jj)
+            zootm_su       = zztmp2 / ptsui(ji,jj)   ! ptsui is in K (it can't be zero ??)
+            pssqi  (ji,jj) = zztmp1 * EXP( zootm_su ) / zrhoa(ji,jj)
+            END DO
+         END DO
+      ENDIF
       !
       IF(ln_ctl)  CALL prt_ctl(tab2d_1=wndm_ice  , clinfo1=' blk_ice: wndm_ice : ')
@@ -1248 +1251 @@
       REAL(wp) ::   z0w, z0i, zfmi, zfmw, zfhi, zfhw
       REAL(wp) ::   zCdn_form_tmp
-      !!----------------------------------------------------------------------
-
+      !!---------------------------------------------------------------------------
+      ! fl: global variable tm_su could be used directly instead of recomputing it 
+     !
       ! mean temperature
       WHERE( at_i_b(:,:) > 1.e-20 )   ;   ztm_su(:,:) = SUM( t_su(:,:,:) * a_i_b(:,:,:) , dim=3 ) / at_i_b(:,:)