Changeset 14789 for NEMO/branches/2021/dev_r13747_HPC-11_mcastril_HPDAonline_DiagGPU/src/ABL/sbcabl.F90
- Timestamp:
- 2021-05-05T13:18:04+02:00 (3 years ago)
- Location:
- NEMO/branches/2021/dev_r13747_HPC-11_mcastril_HPDAonline_DiagGPU
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r13747_HPC-11_mcastril_HPDAonline_DiagGPU
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev _r12970_AGRIF_CMEMSext/AGRIF5 ^/vendors/AGRIF/dev@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 ^/vendors/PPR@HEAD ext/PPR 8 9 9 10 # SETTE 10 ^/utils/CI/sette@1 3559sette11 ^/utils/CI/sette@14244 sette
-
- Property svn:externals
-
NEMO/branches/2021/dev_r13747_HPC-11_mcastril_HPDAonline_DiagGPU/src/ABL/sbcabl.F90
r13214 r14789 22 22 USE sbc_oce ! Surface boundary condition: ocean fields 23 23 USE sbcblk ! Surface boundary condition: bulk formulae 24 USE sbc blk_phy ! Surface boundary condition: bulk formulae24 USE sbc_phy ! Catalog of functions for physical/meteorological parameters in the marine boundary layer 25 25 USE dom_oce, ONLY : tmask 26 26 ! … … 36 36 USE sbc_ice, ONLY : wndm_ice, utau_ice, vtau_ice 37 37 #endif 38 #if ! defined key_ iomput38 #if ! defined key_xios 39 39 USE diawri , ONLY : dia_wri_alloc_abl 40 40 #endif … … 149 149 CALL iom_close( inum ) 150 150 151 #if ! defined key_ iomput151 #if ! defined key_xios 152 152 IF( dia_wri_alloc_abl() /= 0 ) CALL ctl_stop( 'STOP', 'abl_init : unable to allocate arrays' ) 153 153 #endif … … 175 175 IF(nn_amxl==0) WRITE(numout,*) 'Deardorff 80 length-scale ' 176 176 IF(nn_amxl==1) WRITE(numout,*) 'Modified Deardorff 80 length-scale ' 177 IF(nn_amxl==2) WRITE(numout,*) 'Bougeault and Lacarrere length-scale ' 178 IF(nn_amxl==3) WRITE(numout,*) 'Rodier et al. length-scale ' 177 IF(nn_amxl==2) WRITE(numout,*) 'Bougeault and Lacarrere length-scale ' 178 IF(nn_amxl==3) WRITE(numout,*) 'Rodier et al. length-scale ' 179 179 WRITE(numout,*) ' Minimum value of atmospheric TKE = ',tke_min,' m^2 s^-2' 180 180 WRITE(numout,*) ' Minimum value of atmospheric mixing length = ',mxl_min,' m' … … 320 320 INTEGER , INTENT(in) :: kt ! ocean time step 321 321 !! 322 REAL(wp), DIMENSION(jpi,jpj) :: zssq, zcd_du, zsen, z evp322 REAL(wp), DIMENSION(jpi,jpj) :: zssq, zcd_du, zsen, zlat, zevp 323 323 #if defined key_si3 324 324 REAL(wp), DIMENSION(jpi,jpj) :: zssqi, zcd_dui, zseni, zevpi … … 344 344 & sf(jp_uoatm)%fnow(:,:,1), sf(jp_voatm)%fnow(:,:,1), & ! <<= in 345 345 & sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1) , & ! <<= in 346 & tsk_m, zssq, zcd_du, zsen, z evp) ! =>> out346 & tsk_m, zssq, zcd_du, zsen, zlat, zevp ) ! =>> out 347 347 348 348 #if defined key_si3 … … 356 356 !! 3 - Advance ABL variables from now (n) to after (n+1) 357 357 !!------------------------------------------------------------------------------------------- 358 358 359 359 CALL abl_stp( kt, tsk_m, ssu_m, ssv_m, zssq, & ! <<= in 360 360 & sf(jp_wndi)%fnow(:,:,:), sf(jp_wndj)%fnow(:,:,:), & ! <<= in … … 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 … … 375 375 !!------------------------------------------------------------------------------------------- 376 376 377 CALL blk_oce_2( tq_abl(:,:,2,nt_n,jp_ta), & 378 & sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1), & 377 CALL blk_oce_2( tq_abl(:,:,2,nt_n,jp_ta), sf(jp_qlw )%fnow(:,:,1), & 379 378 & sf(jp_prec)%fnow(:,:,1) , sf(jp_snow)%fnow(:,:,1), & 380 & tsk_m, zsen, z evp)381 379 & tsk_m, zsen, zlat, zevp ) 380 382 381 CALL abl_rst_opn( kt ) ! Open abl restart file (if necessary) 383 382 IF( lrst_abl ) CALL abl_rst_write( kt ) ! -- abl restart file … … 385 384 #if defined key_si3 386 385 ! Avoid a USE abl in icesbc module 387 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) 388 387 #endif 389 388 END IF
Note: See TracChangeset
for help on using the changeset viewer.