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

Timestamp:

2019-07-22T17:00:21+02:00 (5 years ago)

Author:

flemarie

Message:

First implementation of ABL (see ticket #2131)

• Update reference and cfg namelists for ORCA2_ICE_ABL
• Run ABL over the ocean and BLK over sea-ice (ABL over sea-ice to come)
• Bug fix in computation of pblh (+ add option to smooth pblh)

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

@@ -477 +477 @@
          DO jj = 1, jpj
             DO ji = 1, jpi
-               zztmp = zU_zu(ji,jj) * tmask(ji,jj,1)
+               zztmp = zU_zu(ji,jj) !* tmask(ji,jj,1)
                wndm(ji,jj)   = zztmp                   ! Store zU_zu in wndm to compute ustar2 in ablmod 
                pcd_du(ji,jj) = zztmp * Cd_atm(ji,jj)
@@ -639 +639 @@
          CALL prt_ctl(tab2d_1=zqla , clinfo1=' blk_oce_2: zqla  : ', tab2d_2=qsr  , clinfo2=' qsr   : ')
          CALL prt_ctl(tab2d_1=emp  , clinfo1=' blk_oce_2: emp   : ')
-      ENDIF
-      !
+      ENDIF 
+     !
    END SUBROUTINE blk_oce_2
 
@@ -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
+ !        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
       !
       IF(ln_ctl)  CALL prt_ctl(tab2d_1=wndm_ice  , clinfo1=' blk_ice: wndm_ice : ')