Changeset 13451 for branches/UKMO/AMM15_v3_6_STABLE_package_collate_utils366_fabmv1/NEMOGCM/NEMO/TOP_SRC/FABM/trcsms_fabm.F90

Timestamp:

2020-09-07T18:50:53+02:00 (4 years ago)

Author:

dford

Message:

Modifications to accommodate interaction between NEMO-FABM coupler only saving diagnostics requested in iodef.xml, and the assimilation and 25h/tmb diagnostic code.

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_collate_utils366_fabmv1/NEMOGCM/NEMO/TOP_SRC/FABM/trcsms_fabm.F90 (modified) (4 diffs)

branches/UKMO/AMM15_v3_6_STABLE_package_collate_utils366_fabmv1/NEMOGCM/NEMO/TOP_SRC/FABM/trcsms_fabm.F90

@@ -52 +52 @@
    PUBLIC   trc_sms_fabm_alloc      ! called by trcini_fabm.F90 module
    PUBLIC   trc_sms_fabm_check_mass ! called by trcwri_fabm.F90
+   PUBLIC   nemo_fabm_start
 
    ! Work arrays
@@ -556 +557 @@
       fabm_st2Db = fabm_st2Dn
 
-      ! Pointer to spectral heating term
-      swrad => model%get_data(model%get_interior_variable_id(standard_variables%net_rate_of_absorption_of_shortwave_energy_in_layer))
-
-
    END FUNCTION trc_sms_fabm_alloc
 
@@ -569 +566 @@
       DO jn=1,size(model%interior_diagnostic_variables)
          model%interior_diagnostic_variables(jn)%save = iom_use(model%interior_diagnostic_variables(jn)%name) &
-            .or. iom_use(TRIM(model%interior_diagnostic_variables(jn)%name)//'_VINT')
-      END DO
+            .or. iom_use(TRIM(model%interior_diagnostic_variables(jn)%name)//'_VINT') &
+            .or. iom_use(TRIM(model%interior_diagnostic_variables(jn)%name)//'25h') &
+            .or. iom_use('top_'//TRIM(model%interior_diagnostic_variables(jn)%name)) &
+            .or. iom_use('mid_'//TRIM(model%interior_diagnostic_variables(jn)%name)) &
+            .or. iom_use('bot_'//TRIM(model%interior_diagnostic_variables(jn)%name))
+      END DO
+      model%interior_diagnostic_variables(jp_fabm_o3ta)%save = .TRUE.
+      model%interior_diagnostic_variables(jp_fabm_o3ph)%save = .TRUE.
+      model%interior_diagnostic_variables(jp_fabm_o3pc)%save = .TRUE.
+      model%interior_diagnostic_variables(jp_fabm_xeps)%save = .TRUE.
+      model%interior_diagnostic_variables(jp_fabm_pgrow)%save = .TRUE.
+      model%interior_diagnostic_variables(jp_fabm_ploss)%save = .TRUE.
       DO jn=1,size(model%horizontal_diagnostic_variables)
-         model%horizontal_diagnostic_variables(jn)%save = iom_use(model%horizontal_diagnostic_variables(jn)%name)
+         model%horizontal_diagnostic_variables(jn)%save = iom_use(model%horizontal_diagnostic_variables(jn)%name) &
+            .or. iom_use(TRIM(model%horizontal_diagnostic_variables(jn)%name)//'25h')
       END DO
 
@@ -580 +588 @@
 
       started = .TRUE.
+
+      IF( ln_qsr_spec ) THEN
+         ! Pointer to spectral heating term
+         swrad => model%get_data(model%get_interior_variable_id(standard_variables%net_rate_of_absorption_of_shortwave_energy_in_layer))
+      ENDIF
+
    END SUBROUTINE