- Timestamp:
- 2017-12-18T13:16:44+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r8427 r9114 277 277 !! as proxy of org matter from the ocean 278 278 !! -- needed for the coupling with atm 279 !! 07-12-2017 -- To make things cleaner, we want to store an 280 !! unscaled Chl field in the restart and only 281 !! scale it when reading it in. 282 279 283 IF( iom_varid( numrtr, 'N_CHL_srf', ldstop = .FALSE. ) > 0 ) THEN 280 IF(lwp) WRITE(numout,*) 'Chl surf concentration - reading in ...'284 IF(lwp) WRITE(numout,*) 'Chl cpl concentration - reading in ... - scale by ', scl_chl 281 285 CALL iom_get( numrtr, jpdom_autoglo, 'N_CHL_srf', zn_chl_srf(:,:) ) 282 286 ELSE 283 IF(lwp) WRITE(numout,*) ' Chl surf concentration - setting to zero ...'284 zn_chl_srf(:,:) = (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * 1.E-6287 IF(lwp) WRITE(numout,*) 'set Chl coupled concentration - scaled by ', scl_chl 288 zn_chl_srf(:,:) = MAX( 0.0, (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * 1.E-6 ) 285 289 ENDIF 286 290 IF (lk_oasis) THEN 287 chloro_out_cpl(:,:) = zn_chl_srf(:,:) !! Coupling variable291 chloro_out_cpl(:,:) = zn_chl_srf(:,:) * scl_chl !! Coupling variable 288 292 END IF 289 293 !! … … 297 301 call trc_rst_dia_stat(zn_dms_srf(:,:), 'DMS surf') 298 302 call trc_rst_dia_stat(zn_co2_flx(:,:), 'CO2 flux') 299 call trc_rst_dia_stat(zn_chl_srf(:,:), 'CHL surf') 303 IF (lk_oasis) THEN 304 call trc_rst_dia_stat(chloro_out_cpl(:,:), 'CHL cpl') 305 END IF 300 306 !! 301 307 !! JPALM 14-06-2016 -- add Carbonate chenistry variables through the restart … … 457 463 CALL iom_rstput( kt, nitrst, numrtw, 'B_CO2_flx', zb_co2_flx(:,:) ) 458 464 CALL iom_rstput( kt, nitrst, numrtw, 'N_CO2_flx', zn_co2_flx(:,:) ) 465 !! JPALM 07-12-2017 -- To make things cleaner, we want to store an 466 !! unscaled Chl field in the restart and only 467 !! scale it when reading it in. 459 468 CALL iom_rstput( kt, nitrst, numrtw, 'N_CHL_srf', zn_chl_srf(:,:) ) 460 469 !! … … 468 477 call trc_rst_dia_stat(zn_dms_srf(:,:), 'DMS surf') 469 478 call trc_rst_dia_stat(zn_co2_flx(:,:), 'CO2 flux') 470 call trc_rst_dia_stat(zn_chl_srf(:,:), ' CHL surf')479 call trc_rst_dia_stat(zn_chl_srf(:,:), 'unscaled CHL cpl') 471 480 !! 472 481 IF(lwp) WRITE(numout,*) ' MEDUSA averaged prop. for dust and iron dep.'
Note: See TracChangeset
for help on using the changeset viewer.