Ignore:
Timestamp:
2017-12-08T12:37:13+01:00 (3 years ago)
Author:
jpalmier
Message:

unscale-scale coupled chl read-write in the dump

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

    r8632 r8955  
    277277      !!                     as proxy of org matter from the ocean 
    278278      !!                  -- needed for the coupling with atm 
     279      !!       07-12-2017 -- To make things cleaner we want to pass ans read an  
     280      !!                     unscaled Chl field to the restart and scale it when 
     281      !!                     reading it. 
    279282      IF( iom_varid( numrtr, 'N_CHL_srf', ldstop = .FALSE. ) > 0 ) THEN 
    280          IF(lwp) WRITE(numout,*) 'Chl surf concentration - reading in ...' 
     283         IF(lwp) WRITE(numout,*) 'Chl cpl concentration - reading in ... - scale by ', scl_chl 
    281284         CALL iom_get( numrtr, jpdom_autoglo, 'N_CHL_srf',  zn_chl_srf(:,:)  ) 
     285         zn_chl_srf(:,:) =  zn_chl_srf(:,:) * scl_chl 
    282286      ELSE 
    283          IF(lwp) WRITE(numout,*) 'Chl surf concentration - setting to zero ...' 
     287         IF(lwp) WRITE(numout,*) 'set Chl coupled concentration - scaled by ', scl_chl 
    284288         zn_chl_srf(:,:)  = MAX( 0.0, (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * scl_chl * 1.E-6 ) 
    285289      ENDIF 
     
    297301      call trc_rst_dia_stat(zn_dms_srf(:,:), 'DMS surf') 
    298302      call trc_rst_dia_stat(zn_co2_flx(:,:), 'CO2 flux') 
    299       call trc_rst_dia_stat(zn_chl_srf(:,:), 'CHL surf') 
     303      call trc_rst_dia_stat(zn_chl_srf(:,:), 'CHL  cpl') 
    300304      !!   
    301305      !! JPALM 14-06-2016 -- add Carbonate chenistry variables through the restart 
     
    457461      CALL iom_rstput( kt, nitrst, numrtw, 'B_CO2_flx',  zb_co2_flx(:,:)  ) 
    458462      CALL iom_rstput( kt, nitrst, numrtw, 'N_CO2_flx',  zn_co2_flx(:,:)  ) 
     463      !! JPALM 07-12-2017 -- To make things cleaner we want to pass ans read an  
     464      !!                     unscaled Chl field to the restart and scale it when 
     465      !!                     reading it. 
     466      zn_chl_srf(:,:) = zn_chl_srf(:,:) / scl_chl  
    459467      CALL iom_rstput( kt, nitrst, numrtw, 'N_CHL_srf',  zn_chl_srf(:,:)  ) 
    460468      !! 
     
    468476      call trc_rst_dia_stat(zn_dms_srf(:,:), 'DMS surf') 
    469477      call trc_rst_dia_stat(zn_co2_flx(:,:), 'CO2 flux') 
    470       call trc_rst_dia_stat(zn_chl_srf(:,:), 'CHL surf') 
     478      call trc_rst_dia_stat(zn_chl_srf(:,:), 'unscaled CHL cpl') 
    471479      !! 
    472480      IF(lwp) WRITE(numout,*) ' MEDUSA averaged prop. for dust and iron dep.' 
Note: See TracChangeset for help on using the changeset viewer.