Changeset 7837 for branches/UKMO/dev_r5785_SSS_obsoper/NEMOGCM/NEMO/OPA_SRC/OBS/obs_sstbias.F90

Timestamp:

2017-03-27T10:50:24+02:00 (7 years ago)

Author:

mattmartin

Message:

First version of generic obs oper which works at NEMO3.6 with other GO6 branches.

File:

• branches/UKMO/dev_r5785_SSS_obsoper/NEMOGCM/NEMO/OPA_SRC/OBS/obs_sstbias.F90 (modified) (12 diffs)

branches/UKMO/dev_r5785_SSS_obsoper/NEMOGCM/NEMO/OPA_SRC/OBS/obs_sstbias.F90

@@ -49 +49 @@
       !! History : 
       !!      ! :  2014-08 (J. While) Bias correction code for SST obs,
-      !!      !                      based on obs_rea_altbias
+      !!      !                       based on obs_rea_altbias
       !!----------------------------------------------------------------------
       !! * Modules used
@@ -69 +69 @@
       INTEGER :: jpisstbias   ! Number of grid point in latitude for the bias
       INTEGER :: jpjsstbias   ! Number of grid point in longitude for the bias
-      INTEGER :: iico         ! Grid point indices
-      INTEGER :: ijco
       INTEGER :: jt
       INTEGER :: i_nx_id      ! Index to read the NetCDF file
@@ -129 +127 @@
          CALL iom_open( cl_bias_files(jtype), numsstbias, ldstop=.FALSE. )
 
-         IF (numsstbias .GT. 0) THEN
+         IF (numsstbias > 0) THEN
      
             !Read the bias type from the file
@@ -137 +135 @@
             iret=NF90_OPEN(TRIM(cl_bias_files(jtype)), NF90_NOWRITE, incfile)
             iret=NF90_GET_ATT( incfile, NF90_GLOBAL, "SST_source", &
-                              ifile_source )
+                               ifile_source )
             ibiastypes(jtype) = ifile_source                 
             iret=NF90_CLOSE(incfile)       
@@ -151 +149 @@
             ! Close the file
             CALL iom_close(numsstbias)
-     
+            IF ( lwp ) WRITE(numout,*) 'Read in bias for type: ',ibiastypes(jtype)
          ELSE
             CALL ctl_stop('obs_read_sstbias: File '// & 
@@ -178 +176 @@
       END DO
 
-      CALL obs_int_comm_2d( 2, 2, sstdata%nsurf, &
+      CALL obs_int_comm_2d( 2, 2, sstdata%nsurf, jpi, jpj, &
          &                  igrdi, igrdj, glamt, zglam )
-      CALL obs_int_comm_2d( 2, 2, sstdata%nsurf, &
+      CALL obs_int_comm_2d( 2, 2, sstdata%nsurf, jpi, jpj, &
          &                  igrdi, igrdj, gphit, zgphi )
-      CALL obs_int_comm_2d( 2, 2, sstdata%nsurf, &
+      CALL obs_int_comm_2d( 2, 2, sstdata%nsurf, jpi, jpj, &
          &                  igrdi, igrdj, tmask(:,:,1), zmask )
 
       DO jtype = 1, knumtypes
          
-         !Find the number observations of type
-         !and alllocate tempory arrays
+         !Find the number observations of this type
+         !and allocate tempory arrays
          inumtype = COUNT( sstdata%ntyp(:) == ibiastypes(jtype) )
 
@@ -199 +197 @@
             & zbias( 2,2,inumtype ) )
 
+         IF ( lwp ) WRITE(numout,*) 'Number of obs of type: ',ibiastypes(jtype),' is: ',inumtype
+         
          jt=1
          DO jobs = 1, sstdata%nsurf 
@@ -207 +207 @@
                zglam_tmp(:,:,jt) = zglam(:,:,jobs)
                zgphi_tmp(:,:,jt) = zgphi(:,:,jobs)
-               zgphi_tmp(:,:,jt) = zgphi(:,:,jobs)
                zmask_tmp(:,:,jt) = zmask(:,:,jobs)
 
@@ -215 +214 @@
          END DO
                          
-         CALL obs_int_comm_2d( 2, 2, inumtype, &
+         CALL obs_int_comm_2d( 2, 2, inumtype, jpi, jpj,       &
                &           igrdi_tmp(:,:,:), igrdj_tmp(:,:,:), &
                &           z_sstbias(:,:,jtype), zbias(:,:,:) )
@@ -225 +224 @@
                zlam = sstdata%rlam(jobs)
                zphi = sstdata%rphi(jobs)
-               iico = sstdata%mi(jobs)
-               ijco = sstdata%mj(jobs)         
-
-               CALL obs_int_h2d_init( 1, 1, k2dint, zlam, zphi,         &
-                  &                   zglam_tmp(:,:,jt), &
-                  &                   zgphi_tmp(:,:,jt), &
+
+               CALL obs_int_h2d_init( 1, 1, k2dint, zlam, zphi,             &
+                  &                   zglam_tmp(:,:,jt), zgphi_tmp(:,:,jt), &
                   &                   zmask_tmp(:,:,jt), zweig, zobsmask )
 
@@ -236 +232 @@
 
                ! adjust sst with bias field
-               sstdata%robs(jobs,1) = &
-                  &    sstdata%robs(jobs,1) - zext(1)
+               sstdata%robs(jobs,1) = sstdata%robs(jobs,1) - zext(1)
 
                jt=jt+1
@@ -243 +238 @@
             ENDIF
          END DO 
-               
+
+         IF ( lwp ) THEN
+            WRITE(numout,*) 'Applied bias to obs of type: ',ibiastypes(jtype)
+            WRITE(numout,*) 'Max/min of bias: ', maxval(zbias(:,:,:)), minval(zbias(:,:,:))
+         ENDIF
+
          !Deallocate arrays
          DEALLOCATE( &