Changeset 15359 for NEMO/branches/UKMO

Timestamp:

2021-10-13T11:49:46+02:00 (3 years ago)

Author:

jenniewaters

Message:

Include code for the pcbias

Location:

NEMO/branches/UKMO/NEMO_4.0.4_FOAM_pcbias/src/OCE

Files:

• ASM/bias.F90 (added)
• ASM/biaspar.F90 (added)
• nemogcm.F90 (modified) (1 diff)
• step.F90 (modified) (4 diffs)
• step_oce.F90 (modified) (1 diff)

NEMO/branches/UKMO/NEMO_4.0.4_FOAM_pcbias/src/OCE/nemogcm.F90

@@ -497 +497 @@
       !                                      ! Assimilation increments
       IF( lk_asminc    )   CALL asm_inc_init    ! Initialize assimilation increments
+                           CALL    bias_init    ! Pressure correction bias
       !
       IF(lwp) WRITE(numout,cform_aaa)           ! Flag AAAAAAA

NEMO/branches/UKMO/NEMO_4.0.4_FOAM_pcbias/src/OCE/step.F90

@@ -105 +105 @@
       IF( ln_crs         )   CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" )   ! tell IOM we are at time step kstp
 
+      IF( ln_bias )          CALL bias_opn( kstp )
+
       !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
       ! Update external forcing (tides, open boundaries, and surface boundary condition (including sea-ice)
@@ -230 +232 @@
       IF( ln_trabbl  )   CALL tra_bbl       ( kstp )  ! advective (and/or diffusive) bottom boundary layer scheme
       IF( ln_tradmp  )   CALL tra_dmp       ( kstp )  ! internal damping trends
+      IF( ln_bias    )   CALL tra_bias      ( kstp )  ! tracer bias
       IF( ln_bdy     )   CALL bdy_tra_dmp   ( kstp )  ! bdy damping trends
 #if defined key_agrif
@@ -262 +265 @@
 !!jc2: dynnxt must be the latest call. e3t_b are indeed updated in that routine
                          CALL tra_nxt       ( kstp )  ! finalize (bcs) tracer fields at next time step and swap
+      IF( ln_bias    )   CALL dyn_bias      ( kstp )
                          CALL dyn_nxt       ( kstp )  ! finalize (bcs) velocities at next time step and swap (always called after tra_nxt)
                          CALL ssh_swp       ( kstp )  ! swap of sea surface height
@@ -297 +301 @@
 
       IF( ln_diaobs  )   CALL dia_obs      ( kstp )      ! obs-minus-model (assimilation) diagnostics (call after dynamics update)
+      IF( lrst_bias )    CALL bias_wrt     ( kstp )      ! write bias fields
 
       !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

NEMO/branches/UKMO/NEMO_4.0.4_FOAM_pcbias/src/OCE/step_oce.F90

@@ -88 +88 @@
 
    USE asminc          ! assimilation increments      (tra_asm_inc routine)
-   !                                                   (dyn_asm_inc routine)
+                       !                              (dyn_asm_inc routine)
+   USE biaspar         ! bias param
+   USE bias            ! bias routines                    (tra_bias routine
+                       !                                  (dyn_bias routine)
    USE asmbkg
    USE stpctl          ! time stepping control            (stp_ctl routine)