Changeset 2001 for branches/dev_1784_OBS/NEMO/OPA_SRC/opa.F90

Timestamp:

2010-07-09T14:18:30+02:00 (14 years ago)

Author:

djlea

Message:

Adding observation operator code

File:

• branches/dev_1784_OBS/NEMO/OPA_SRC/opa.F90 (modified) (3 diffs)

branches/dev_1784_OBS/NEMO/OPA_SRC/opa.F90

@@ -54 +54 @@
    USE trdmod          ! momentum/tracers trends       (trd_mod_init routine)
    USE diaptr          ! poleward transports           (dia_ptr_init routine)
+   USE diaobs          ! Observation diagnostics       (dia_obs_init routine)
    USE step            ! OPA time-stepping                  (stp     routine)
 #if defined key_oasis3
@@ -143 +144 @@
          END DO
       ENDIF
+
+      IF( lk_diaobs ) CALL dia_obs_wri
        
       !                            !------------------------!
@@ -269 +272 @@
       !                                     ! diagnostics
       CALL iom_init                             ! iom_put initialization
+
+      IF( lk_diaobs ) THEN
+         CALL dia_obs_init                  ! Initialize observational data
+         CALL dia_obs( nit000 - 1 )         ! Observation operator for restart
+      ENDIF      
+
       CALL dia_ptr_init                         ! Poleward TRansports initialization
       CALL trd_mod_init                         ! Mixed-layer/Vorticity/Integral constraints trends