- Timestamp:
- 2020-05-15T19:41:01+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo
- Files:
-
- 2 edited
-
. (modified) (1 prop)
-
src/ABL/sbcabl.F90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo
- Property svn:externals
-
old new 8 8 9 9 # SETTE 10 ^/utils/CI/sette@ HEADsette10 ^/utils/CI/sette@12931 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/ABL/sbcabl.F90
r12549 r12939 68 68 LOGICAL :: lluldl 69 69 NAMELIST/namsbc_abl/ cn_dir, cn_dom, cn_ablrst_in, cn_ablrst_out, & 70 & cn_ablrst_indir, cn_ablrst_outdir, 70 & cn_ablrst_indir, cn_ablrst_outdir, ln_rstart_abl, & 71 71 & ln_hpgls_frc, ln_geos_winds, nn_dyn_restore, & 72 72 & rn_ldyn_min , rn_ldyn_max, rn_ltra_min, rn_ltra_max, & … … 263 263 264 264 ! Initialize the time index for now time (nt_n) and after time (nt_a) 265 nt_n = 1 + MOD( nit000 , 2) 266 nt_a = 1 + MOD( nit000+1, 2) 265 nt_n = 1; nt_a = 2 267 266 268 267 ! initialize ABL from data or restart 269 IF( ln_rstart ) THEN268 IF( ln_rstart_abl ) THEN 270 269 CALL abl_rst_read 271 270 ELSE … … 288 287 ENDIF 289 288 290 rhoa(:,:) = rho_air( tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), sf(jp_slp)%fnow(:,:,1) ) !!GS: rhoa must be (re)computed here here to avoid division by zero in blk_ice_1 (TBI)291 292 289 END SUBROUTINE sbc_abl_init 293 290 … … 329 326 CALL fld_read( kt, nn_fsbc, sf ) ! input fields provided at the current time-step 330 327 331 !!------------------------------------------------------------------------------------------- 332 !! 2 - Compute Cd x ||U||, Ch x ||U||, Ce x ||U||, and SSQ using now fields 333 !!------------------------------------------------------------------------------------------- 334 335 CALL blk_oce_1( kt, u_abl(:,:,2,nt_n ), v_abl(:,:,2,nt_n ), & ! <<= in 336 & tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), & ! <<= in 337 & sf(jp_slp )%fnow(:,:,1) , sst_m, ssu_m, ssv_m , & ! <<= in 338 & sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1) , & ! <<= in 339 & tsk_m, zssq, zcd_du, zsen, zevp ) ! =>> out 340 341 #if defined key_si3 342 CALL blk_ice_1( u_abl(:,:,2,nt_n ), v_abl(:,:,2,nt_n ), & ! <<= in 343 & tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), & ! <<= in 344 & sf(jp_slp)%fnow(:,:,1) , u_ice, v_ice, tm_su , & ! <<= in 345 & pseni=zseni, pevpi=zevpi, pssqi=zssqi, pcd_dui=zcd_dui ) ! <<= out 346 #endif 347 348 !!------------------------------------------------------------------------------------------- 349 !! 3 - Advance ABL variables from now (n) to after (n+1) 350 !!------------------------------------------------------------------------------------------- 351 352 CALL abl_stp( kt, tsk_m, ssu_m, ssv_m, zssq, & ! <<= in 353 & sf(jp_wndi)%fnow(:,:,:), sf(jp_wndj)%fnow(:,:,:), & ! <<= in 354 & sf(jp_tair)%fnow(:,:,:), sf(jp_humi)%fnow(:,:,:), & ! <<= in 355 & sf(jp_slp )%fnow(:,:,1), & ! <<= in 356 & sf(jp_hpgi)%fnow(:,:,:), sf(jp_hpgj)%fnow(:,:,:), & ! <<= in 357 & zcd_du, zsen, zevp, & ! <=> in/out 358 & wndm, utau, vtau, taum & ! =>> out 359 #if defined key_si3 360 & , tm_su, u_ice, v_ice, zssqi, zcd_dui & ! <<= in 361 & , zseni, zevpi, wndm_ice, ato_i & ! <<= in 362 & , utau_ice, vtau_ice & ! =>> out 363 #endif 364 & ) 365 !!------------------------------------------------------------------------------------------- 366 !! 4 - Finalize flux computation using ABL variables at (n+1), nt_n corresponds to (n+1) since 367 !! time swap is done in abl_stp 368 !!------------------------------------------------------------------------------------------- 369 370 CALL blk_oce_2( tq_abl(:,:,2,nt_n,jp_ta), & 371 & sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1), & 372 & sf(jp_prec)%fnow(:,:,1) , sf(jp_snow)%fnow(:,:,1), & 373 & tsk_m, zsen, zevp ) 374 375 CALL abl_rst_opn( kt ) ! Open abl restart file (if necessary) 376 IF( lrst_abl ) CALL abl_rst_write( kt ) ! -- abl restart file 377 378 #if defined key_si3 379 ! Avoid a USE abl in icesbc module 380 sf(jp_tair)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_ta); sf(jp_humi)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_qa) 381 #endif 328 IF( MOD( kt - 1, nn_fsbc ) == 0 ) THEN 329 330 !!------------------------------------------------------------------------------------------- 331 !! 2 - Compute Cd x ||U||, Ch x ||U||, Ce x ||U||, and SSQ using now fields 332 !!------------------------------------------------------------------------------------------- 333 334 CALL blk_oce_1( kt, u_abl(:,:,2,nt_n ), v_abl(:,:,2,nt_n ), & ! <<= in 335 & tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), & ! <<= in 336 & sf(jp_slp )%fnow(:,:,1) , sst_m, ssu_m, ssv_m , & ! <<= in 337 & sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1) , & ! <<= in 338 & tsk_m, zssq, zcd_du, zsen, zevp ) ! =>> out 339 340 #if defined key_si3 341 CALL blk_ice_1( u_abl(:,:,2,nt_n ), v_abl(:,:,2,nt_n ), & ! <<= in 342 & tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), & ! <<= in 343 & sf(jp_slp)%fnow(:,:,1) , u_ice, v_ice, tm_su , & ! <<= in 344 & pseni=zseni, pevpi=zevpi, pssqi=zssqi, pcd_dui=zcd_dui ) ! <<= out 345 #endif 346 347 !!------------------------------------------------------------------------------------------- 348 !! 3 - Advance ABL variables from now (n) to after (n+1) 349 !!------------------------------------------------------------------------------------------- 350 351 CALL abl_stp( kt, tsk_m, ssu_m, ssv_m, zssq, & ! <<= in 352 & sf(jp_wndi)%fnow(:,:,:), sf(jp_wndj)%fnow(:,:,:), & ! <<= in 353 & sf(jp_tair)%fnow(:,:,:), sf(jp_humi)%fnow(:,:,:), & ! <<= in 354 & sf(jp_slp )%fnow(:,:,1), & ! <<= in 355 & sf(jp_hpgi)%fnow(:,:,:), sf(jp_hpgj)%fnow(:,:,:), & ! <<= in 356 & zcd_du, zsen, zevp, & ! <=> in/out 357 & wndm, utau, vtau, taum & ! =>> out 358 #if defined key_si3 359 & , tm_su, u_ice, v_ice, zssqi, zcd_dui & ! <<= in 360 & , zseni, zevpi, wndm_ice, ato_i & ! <<= in 361 & , utau_ice, vtau_ice & ! =>> out 362 #endif 363 & ) 364 !!------------------------------------------------------------------------------------------- 365 !! 4 - Finalize flux computation using ABL variables at (n+1), nt_n corresponds to (n+1) since 366 !! time swap is done in abl_stp 367 !!------------------------------------------------------------------------------------------- 368 369 CALL blk_oce_2( tq_abl(:,:,2,nt_n,jp_ta), & 370 & sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1), & 371 & sf(jp_prec)%fnow(:,:,1) , sf(jp_snow)%fnow(:,:,1), & 372 & tsk_m, zsen, zevp ) 373 374 CALL abl_rst_opn( kt ) ! Open abl restart file (if necessary) 375 IF( lrst_abl ) CALL abl_rst_write( kt ) ! -- abl restart file 376 377 #if defined key_si3 378 ! Avoid a USE abl in icesbc module 379 sf(jp_tair)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_ta); sf(jp_humi)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_qa) 380 #endif 381 END IF 382 382 383 383 END SUBROUTINE sbc_abl
Note: See TracChangeset
for help on using the changeset viewer.
