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