Changeset 15804

Timestamp:

2022-04-27T14:02:36+02:00 (2 years ago)

Author:

jenniewaters

Message:

initial testing of inertial oscillation assimilation.

File:

• branches/UKMO/dev_r5518_assim_inertial_osci/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90 (modified) (4 diffs)

branches/UKMO/dev_r5518_assim_inertial_osci/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

@@ -35 +35 @@
    USE asmpar           ! Parameters for the assmilation interface
    USE c1d              ! 1D initialization
+   USE phycst,  ONLY: &
+      & omega
    USE in_out_manager   ! I/O manager
    USE lib_mpp          ! MPP library
@@ -876 +878 @@
       INTEGER :: jk
       INTEGER :: it
-      REAL(wp) :: zincwgt  ! IAU weight for current time step
+      REAL(wp) :: zincwgt    ! IAU weight for current time step
+      REAL(wp) :: z_delta_t  ! time shift from correct obs time
+      REAL (wp), DIMENSION(jpi,jpj,jpk) :: z_u_bkginc_rot
+      REAL (wp), DIMENSION(jpi,jpj,jpk) :: z_v_bkginc_rot
       !!----------------------------------------------------------------------
 
@@ -889 +894 @@
             it = kt - nit000 + 1
             zincwgt = wgtiau(it) / rdt   ! IAU weight for the current time step
+
+            z_delta_t=(real(kt) -49.)*900.
 
             IF(lwp) THEN
@@ -894 +901 @@
                WRITE(numout,*) 'dyn_asm_inc : Dynamics IAU at time step = ', &
                   &  kt,' with IAU weight = ', wgtiau(it)
+               WRITE(numout,*) 
+               WRITE(numout,*) ' z_delta_t = ',z_delta_t
                WRITE(numout,*) '~~~~~~~~~~~~'
             ENDIF
 
+            DO jk = 1, jpkm1
+               z_u_bkginc_rot(:,:,jk) = ( u_bkginc(:,:,jk)*COS(ff(:,:)*z_delta_t) ) + ( v_bkginc(:,:,jk)*SIN(ff(:,:)*z_delta_t) )
+               z_v_bkginc_rot(:,:,jk) = ( v_bkginc(:,:,jk)*COS(ff(:,:)*z_delta_t) ) - ( u_bkginc(:,:,jk)*SIN(ff(:,:)*z_delta_t) )
+            END DO
+
+
             ! Update the dynamic tendencies
             DO jk = 1, jpkm1
-               ua(:,:,jk) = ua(:,:,jk) + u_bkginc(:,:,jk) * zincwgt
-               va(:,:,jk) = va(:,:,jk) + v_bkginc(:,:,jk) * zincwgt
+               ua(:,:,jk) = ua(:,:,jk) + z_u_bkginc_rot(:,:,jk) * zincwgt
+               va(:,:,jk) = va(:,:,jk) + z_v_bkginc_rot(:,:,jk) * zincwgt
             END DO