- Timestamp:
- 2020-02-27T12:27:50+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/OCE/SBC/sbcssm.F90
r10888 r12475 59 59 !!--------------------------------------------------------------------- 60 60 ! 61 write(numout,*)'sbcssm 1: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 61 62 ! !* surface T-, U-, V- ocean level variables (T, S, depth, velocity) 62 63 DO jj = 1, jpj … … 83 84 ! 84 85 frq_m(:,:) = fraqsr_1lev(:,:) 86 write(numout,*)'sbcssm 2: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 85 87 ! 86 88 ELSE … … 106 108 ! 107 109 frq_m(:,:) = zcoef * fraqsr_1lev(:,:) 110 write(numout,*)'sbcssm 3: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 108 111 ! ! ---------------------------------------- ! 109 112 ELSEIF( MOD( kt - 2 , nn_fsbc ) == 0 ) THEN ! Initialisation: New mean computation ! … … 116 119 e3t_m(:,:) = 0._wp 117 120 frq_m(:,:) = 0._wp 121 write(numout,*)'sbcssm 4: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 118 122 ENDIF 119 123 ! ! ---------------------------------------- ! … … 121 125 ! ! ---------------------------------------- ! 122 126 ssu_m(:,:) = ssu_m(:,:) + ub(:,:,1) 127 write(numout,*)'sbcssm 5: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 123 128 ssv_m(:,:) = ssv_m(:,:) + vb(:,:,1) 124 129 IF( l_useCT ) THEN ; sst_m(:,:) = sst_m(:,:) + eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) … … 146 151 e3t_m(:,:) = e3t_m(:,:) * zcoef ! mean vertical scale factor [m] 147 152 frq_m(:,:) = frq_m(:,:) * zcoef ! mean fraction of solar net radiation absorbed in the 1st T level [-] 153 write(numout,*)'sbcssm 6: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 148 154 ! 149 155 ENDIF … … 167 173 ! 168 174 IF( lwxios ) CALL iom_swap( cxios_context ) 175 write(numout,*)'sbcssm 7: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 169 176 ENDIF 170 177 ! … … 180 187 CALL iom_put( 'frq_m', frq_m ) 181 188 ENDIF 189 190 write(numout,*)'sbcssm 8: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 191 182 192 ! 183 193 END SUBROUTINE sbc_ssm … … 195 205 !!---------------------------------------------------------------------- 196 206 ! 207 208 write(numout,*)'sbcssm 1: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 209 197 210 IF( nn_fsbc == 1 ) THEN 198 211 ! … … 216 229 CALL iom_get( numror, jpdom_autoglo, 'ssh_m' , ssh_m, ldxios = lrxios ) ! " " height (T-point) 217 230 CALL iom_get( numror, jpdom_autoglo, 'e3t_m' , e3t_m, ldxios = lrxios ) ! 1st level thickness (T-point) 231 232 write(numout,*)'sbcssm 2: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 233 218 234 ! fraction of solar net radiation absorbed in 1st T level 219 235 IF( iom_varid( numror, 'frq_m', ldstop = .FALSE. ) > 0 ) THEN … … 233 249 e3t_m(:,:) = zcoef * e3t_m(:,:) 234 250 frq_m(:,:) = zcoef * frq_m(:,:) 251 write(numout,*)'sbcssm 3: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 235 252 ELSE 236 253 IF(lwp) WRITE(numout,*) ' mean fields read in the ocean restart file' … … 238 255 ENDIF 239 256 ENDIF 257 258 write(numout,*)'sbcssm 4: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 259 240 260 ! 241 261 IF( .NOT.l_ssm_mean ) THEN ! default initialisation. needed by iceistate … … 253 273 ! 254 274 ENDIF 275 276 write(numout,*)'sbcssm 5: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 277 255 278 ! 256 279 IF( .NOT. ln_traqsr ) fraqsr_1lev(:,:) = 1._wp ! default definition: qsr 100% in the fisrt level … … 267 290 ENDIF 268 291 292 write(numout,*)'sbcssm 6: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 293 269 294 END SUBROUTINE sbc_ssm_init 270 295
Note: See TracChangeset
for help on using the changeset viewer.