Changeset 15799 for NEMO/branches/UKMO/NEMO_4.0.4_FOAM_package/src/OCE/TRA/tradmp.F90

Timestamp:

2022-04-25T17:15:21+02:00 (2 years ago)

Author:

dford

Message:

More generic interface and structure for OBS code. See Met Office utils tickets 471 and 530.

File:

• NEMO/branches/UKMO/NEMO_4.0.4_FOAM_package/src/OCE/TRA/tradmp.F90 (modified) (3 diffs)

NEMO/branches/UKMO/NEMO_4.0.4_FOAM_package/src/OCE/TRA/tradmp.F90

@@ -50 +50 @@
    !
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   resto    !: restoring coeff. on T and S (s-1)
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   tclim    !: temperature climatology on each time step(Celcius)
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   sclim    !: salinity climatology on each time step (psu)
 
    !! * Substitutions
@@ -64 +66 @@
       !!                ***  FUNCTION tra_dmp_alloc  ***
       !!----------------------------------------------------------------------
-      ALLOCATE( resto(jpi,jpj,jpk), STAT= tra_dmp_alloc )
+      ALLOCATE( resto(jpi,jpj,jpk), &
+         &      tclim(jpi,jpj,jpk), &
+         &      sclim(jpi,jpj,jpk), STAT= tra_dmp_alloc )
       !
       CALL mpp_sum ( 'tradmp', tra_dmp_alloc )
@@ -106 +110 @@
       CALL dta_tsd( kt, zts_dta )            ! read and interpolates T-S data at kt
       !
+      tclim(:,:,:) = zts_dta(:,:,:,jp_tem)
+      sclim(:,:,:) = zts_dta(:,:,:,jp_sal)
+      !
       SELECT CASE ( nn_zdmp )     !==  type of damping  ==!
       !