New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 15804 – NEMO

Changeset 15804


Ignore:
Timestamp:
2022-04-27T14:02:36+02:00 (5 months ago)
Author:
jenniewaters
Message:

initial testing of inertial oscillation assimilation.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_assim_inertial_osci/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r14194 r15804  
    3535   USE asmpar           ! Parameters for the assmilation interface 
    3636   USE c1d              ! 1D initialization 
     37   USE phycst,  ONLY: & 
     38      & omega 
    3739   USE in_out_manager   ! I/O manager 
    3840   USE lib_mpp          ! MPP library 
     
    876878      INTEGER :: jk 
    877879      INTEGER :: it 
    878       REAL(wp) :: zincwgt  ! IAU weight for current time step 
     880      REAL(wp) :: zincwgt    ! IAU weight for current time step 
     881      REAL(wp) :: z_delta_t  ! time shift from correct obs time 
     882      REAL (wp), DIMENSION(jpi,jpj,jpk) :: z_u_bkginc_rot 
     883      REAL (wp), DIMENSION(jpi,jpj,jpk) :: z_v_bkginc_rot 
    879884      !!---------------------------------------------------------------------- 
    880885 
     
    889894            it = kt - nit000 + 1 
    890895            zincwgt = wgtiau(it) / rdt   ! IAU weight for the current time step 
     896 
     897            z_delta_t=(real(kt) -49.)*900. 
    891898 
    892899            IF(lwp) THEN 
     
    894901               WRITE(numout,*) 'dyn_asm_inc : Dynamics IAU at time step = ', & 
    895902                  &  kt,' with IAU weight = ', wgtiau(it) 
     903               WRITE(numout,*)  
     904               WRITE(numout,*) ' z_delta_t = ',z_delta_t 
    896905               WRITE(numout,*) '~~~~~~~~~~~~' 
    897906            ENDIF 
    898907 
     908            DO jk = 1, jpkm1 
     909               z_u_bkginc_rot(:,:,jk) = ( u_bkginc(:,:,jk)*COS(ff(:,:)*z_delta_t) ) + ( v_bkginc(:,:,jk)*SIN(ff(:,:)*z_delta_t) ) 
     910               z_v_bkginc_rot(:,:,jk) = ( v_bkginc(:,:,jk)*COS(ff(:,:)*z_delta_t) ) - ( u_bkginc(:,:,jk)*SIN(ff(:,:)*z_delta_t) ) 
     911            END DO 
     912 
     913 
    899914            ! Update the dynamic tendencies 
    900915            DO jk = 1, jpkm1 
    901                ua(:,:,jk) = ua(:,:,jk) + u_bkginc(:,:,jk) * zincwgt 
    902                va(:,:,jk) = va(:,:,jk) + v_bkginc(:,:,jk) * zincwgt 
     916               ua(:,:,jk) = ua(:,:,jk) + z_u_bkginc_rot(:,:,jk) * zincwgt 
     917               va(:,:,jk) = va(:,:,jk) + z_v_bkginc_rot(:,:,jk) * zincwgt 
    903918            END DO 
    904919            
Note: See TracChangeset for help on using the changeset viewer.