Changeset 14235 for NEMO/trunk/src/ABL/ablmod.F90

Timestamp:

2020-12-21T14:22:03+01:00 (4 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)

File:

• NEMO/trunk/src/ABL/ablmod.F90 (modified) (4 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