Changeset 10664 for branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/OPA_SRC/ASM/asmphyto2dbal_ersem.F90

Timestamp:

2019-02-13T11:23:00+01:00 (5 years ago)

Author:

dford

Message:

Update comment.

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/OPA_SRC/ASM/asmphyto2dbal_ersem.F90 (modified) (1 diff)

branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/OPA_SRC/ASM/asmphyto2dbal_ersem.F90

@@ -197 +197 @@
       !!---------------------------------------------------------------------------
       
-      ! Set parameters
+      ! Set parameters, firstly molar mass of carbon and nitrogen
       zmassc = 12.01
       zmassn = 14.01
+      
+      ! Then mesozooplankton nitrogen(mmol):carbon(mg), ERSEM parameter Z4/qnc
+      ! Hardwire it for now due to difficulty of getting at parameters through FABM
+      ! I think the following should work:
+      ! z4qnc = model%state_variables(jp_fabm_z4c)%properties%get_property_by_name('qnc')%value
+      ! but, get_property_by_name is in the module fabm_properties, which we can't get at directly
+      ! We can do a "USE fabm" which itself does a "USE fabm_properties", but doing it indirectly
+      ! like that means get_property_by_name is treated as private so we can't use it
+      ! So hardwire to value used in SSB runs and v4/5 CMEMS reanalysis
       z4qnc  = 0.0126
-      !z4qnc  = model%state_variables(jp_fabm_z4c)%parameters%qnc%value
-      !z4qnc  = get_property_by_name(model%state_variables(jp_fabm_z4c)%parameters, 'qnc')
-      IF (lwp) WRITE(numout,*) 'z4qnc = ', z4qnc
 
       ! If p_maxchlinc > 0 then cap total absolute chlorophyll increment at that value