Changeset 11360 for NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/src/ABL/sbcabl.F90

Timestamp:

2019-07-26T20:07:25+02:00 (5 years ago)

Author:

gsamson

Message:

dev_r11265_ABL : see #2131

• use the same tm_su everywhere (Lupkes 2015 bulk and abl)
• cosmetic changes in sbcblk.F90

• change the results when using L15 because of tm_su calculation differences
• identical results when using constant sea-ice drag coefficient

File:

• NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/src/ABL/sbcabl.F90 (modified) (2 diffs)

NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/src/ABL/sbcabl.F90

@@ -335 +335 @@
   
 #if defined key_si3
-      CALL blk_ice_1(  u_abl(:,:,2,nt_n      ),  v_abl(:,:,2,nt_n      ),   &   !   <<= in
-         &            tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa),   &   !   <<= in
-         &            sf(jp_slp)%fnow(:,:,1)  ,  u_ice, v_ice           ,   &   !   <<= in
-         &            pseni=zseni, pevpi=zevpi, ptsui=tm_su             ,   &   !   <<= out
-         &            pssqi=zssqi, pcd_dui=zcd_dui                      )       !   =>> out
+      CALL blk_ice_1(  u_abl(:,:,2,nt_n      ),  v_abl(:,:,2,nt_n      ),    &   !   <<= in
+         &            tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa),    &   !   <<= in
+         &            sf(jp_slp)%fnow(:,:,1)  ,  u_ice, v_ice, tm_su    ,    &   !   <<= in
+         &            pseni=zseni, pevpi=zevpi, pssqi=zssqi, pcd_dui=zcd_dui )   !   <<= out
 #endif  
 
@@ -364 +363 @@
       !!-------------------------------------------------------------------------------------------   
 
-      CALL blk_oce_2( kt, tq_abl(:,:,2,nt_n,jp_ta),                            &
-         &                sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1),   &
-         &                sf(jp_prec)%fnow(:,:,1) , sf(jp_snow)%fnow(:,:,1),   &
-         &                sst_m, zsen, zevp                                )
+      CALL blk_oce_2( tq_abl(:,:,2,nt_n,jp_ta),                            &
+         &            sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1),   &
+         &            sf(jp_prec)%fnow(:,:,1) , sf(jp_snow)%fnow(:,:,1),   &
+         &            sst_m, zsen, zevp                                )
 
 #if defined key_si3