- Timestamp:
- 2019-05-16T15:23:56+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/closea.F90
r10888 r10986 97 97 !!---------------------------------------------------------------------- 98 98 ! 99 IF(lwp) WRITE(numout,*) 100 IF(lwp) WRITE(numout,*)'dom_clo : read in masks to define closed seas ' 101 IF(lwp) WRITE(numout,*)'~~~~~~~' 99 IF(lwp) THEN 100 WRITE(numout,*) 101 WRITE(numout,*)'dom_clo : read in masks to define closed seas ' 102 WRITE(numout,*)'~~~~~~~' 103 IF(lflush) CALL FLUSH(numout) 104 ENDIF 102 105 ! 103 106 ! read the closed seas masks (if they exist) from domain_cfg file (if it exists) … … 120 123 CALL mpp_max('closea', jncs) 121 124 IF( jncs > 0 ) THEN 122 IF( lwp ) WRITE(numout,*) 'Number of closed seas : ',jncs 125 IF( lwp ) THEN 126 WRITE(numout,*) 'Number of closed seas : ',jncs 127 IF(lflush) CALL FLUSH(numout) 128 ENDIF 123 129 ELSE 124 130 CALL ctl_stop( 'Problem with closea_mask field in domain_cfg file. Has no values > 0 so no closed seas defined.') 125 131 ENDIF 126 132 ELSE 127 IF( lwp ) WRITE(numout,*) 128 IF( lwp ) WRITE(numout,*) ' ==>>> closea_mask field not found in domain_cfg file.' 129 IF( lwp ) WRITE(numout,*) ' No closed seas defined.' 130 IF( lwp ) WRITE(numout,*) 133 IF( lwp ) THEN 134 WRITE(numout,*) 135 WRITE(numout,*) ' ==>>> closea_mask field not found in domain_cfg file.' 136 WRITE(numout,*) ' No closed seas defined.' 137 WRITE(numout,*) 138 IF(lflush) CALL FLUSH(numout) 139 ENDIF 131 140 l_sbc_clo = .false. 132 141 jncs = 0 … … 148 157 CALL mpp_max('closea', jncsr) 149 158 IF( jncsr > 0 ) THEN 150 IF( lwp ) WRITE(numout,*) 'Number of closed seas rnf mappings : ',jncsr 159 IF( lwp ) THEN 160 WRITE(numout,*) 'Number of closed seas rnf mappings : ',jncsr 161 IF(lflush) CALL FLUSH(numout) 162 ENDIF 151 163 ELSE 152 164 CALL ctl_stop( 'Problem with closea_mask_rnf field in domain_cfg file. Has no values > 0 so no closed seas rnf mappings defined.') 153 165 ENDIF 154 166 ELSE 155 IF( lwp ) WRITE(numout,*) 'closea_mask_rnf field not found in domain_cfg file. No closed seas rnf mappings defined.' 167 IF( lwp ) THEN 168 WRITE(numout,*) 'closea_mask_rnf field not found in domain_cfg file. No closed seas rnf mappings defined.' 169 IF(lflush) CALL FLUSH(numout) 170 ENDIF 156 171 jncsr = 0 157 172 ENDIF … … 168 183 CALL mpp_max('closea', jncse) 169 184 IF( jncse > 0 ) THEN 170 IF( lwp ) WRITE(numout,*) 'Number of closed seas empmr mappings : ',jncse 185 IF( lwp ) THEN 186 WRITE(numout,*) 'Number of closed seas empmr mappings : ',jncse 187 IF(lflush) CALL FLUSH(numout) 188 ENDIF 171 189 ELSE 172 190 CALL ctl_stop( 'Problem with closea_mask_empmr field in domain_cfg file. Has no values > 0 so no closed seas empmr mappings defined.') 173 191 ENDIF 174 192 ELSE 175 IF( lwp ) WRITE(numout,*) 'closea_mask_empmr field not found in domain_cfg file. No closed seas empmr mappings defined.' 193 IF( lwp ) THEN 194 WRITE(numout,*) 'closea_mask_empmr field not found in domain_cfg file. No closed seas empmr mappings defined.' 195 IF(lflush) CALL FLUSH(numout) 196 ENDIF 176 197 jncse = 0 177 198 ENDIF … … 182 203 ! 183 204 ELSE ! ln_read_cfg = .false. so no domain_cfg file 184 IF( lwp ) WRITE(numout,*) 'No domain_cfg file so no closed seas defined.' 205 IF( lwp ) THEN 206 WRITE(numout,*) 'No domain_cfg file so no closed seas defined.' 207 IF(lflush) CALL FLUSH(numout) 208 ENDIF 185 209 l_sbc_clo = .false. 186 210 l_clo_rnf = .false. … … 219 243 IF( kt == nit000 ) THEN ! Initialisation ! 220 244 ! !------------------! 221 IF(lwp) WRITE(numout,*) 222 IF(lwp) WRITE(numout,*)'sbc_clo : closed seas ' 223 IF(lwp) WRITE(numout,*)'~~~~~~~' 245 IF(lwp) THEN 246 WRITE(numout,*) 247 WRITE(numout,*)'sbc_clo : closed seas ' 248 WRITE(numout,*)'~~~~~~~' 249 IF(lflush) CALL FLUSH(numout) 250 ENDIF 224 251 225 252 ALLOCATE( surf(jncs+1) , STAT=ierr ) … … 267 294 ENDIF 268 295 ! 269 IF(lwp) WRITE(numout,*)' Closed sea surface areas (km2)' 270 DO jc = 1, jncs 271 IF(lwp) WRITE(numout,FMT='(1I3,5X,ES12.2)') jc, surf(jc) * 1.0e-6 272 END DO 273 IF(lwp) WRITE(numout,FMT='(A,ES12.2)') 'Global surface area excluding closed seas (km2): ', surf(jncs+1) * 1.0e-6 274 ! 275 IF(jncsr > 0) THEN 276 IF(lwp) WRITE(numout,*)' Closed sea target rnf surface areas (km2)' 277 DO jcr = 1, jncsr 278 IF(lwp) WRITE(numout,FMT='(1I3,5X,ES12.2)') jcr, surfr(jcr) * 1.0e-6 296 IF(lwp) THEN 297 WRITE(numout,*)' Closed sea surface areas (km2)' 298 DO jc = 1, jncs 299 WRITE(numout,FMT='(1I3,5X,ES12.2)') jc, surf(jc) * 1.0e-6 279 300 END DO 280 ENDIF 281 ! 282 IF(jncse > 0) THEN 283 IF(lwp) WRITE(numout,*)' Closed sea target empmr surface areas (km2)' 284 DO jce = 1, jncse 285 IF(lwp) WRITE(numout,FMT='(1I3,5X,ES12.2)') jce, surfe(jce) * 1.0e-6 286 END DO 301 WRITE(numout,FMT='(A,ES12.2)') 'Global surface area excluding closed seas (km2): ', surf(jncs+1) * 1.0e-6 302 ! 303 IF(jncsr > 0) THEN 304 WRITE(numout,*)' Closed sea target rnf surface areas (km2)' 305 DO jcr = 1, jncsr 306 WRITE(numout,FMT='(1I3,5X,ES12.2)') jcr, surfr(jcr) * 1.0e-6 307 END DO 308 ENDIF 309 ! 310 IF(jncse > 0) THEN 311 WRITE(numout,*)' Closed sea target empmr surface areas (km2)' 312 DO jce = 1, jncse 313 WRITE(numout,FMT='(1I3,5X,ES12.2)') jce, surfe(jce) * 1.0e-6 314 END DO 315 ENDIF 316 IF(lflush) CALL FLUSH(numout) 287 317 ENDIF 288 318 ENDIF … … 453 483 WRITE(numout,*) 'clo_bat : suppression of closed seas' 454 484 WRITE(numout,*) '~~~~~~~' 485 IF(lflush) CALL FLUSH(numout) 455 486 ENDIF 456 487 ! … … 461 492 id = iom_varid(inum, 'closea_mask', ldstop = .false.) 462 493 IF( id > 0 ) THEN 463 IF( lwp ) WRITE(numout,*) 'Suppressing closed seas in bathymetry based on closea_mask field,' 494 IF( lwp ) THEN 495 WRITE(numout,*) 'Suppressing closed seas in bathymetry based on closea_mask field,' 496 IF(lflush) CALL FLUSH(numout) 497 ENDIF 464 498 CALL iom_get ( inum, jpdom_data, 'closea_mask', zdata_in ) 465 499 closea_mask(:,:) = NINT(zdata_in(:,:)) … … 469 503 ENDWHERE 470 504 ELSE 471 IF( lwp ) WRITE(numout,*) 'No closea_mask field found in domain_cfg file. No suppression of closed seas.' 505 IF( lwp ) THEN 506 WRITE(numout,*) 'No closea_mask field found in domain_cfg file. No suppression of closed seas.' 507 IF(lflush) CALL FLUSH(numout) 508 ENDIF 472 509 ENDIF 473 510 ! … … 475 512 ! 476 513 ELSE 477 IF( lwp ) WRITE(numout,*) 'No domain_cfg file => no suppression of closed seas.' 514 IF( lwp ) THEN 515 WRITE(numout,*) 'No domain_cfg file => no suppression of closed seas.' 516 IF(lflush) CALL FLUSH(numout) 517 ENDIF 478 518 ENDIF 479 519 !
Note: See TracChangeset
for help on using the changeset viewer.