- Timestamp:
- 2016-10-14T19:06:24+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_haney_arctic_mask_for_tests/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx_adj.F90
r6913 r7031 32 32 PUBLIC sbc_flx_adj_init ! routine called in sbcmod 33 33 34 !ky 06/09/2016 uncommented below two lines for FA test 34 35 !ky 11/12/2015 recommented below two lines 35 36 !!ky 3/12/2015 uncommented below two lines for FA test! 36 !REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: erp !: evaporation damping [kg/m2/s]37 !REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qrp !: heat flux damping [w/m2]37 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: erp !: evaporation damping [kg/m2/s] 38 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: qrp !: heat flux damping [w/m2] 38 39 39 40 ! !!* Namelist namsbc_flx_adj * … … 102 103 & * sf_qrp(1)%fnow(ji,jj,1) 103 104 !zqrp = sf_qrp(1)%fnow(ji,jj,1) 105 !ky 07/09/2016 copied 1 line below for FA test 104 106 !ky 11/12/2015 commented out 1 line below 105 107 !!ky 3/12/2015 1 line below for FA test! 106 108 !qrp(ji,jj) = sf_qrp(1)%fnow(ji,jj,1) 109 !ky 07/09/2016 copied 1 line above and modified as below for FA test 110 qrp(ji,jj) = zqrp 107 111 qns(ji,jj) = qns(ji,jj) + zqrp 108 112 END DO 109 113 END DO 114 !ky 06/09/2016 uncommented below 1 line below for FA test! 110 115 !ky 11/12/2015 recommented out below 1 line below 111 116 !!ky 3/12/2015 uncommented below 1 line below for FA test! 112 !CALL iom_put( "qrp", qrp ) ! heat flux damping117 CALL iom_put( "qrp", qrp ) ! heat flux damping 113 118 ENDIF 114 119 ! … … 127 132 !!ky 3/12/2015 1 line below for FA test! 128 133 !erp(ji,jj) = sf_erp(1)%fnow(ji,jj,1) 134 !ky 06/09/2016 copied line above and changed to below two lines for FA test! 135 erp(ji,jj) = ( 1. - 2.*rnfmsk_arcmsk(ji,jj) ) & ! No damping in vicinity of river mouths 136 & * sf_erp(1)%fnow(ji,jj,1) 129 137 sfx(ji,jj) = sfx(ji,jj) + zerp ! salt flux 130 138 END DO 131 139 END DO 140 !ky 06/09/2016 uncommented below one line for FA test! 132 141 !ky 11/12/2015 recommented out below one line 133 142 !!ky 3/12/2015 uncommented below one line for FA test! 134 !CALL iom_put( "erp", erp ) ! freshwater flux damping143 CALL iom_put( "erp", erp ) ! freshwater flux damping 135 144 ! 136 145 ELSEIF( nn_flxadjfw == 2 ) THEN !* Salinity damping term (volume flux (emp) and associated heat flux (qns) … … 145 154 !!ky 3/12/2015 1 line below for FA test! 146 155 !erp(ji,jj) = sf_erp(1)%fnow(ji,jj,1) 156 !ky 06/09/2016 copied line above and changed to below line for FA test! 157 erp(ji,jj) = zerp 147 158 emp(ji,jj) = emp (ji,jj) + zerp 148 159 qns(ji,jj) = qns(ji,jj) - zerp * rcp * sst_m(ji,jj) 149 160 END DO 150 161 END DO 162 !ky 06/09/2016 uncommented below 1 line for FA test! 151 163 !ky 11/12/2015 recommented out below 1 line 152 164 !!ky 3/12/2015 uncommented below 1 line for FA test! 153 !CALL iom_put( "erp", erp ) ! freshwater flux damping165 CALL iom_put( "erp", erp ) ! freshwater flux damping 154 166 ENDIF 155 167 ! … … 210 222 ! 211 223 ! !* Allocate erp and qrp array 224 !ky 06/09/2016 uncommented below two lines for FA test! 212 225 !ky 11/12/2015 recommented out below two lines 213 226 !!ky 3/12/2015 uncommented below two lines for FA test! 214 !ALLOCATE( qrp(jpi,jpj), erp(jpi,jpj), STAT=ierror )215 !IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_ssr: unable to allocate erp and qrp array' )227 ALLOCATE( qrp(jpi,jpj), erp(jpi,jpj), STAT=ierror ) 228 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_ssr: unable to allocate erp and qrp array' ) 216 229 ! 217 230 IF( nn_flxadjht == 1 ) THEN !* set sf_qrp structure & allocate arrays … … 244 257 ! 245 258 !ky!! !* Initialize qrp and erp if no restoring 259 !ky 06/09/2016 below two lines for FA test! 246 260 !ky 11/12/2015 commented out below two lines 247 261 !!ky 3/12/2015 below two lines for FA test! 248 !IF( nn_flxadjht /= 1 ) qrp(:,:) = 0._wp249 !IF( nn_flxadjfw /= 1 .OR. nn_flxadjfw /= 2 ) erp(:,:) = 0._wp262 IF( nn_flxadjht /= 1 ) qrp(:,:) = 0._wp 263 IF( nn_flxadjfw /= 1 .OR. nn_flxadjfw /= 2 ) erp(:,:) = 0._wp 250 264 !!ky!!IF( nn_sstr /= 1 ) qrp(:,:) = 0._wp 251 265 !!ky!!IF( nn_sssr /= 1 .OR. nn_sssr /= 2 ) erp(:,:) = 0._wp
Note: See TracChangeset
for help on using the changeset viewer.