Changeset 3317 for branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90

Timestamp:

2012-02-23T12:21:08+01:00 (12 years ago)

Author:

gm

Message:

Ediag branche: #927 restructuration of the trdicp computation - part I

File:

• branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90 (modified) (2 diffs)

branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90

@@ -31 +31 @@
    USE dynadv          ! dynamics: vector invariant versus flux form
    USE domvvl          ! variable volume
+   USE trdmod         ! ocean dynamics trends 
+   USE trdmod_oce     ! ocean variables trends
    USE obc_oce         ! ocean open boundary conditions
    USE obcdyn          ! open boundary condition for momentum (obc_dyn routine)
@@ -292 +294 @@
 
       IF( ln_3D_trd_d ) THEN             ! 3D output: asselin filter trends on momentum
-        zua(:,:,:) = ( ub(:,:,:) - zua(:,:,:) ) * z1_2dt
-        zva(:,:,:) = ( vb(:,:,:) - zva(:,:,:) ) * z1_2dt
-        CALL iom_put( "utrd_atf", zua )        ! asselin filter trends on momentum 
-        CALL iom_put( "vtrd_atf", zva )
+         zua(:,:,:) = ( ub(:,:,:) - zua(:,:,:) ) * z1_2dt
+         zva(:,:,:) = ( vb(:,:,:) - zva(:,:,:) ) * z1_2dt
+         CALL trd_mod( zua, zva, jpdyn_trd_atf, 'DYN', kt )
       ENDIF