- Timestamp:
- 2020-12-21T14:22:03+01:00 (4 years ago)
- Location:
- NEMO/trunk/src/ABL
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ABL/ablmod.F90
r14215 r14235 41 41 & pu_dta, pv_dta, pt_dta, pq_dta, & 42 42 & pslp_dta, pgu_dta, pgv_dta, & 43 & pcd_du, psen, pevp, 43 & pcd_du, psen, pevp, plat, & ! in/out 44 44 & pwndm, ptaui, ptauj, ptaum & 45 45 #if defined key_si3 … … 84 84 REAL(wp) , INTENT(inout), DIMENSION(:,: ) :: pevp ! Ce x Du 85 85 REAL(wp) , INTENT(inout), DIMENSION(:,: ) :: pwndm ! ||uwnd|| 86 REAL(wp) , INTENT( out), DIMENSION(:,: ) :: plat ! latent heat flux 86 87 REAL(wp) , INTENT( out), DIMENSION(:,: ) :: ptaui ! taux 87 88 REAL(wp) , INTENT( out), DIMENSION(:,: ) :: ptauj ! tauy … … 530 531 !------------- 531 532 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 532 ! ! 6 *** MPI exchanges 533 ! ! 6 *** MPI exchanges & IOM outputs 533 534 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 534 535 ! … … 590 591 psen( ji, jj ) = - cp_air(zhumi) * zcff * psen(ji,jj) * ( psst(ji,jj) + rt0 - ztemp ) !GS: negative sign to respect aerobulk convention 591 592 pevp( ji, jj ) = rn_efac*MAX( 0._wp, zcff * pevp(ji,jj) * ( pssq(ji,jj) - zhumi ) ) 593 plat( ji, jj ) = - L_vap( psst(ji,jj) ) * pevp( ji, jj ) 592 594 rhoa( ji, jj ) = zcff 593 595 END_2D -
NEMO/trunk/src/ABL/sbcabl.F90
r14072 r14235 363 363 & sf(jp_hpgi)%fnow(:,:,:), sf(jp_hpgj)%fnow(:,:,:), & ! <<= in 364 364 & zcd_du, zsen, zevp, & ! <=> in/out 365 & wndm, utau, vtau, taum& ! =>> out365 & zlat, wndm, utau, vtau, taum & ! =>> out 366 366 #if defined key_si3 367 367 & , 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), & 378 378 & sf(jp_prec)%fnow(:,:,1) , sf(jp_snow)%fnow(:,:,1), & 379 & tsk_m, zsen, zlat, zevp 379 & tsk_m, zsen, zlat, zevp ) 380 380 381 381 CALL abl_rst_opn( kt ) ! Open abl restart file (if necessary) … … 384 384 #if defined key_si3 385 385 ! Avoid a USE abl in icesbc module 386 sf(jp_tair)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_ta); sf(jp_humi)%fnow(:,:,1)= tq_abl(:,:,2,nt_n,jp_qa)386 theta_air_zt = tq_abl(:,:,2,nt_n,jp_ta) ; q_air_zt = tq_abl(:,:,2,nt_n,jp_qa) 387 387 #endif 388 388 END IF
Note: See TracChangeset
for help on using the changeset viewer.