Changeset 2002 for branches/dev_1784_ASM/NEMO/OPA_SRC/opa.F90

Timestamp:

2010-07-09T14:19:55+02:00 (14 years ago)

Author:

djlea

Message:

Adding ASM assimilation code

File:

• branches/dev_1784_ASM/NEMO/OPA_SRC/opa.F90 (modified) (4 diffs)

branches/dev_1784_ASM/NEMO/OPA_SRC/opa.F90

@@ -25 +25 @@
    !!             -   ! 2007-07  (J. Chanut, A. Sellar) Unstructured open boundaries (BDY)
    !!            3.2  ! 2009-08  (S. Masson)  open/write in the listing file in mpp
+   !!            3.3  ! 2010-05  (K. Mogensen, A. Weaver, M. Martin, D. Lea) Assimilation interface 
    !!----------------------------------------------------------------------
 
@@ -53 +54 @@
    USE phycst          ! physical constant                  (par_cst routine)
    USE trdmod          ! momentum/tracers trends       (trd_mod_init routine)
+   USE asminc          ! assimilation increments       (asm_inc_init routine)
+   USE asmtrj          ! writing out state trajectory
+   USE sshwzv          ! vertical velocity used in asm
    USE diaptr          ! poleward transports           (dia_ptr_init routine)
    USE step            ! OPA time-stepping                  (stp     routine)
@@ -133 +137 @@
          END DO
       ELSE                              !==  3D ocean with  ==!
+          IF( lk_asminc ) THEN
+             IF( ln_bkgwri ) CALL asm_bkg_wri( nit000 - 1 )    ! Output background fields
+             IF( ln_trjwri ) CALL asm_trj_wri( nit000 - 1 )    ! Output trajectory fields
+             IF( ln_asmdin ) THEN                        ! Direct initialization
+                IF( ln_trainc ) CALL tra_asm_inc( nit000 - 1 )    ! Tracers
+                IF( ln_dyninc ) THEN 
+                   CALL dyn_asm_inc( nit000 - 1 )    ! Dynamics
+                   IF ( ln_asmdin ) CALL ssh_wzv ( nit000 - 1 )      ! update vertical velocity 
+                ENDIF
+                IF( ln_sshinc ) CALL ssh_asm_inc( nit000 - 1 )    ! SSH
+             ENDIF
+          ENDIF
+        
          DO WHILE ( istp <= nitend .AND. nstop == 0 )
 #if defined key_agrif
@@ -271 +288 @@
       CALL dia_ptr_init                         ! Poleward TRansports initialization
       CALL trd_mod_init                         ! Mixed-layer/Vorticity/Integral constraints trends
+
+      IF( lk_asminc ) CALL asm_inc_init     ! Initialize assimilation increments
+      IF(lwp) WRITE(numout,*)'Euler time step switch is ', neuler
+
       !
    END SUBROUTINE opa_init