New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 9032 for branches/NERC – NEMO

Changeset 9032 for branches/NERC


Ignore:
Timestamp:
2017-12-14T11:26:54+01:00 (5 years ago)
Author:
jpalmier
Message:

JPALM -- only the coupled Chl sent to OASIS is scaled - avoid all kind of scaling/unscaling on MEDUSA variable and associated errors

Location:
branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC
Files:
3 edited

Legend:

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

    r8671 r9032  
    244244            IF (chl_out.eq.1) THEN 
    245245               !! export and scale surface chl 
    246                zn_chl_srf(:,:) = MAX( 0.0, (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * scl_chl * 1.0E-6 ) 
     246               zn_chl_srf(:,:) = MAX( 0.0, (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * 1.0E-6 ) 
    247247                                 !! surf Chl in Kg-chl/m3 as needed for cpl 
    248248            ELSEIF (chl_out.eq.2) THEN 
    249249               !! export and scale mld chl 
    250                zn_chl_srf(:,:) = MAX( 0.0, fchl_ml(:,:) * scl_chl * 1.0E-6 ) 
     250               zn_chl_srf(:,:) = MAX( 0.0, fchl_ml(:,:) * 1.0E-6 ) 
    251251                                 !! mld Chl in Kg-chl/m3 as needed for cpl 
    252252            ENDIF 
    253             chloro_out_cpl(:,:) = zn_chl_srf(:,:)        !! Coupling Chl 
     253            chloro_out_cpl(:,:) = zn_chl_srf(:,:) * scl_chl        !! Coupling Chl 
    254254         END IF 
    255255 
     
    285285         ENDIF 
    286286         IF ( med_diag%CHL_CPL%dgsave ) THEN 
    287             CALL iom_put( "CHL_CPL"  , zn_chl_srf ) 
     287            CALL iom_put( "CHL_CPL"  , chloro_out_cpl ) 
    288288         ENDIF 
    289289         IF ( med_diag%PN_JLIM%dgsave ) THEN 
  • branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcini_medusa.F90

    r8955 r9032  
    271271      zn_dms_srf(:,:)  = 0.0 
    272272      za_dms_srf(:,:)  = 0.0 
    273       zn_chl_srf(:,:)  = 2.0E-8 * scl_chl !! Chl cpl - set first as surf 
     273      zn_chl_srf(:,:)  = 2.0E-8 !! Chl cpl - set first as surf 
    274274      !! 
    275275      IF(lwp) WRITE(numout,*) ' trc_ini_medusa: DMS and CO2 flux (UKESM) initialised to zero' 
     
    278278         CO2Flux_out_cpl(:,:) =  zn_co2_flx(:,:)   !! Coupling variable 
    279279         DMS_out_cpl(:,:)     =  zn_dms_srf(:,:)   !! Coupling variable 
    280          chloro_out_cpl(:,:)  =  zn_chl_srf(:,:)   !! Coupling variable 
     280         chloro_out_cpl(:,:)  =  zn_chl_srf(:,:) * scl_chl   !! Coupling variable 
    281281      END IF 
    282282      !! 
  • branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

    r9009 r9032  
    284284         IF(lwp) WRITE(numout,*) 'Chl cpl concentration - reading in ... - scale by ', scl_chl 
    285285         CALL iom_get( numrtr, jpdom_autoglo, 'N_CHL_srf',  zn_chl_srf(:,:)  ) 
    286          zn_chl_srf(:,:) =  zn_chl_srf(:,:) * scl_chl 
    287286      ELSE 
    288287         IF(lwp) WRITE(numout,*) 'set Chl coupled concentration - scaled by ', scl_chl 
    289          zn_chl_srf(:,:)  = MAX( 0.0, (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * scl_chl * 1.E-6 ) 
     288         zn_chl_srf(:,:)  = MAX( 0.0, (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * 1.E-6 ) 
    290289      ENDIF 
    291290      IF (lk_oasis) THEN 
    292          chloro_out_cpl(:,:) = zn_chl_srf(:,:)        !! Coupling variable 
     291         chloro_out_cpl(:,:) = zn_chl_srf(:,:) * scl_chl        !! Coupling variable 
    293292      END IF 
    294293      !! 
     
    302301      call trc_rst_dia_stat(zn_dms_srf(:,:), 'DMS surf') 
    303302      call trc_rst_dia_stat(zn_co2_flx(:,:), 'CO2 flux') 
    304       call trc_rst_dia_stat(zn_chl_srf(:,:), 'CHL  cpl') 
     303      call trc_rst_dia_stat(chloro_out_cpl(:,:), 'CHL  cpl') 
    305304      !!   
    306305      !! JPALM 14-06-2016 -- add Carbonate chenistry variables through the restart 
     
    465464      !!                     unscaled Chl field in the restart and only  
    466465      !!                     scale it when reading it in. 
    467       zn_chl_srf(:,:) = zn_chl_srf(:,:) / scl_chl  
    468466      CALL iom_rstput( kt, nitrst, numrtw, 'N_CHL_srf',  zn_chl_srf(:,:)  ) 
    469467      !! 
Note: See TracChangeset for help on using the changeset viewer.