Changeset 6913
- Timestamp:
- 2016-09-06T13:00:48+02:00 (8 years ago)
- Location:
- branches/UKMO/dev_r5518_haney_arctic_mask_for_tests/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 2 added
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_haney_arctic_mask_for_tests/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90
r5407 r6913 44 44 LOGICAL , PUBLIC :: ln_dm2dc !: Daily mean to Diurnal Cycle short wave (qsr) 45 45 LOGICAL , PUBLIC :: ln_rnf !: runoffs / runoff mouths 46 LOGICAL , PUBLIC :: ln_rnf_arcmsk !: runoffs / runoff mouths + Arctic mask 46 47 LOGICAL , PUBLIC :: ln_ssr !: Sea Surface restoring on SST and/or SSS 48 LOGICAL , PUBLIC :: ln_flx_adj !: Sea Surface restoring on SST and/or SSS 47 49 LOGICAL , PUBLIC :: ln_apr_dyn !: Atmospheric pressure forcing used on dynamics (ocean & ice) 48 50 INTEGER , PUBLIC :: nn_ice !: flag for ice in the surface boundary condition (=0/1/2/3) -
branches/UKMO/dev_r5518_haney_arctic_mask_for_tests/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90
r5501 r6913 41 41 USE cpl_oasis3 ! OASIS routines for coupling 42 42 USE sbcssr ! surface boundary condition: sea surface restoring 43 USE sbcflx_adj ! surface boundary condition: sea surface flux adjustment 43 44 USE sbcrnf ! surface boundary condition: runoffs 45 USE sbc_arcmsk ! surface boundary condition: runoffs & Arctic mask 44 46 USE sbcisf ! surface boundary condition: ice shelf 45 47 USE sbcfwb ! surface boundary condition: freshwater budget … … 89 91 & ln_blk_mfs, ln_apr_dyn, nn_ice, nn_ice_embd, ln_dm2dc , ln_rnf , & 90 92 & ln_ssr , nn_isf , nn_fwb, ln_cdgw , ln_wave , ln_sdw , & 91 & nn_lsm , nn_limflx , nn_components, ln_cpl 93 & nn_lsm , nn_limflx , nn_components, ln_cpl, ln_flx_adj,ln_rnf_arcmsk 92 94 INTEGER :: ios 93 95 INTEGER :: ierr, ierr0, ierr1, ierr2, ierr3, jpm … … 141 143 WRITE(numout,*) ' daily mean to diurnal cycle qsr ln_dm2dc = ', ln_dm2dc 142 144 WRITE(numout,*) ' runoff / runoff mouths ln_rnf = ', ln_rnf 145 WRITE(numout,*) ' runoff / runoff mouths & Arctic mask ln_rnf_arcmsk= ', ln_rnf_arcmsk 143 146 WRITE(numout,*) ' iceshelf formulation nn_isf = ', nn_isf 144 147 WRITE(numout,*) ' Sea Surface Restoring on SST and/or SSS ln_ssr = ', ln_ssr 148 WRITE(numout,*) ' Sea Surface Flux adjustment on heat and/or freshwater ln_flx_adj = ', ln_flx_adj 145 149 WRITE(numout,*) ' FreshWater Budget control (=0/1/2) nn_fwb = ', nn_fwb 146 150 WRITE(numout,*) ' closed sea (=0/1) (set in namdom) nn_closea = ', nn_closea … … 297 301 CALL sbc_ssm_init ! Sea-surface mean fields initialisation 298 302 ! 303 ! KY 12/11/2015 304 IF( ln_flx_adj ) CALL sbc_flx_adj_init ! Sea-Surface Flux adjustment 305 ! 299 306 IF( ln_ssr ) CALL sbc_ssr_init ! Sea-Surface Restoring initialisation 300 307 ! 301 308 CALL sbc_rnf_init ! Runof initialisation 309 ! 310 ! KY 16/03/2016 311 IF( ln_rnf_arcmsk ) CALL sbc_rnf_init_arcmsk ! Runoff & Arctic mask initialisation 302 312 ! 303 313 IF( nn_ice == 3 ) CALL sbc_lim_init ! LIM3 initialisation … … 399 409 IF( ln_rnf ) CALL sbc_rnf( kt ) ! add runoffs to fresh water fluxes 400 410 411 ! KY 12/11/2015 412 IF( ln_flx_adj ) CALL sbc_flx_adj( kt ) ! add flux adjustment term 413 401 414 IF( ln_ssr ) CALL sbc_ssr( kt ) ! add SST/SSS damping term 402 415 -
branches/UKMO/dev_r5518_haney_arctic_mask_for_tests/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssr.F90
r4990 r6913 17 17 USE phycst ! physical constants 18 18 USE sbcrnf ! surface boundary condition : runoffs 19 USE sbc_arcmsk ! surface boundary condition : runoffs & Arctic mask 19 20 ! 20 21 USE fldread ! read input fields … … 97 98 DO jj = 1, jpj 98 99 DO ji = 1, jpi 99 zqrp = rn_dqdt * ( sst_m(ji,jj) - sf_sst(1)%fnow(ji,jj,1) ) 100 ! K.Y. 16/03/2016 Apply Arctic mask to SST restoring 101 zqrp = rn_dqdt * ( 1. - 2.*only_arcmsk(ji,jj) ) & 102 & * ( sst_m(ji,jj) - sf_sst(1)%fnow(ji,jj,1) ) 103 !zqrp = rn_dqdt * ( sst_m(ji,jj) - sf_sst(1)%fnow(ji,jj,1) ) 100 104 qns(ji,jj) = qns(ji,jj) + zqrp 101 105 qrp(ji,jj) = zqrp … … 103 107 END DO 104 108 CALL iom_put( "qrp", qrp ) ! heat flux damping 109 CALL iom_put( "only_arcmsk", only_arcmsk ) ! arctic mask 105 110 ENDIF 106 111 ! … … 110 115 DO jj = 1, jpj 111 116 DO ji = 1, jpi 112 zerp = zsrp * ( 1. - 2.*rnfmsk (ji,jj) ) & ! No damping in vicinity of river mouths117 zerp = zsrp * ( 1. - 2.*rnfmsk_arcmsk(ji,jj) ) & ! No damping in vicinity of river mouths 113 118 & * ( sss_m(ji,jj) - sf_sss(1)%fnow(ji,jj,1) ) 114 119 sfx(ji,jj) = sfx(ji,jj) + zerp ! salt flux 115 120 erp(ji,jj) = zerp / MAX( sss_m(ji,jj), 1.e-20 ) ! converted into an equivalent volume flux (diagnostic only) 116 121 END DO 122 !!write(numout,*) 'erp(30,',jj,')=',erp(30,jj), & 123 !! & 'sf_sss(1)%fnow(30,',jj,',1)=',sf_sss(1)%fnow(30,jj,1), & 124 !! & 'sss_m(30,',jj,')=',sss_m(30,jj),'sst_m(30,',jj,')=',sst_m(30,jj), & 125 !! & 'rnfmsk_arcmsk(30,',jj,')=',rnfmsk_arcmsk(30,jj) 117 126 END DO 118 127 CALL iom_put( "erp", erp ) ! freshwater flux damping 128 CALL iom_put( "rnfmsk_arcmsk", rnfmsk_arcmsk ) ! river mough & arctic 119 129 ! 120 130 ELSEIF( nn_sssr == 2 ) THEN !* Salinity damping term (volume flux (emp) and associated heat flux (qns) … … 124 134 DO jj = 1, jpj 125 135 DO ji = 1, jpi 126 zerp = zsrp * ( 1. - 2.*rnfmsk (ji,jj) ) & ! No damping in vicinity of river mouths136 zerp = zsrp * ( 1. - 2.*rnfmsk_arcmsk(ji,jj) ) & ! No damping in vicinity of river mouths 127 137 & * ( sss_m(ji,jj) - sf_sss(1)%fnow(ji,jj,1) ) & 128 138 & / MAX( sss_m(ji,jj), 1.e-20 ) … … 134 144 END DO 135 145 CALL iom_put( "erp", erp ) ! freshwater flux damping 146 CALL iom_put( "rnfmsk_arcmsk", rnfmsk_arcmsk ) ! river mough & arctic 136 147 ENDIF 137 148 ! -
branches/UKMO/dev_r5518_haney_arctic_mask_for_tests/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r5431 r6913 21 21 USE sbcmod ! ln_rnf 22 22 USE sbcrnf ! River runoff 23 USE sbc_arcmsk ! River runoff & Arctic mask 23 24 USE sbcisf ! Ice shelf 24 25 USE traqsr ! solar radiation penetration -
branches/UKMO/dev_r5518_haney_arctic_mask_for_tests/NEMOGCM/NEMO/OPA_SRC/step_oce.F90
r5501 r6913 24 24 USE sbcmod ! surface boundary condition (sbc routine) 25 25 USE sbcrnf ! surface boundary condition: runoff variables 26 USE sbc_arcmsk ! surface boundary condition: runoff variables 26 27 USE sbccpl ! surface boundary condition: coupled formulation (call send at end of step) 27 28 USE sbc_oce ! surface boundary condition: ocean
Note: See TracChangeset
for help on using the changeset viewer.