Changeset 4245 for branches/2013/dev_LOCEAN_CMCC_INGV_MERC_UKMO_2013/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

Timestamp:

2013-11-19T12:19:21+01:00 (10 years ago)

Author:

cetlod

Message:

dev_locean_cmcc_ingv_ukmo_merc : merge in the MERC_UKMO dev branch with trunk rev 4119

File:

• branches/2013/dev_LOCEAN_CMCC_INGV_MERC_UKMO_2013/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90 (modified) (10 diffs)

branches/2013/dev_LOCEAN_CMCC_INGV_MERC_UKMO_2013/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

@@ -21 +21 @@
    USE par_oce
    USE dom_oce                  ! Ocean space and time domain variables
+   USE obs_fbm, ONLY: ln_cl4    ! Class 4 diagnostic switch
    USE obs_read_prof            ! Reading and allocation of observations (Coriolis)
    USE obs_read_sla             ! Reading and allocation of SLA observations  
@@ -48 +49 @@
    PUBLIC dia_obs_init, &  ! Initialize and read observations
       &   dia_obs,      &  ! Compute model equivalent to observations
-      &   dia_obs_wri      ! Write model equivalent to observations
+      &   dia_obs_wri,  &  ! Write model equivalent to observations
+      &   dia_obs_dealloc  ! Deallocate dia_obs data
 
    !! * Shared Module variables
@@ -80 +82 @@
    LOGICAL, PUBLIC :: ln_ssh         !: Logical switch for sea surface height
    LOGICAL, PUBLIC :: ln_sss         !: Logical switch for sea surface salinity
+   LOGICAL, PUBLIC :: ln_sstnight    !: Logical switch for night mean SST observations
    LOGICAL, PUBLIC :: ln_nea         !: Remove observations near land
    LOGICAL, PUBLIC :: ln_altbias     !: Logical switch for altimeter bias  
@@ -167 +170 @@
          &            nmsshc, mdtcorr, mdtcutoff,                     &
          &            ln_reysst, ln_ghrsst, reysstname, reysstfmt,    &
+         &            ln_sstnight,                                    &
          &            ln_grid_search_lookup,                          &
          &            grid_search_file, grid_search_res,              &
@@ -176 +180 @@
          &            ln_velhradcp, velhradcpfiles,                   &
          &            ln_velfb, velfbfiles, ln_velfb_av,              &
-         &            ln_profb_enatim, ln_ignmis
+         &            ln_profb_enatim, ln_ignmis, ln_cl4
 
       INTEGER :: jprofset
@@ -311 +315 @@
          WRITE(numout,*) '             Logical switch for GHRSST observations          ln_ghrsst = ', ln_ghrsst
          WRITE(numout,*) '             Logical switch for feedback SST data             ln_sstfb = ', ln_sstfb
+         WRITE(numout,*) '             Logical switch for night-time SST obs         ln_sstnight = ', ln_sstnight
          WRITE(numout,*) '             Logical switch for SSS observations                ln_sss = ', ln_sss
          WRITE(numout,*) '             Logical switch for Sea Ice observations         ln_seaice = ', ln_seaice
@@ -704 +709 @@
             nsstsets = nsstsets + 1
 
-            ld_sstnight(nsstsets) = .TRUE.
+            ld_sstnight(nsstsets) = ln_sstnight
 
             CALL obs_rea_sst_rey( reysstname, reysstfmt, sstdata(nsstsets), &
@@ -718 +723 @@
             nsstsets = nsstsets + 1
 
-            ld_sstnight(nsstsets) = .TRUE.
+            ld_sstnight(nsstsets) = ln_sstnight
           
             CALL obs_rea_sst( 1, sstdata(nsstsets), jnumsst, &
@@ -737 +742 @@
                nsstsets = nsstsets + 1
 
-               ld_sstnight(nsstsets) = .TRUE.
+               ld_sstnight(nsstsets) = ln_sstnight
             
                CALL obs_rea_sst( 0, sstdata(nsstsets), 1, &
@@ -1368 +1373 @@
    END SUBROUTINE dia_obs_wri
 
+   SUBROUTINE dia_obs_dealloc
+      IMPLICIT NONE
+      !!----------------------------------------------------------------------
+      !!                    *** ROUTINE dia_obs_dealloc ***
+      !!
+      !!  ** Purpose : To deallocate data to enable the obs_oper online loop.
+      !!               Specifically: dia_obs_init --> dia_obs --> dia_obs_wri
+      !!
+      !!  ** Method : Clean up various arrays left behind by the obs_oper.
+      !!
+      !!  ** Action :
+      !!
+      !!----------------------------------------------------------------------
+      !! obs_grid deallocation
+      CALL obs_grid_deallocate
+
+      !! diaobs deallocation
+      IF ( nprofsets > 0 ) THEN
+          DEALLOCATE(ld_enact, &
+                  &  profdata, &
+                  &  prodatqc)
+      END IF
+      IF ( ln_sla ) THEN
+          DEALLOCATE(sladata, &
+                  &  sladatqc)
+      END IF
+      IF ( ln_seaice ) THEN
+          DEALLOCATE(sladata, &
+                  &  sladatqc)
+      END IF
+      IF ( ln_sst ) THEN
+          DEALLOCATE(sstdata, &
+                  &  sstdatqc)
+      END IF
+      IF ( ln_vel3d ) THEN
+          DEALLOCATE(ld_velav, &
+                  &  velodata, &
+                  &  veldatqc)
+      END IF
+   END SUBROUTINE dia_obs_dealloc
+
    SUBROUTINE ini_date( ddobsini )
       !!----------------------------------------------------------------------