Changeset 14235 for NEMO/trunk/src/ABL

Timestamp:

2020-12-21T14:22:03+01:00 (3 years ago)

Author:

gsamson

Message:

make ABL work again with new bulks (output LHF from ablmod and rename air theta+humidity in sbcabl for ice bulks)

Location:

NEMO/trunk/src/ABL

Files:

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

NEMO/trunk/src/ABL/ablmod.F90

@@ -41 +41 @@
               &            pu_dta, pv_dta, pt_dta, pq_dta,    &
               &            pslp_dta, pgu_dta, pgv_dta,        &
-              &            pcd_du, psen, pevp,                &     ! in/out
+              &            pcd_du, psen, pevp, plat,          &     ! in/out
               &            pwndm, ptaui, ptauj, ptaum         &
 #if defined key_si3
@@ -84 +84 @@
       REAL(wp) , INTENT(inout), DIMENSION(:,:  ) ::   pevp       ! Ce x Du
       REAL(wp) , INTENT(inout), DIMENSION(:,:  ) ::   pwndm      ! ||uwnd||
+      REAL(wp) , INTENT(  out), DIMENSION(:,:  ) ::   plat       ! latent heat flux
       REAL(wp) , INTENT(  out), DIMENSION(:,:  ) ::   ptaui      ! taux
       REAL(wp) , INTENT(  out), DIMENSION(:,:  ) ::   ptauj      ! tauy
@@ -530 +531 @@
       !-------------
       !                            !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
-      !                            !  6 *** MPI exchanges
+      !                            !  6 *** MPI exchanges & IOM outputs
       !                            !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
       !
@@ -590 +591 @@
          psen( ji, jj ) =    - cp_air(zhumi) * zcff * psen(ji,jj) * ( psst(ji,jj) + rt0 - ztemp )   !GS: negative sign to respect aerobulk convention
          pevp( ji, jj ) = rn_efac*MAX( 0._wp,  zcff * pevp(ji,jj) * ( pssq(ji,jj)       - zhumi ) )
+         plat( ji, jj ) = - L_vap( psst(ji,jj) ) * pevp( ji, jj )
          rhoa( ji, jj ) = zcff
       END_2D

NEMO/trunk/src/ABL/sbcabl.F90

@@ -363 +363 @@
             &              sf(jp_hpgi)%fnow(:,:,:), sf(jp_hpgj)%fnow(:,:,:),   &   !   <<= in
             &              zcd_du, zsen, zevp,                                 &   !   <=> in/out
-            &              wndm, utau, vtau, taum                              &   !   =>> out
+            &              zlat, wndm, utau, vtau, taum                        &   !   =>> out
 #if defined key_si3
             &            , tm_su, u_ice, v_ice, zssqi, zcd_dui                 &   !   <<= in
@@ -377 +377 @@
          CALL blk_oce_2( tq_abl(:,:,2,nt_n,jp_ta), sf(jp_qlw )%fnow(:,:,1),   &
             &            sf(jp_prec)%fnow(:,:,1) , sf(jp_snow)%fnow(:,:,1),   &
-            &            tsk_m, zsen, zlat, zevp                                )
+            &            tsk_m, zsen, zlat, zevp                              )
 
          CALL abl_rst_opn( kt )                       ! Open abl restart file (if necessary)
@@ -384 +384 @@
 #if defined key_si3
          ! Avoid a USE abl in icesbc module
-         sf(jp_tair)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_ta);  sf(jp_humi)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_qa)
+         theta_air_zt = tq_abl(:,:,2,nt_n,jp_ta) ;   q_air_zt = tq_abl(:,:,2,nt_n,jp_qa)
 #endif
       END IF