Changeset 13576 for branches/UKMO/AMM15_v3_6_STABLE_package_collate/NEMOGCM/NEMO/OPA_SRC/step.F90

Timestamp:

2020-10-09T12:35:11+02:00 (4 years ago)

Author:

dford

Message:

Update NEMO-FABM coupler for FABM v1, and introduce two-way NEMO-ERSEM coupling options. See ​https://code.metoffice.gov.uk/trac/utils/ticket/366.

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_collate/NEMOGCM/NEMO/OPA_SRC/step.F90 (modified) (1 diff)

branches/UKMO/AMM15_v3_6_STABLE_package_collate/NEMOGCM/NEMO/OPA_SRC/step.F90

@@ -96 +96 @@
                       CALL iom_init(      cxios_context          )  ! iom_put initialization
          IF( ln_crs ) CALL iom_init( TRIM(cxios_context)//"_crs" )  ! initialize context for coarse grid
+         
+#if defined key_fabm
+         ! FABM can only finish initialising once IOM has
+         IF ( lk_fabm ) CALL nemo_fabm_start
+#endif
+         
+         ! First call to dia_* and asm_inc_init must wait for FABM to be initialised
+         ! (if running with FABM, but no harm moving the calls to here from nemo_init either way)
+         CALL dia_obs( nit000 - 1 )   ! Observation operator for restart
+         
+         IF( lk_asminc ) THEN 
+#if defined key_shelf
+            CALL  zdf_mxl_tref()     ! Initialization of hmld_tref
+#endif 
+            CALL asm_inc_init     ! Initialize assimilation increments 
+            IF( ln_bkgwri ) CALL asm_bkg_wri( nit000 - 1 )    ! Output background fields
+            IF( ln_asmdin ) THEN                        ! Direct initialization
+               IF( ln_trainc ) CALL tra_asm_inc( nit000 - 1 )    ! Tracers
+               IF( ln_dyninc ) CALL dyn_asm_inc( nit000 - 1 )    ! Dynamics
+               IF( ln_sshinc ) CALL ssh_asm_inc( nit000 - 1 )    ! SSH
+               IF( lk_bgcinc ) CALL bgc_asm_inc( nit000 - 1 )    ! BGC
+            ENDIF
+         ENDIF
+         
+         CALL dia_tmb_init        ! TMB outputs
+         CALL dia_25h_init        ! 25h mean  outputs
+         CALL dia_diaopfoam_init  ! FOAM operational output
       ENDIF