Changeset 11468 for branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_prof.F90

Timestamp:

2019-08-23T10:37:22+02:00 (5 years ago)

Author:

mattmartin

Message:

Merged changes to allow writing of climatological information to feedback files.

File:

• branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_prof.F90 (modified) (4 diffs)

branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_prof.F90

@@ -46 +46 @@
       &                     kvars, kextr, kstp, ddobsini, ddobsend, &
       &                     ldvar, ldignmis, ldsatt, &
-      &                     ldmod, kdailyavtypes )
+      &                     ldmod, ldclim, kdailyavtypes )
       !!---------------------------------------------------------------------
       !!
@@ -78 +78 @@
       LOGICAL, INTENT(IN) :: ldsatt     ! Compute salinity at all temperature points
       LOGICAL, INTENT(IN) :: ldmod      ! Initialize model from input data
+      LOGICAL, INTENT(IN) :: ldclim     ! Set flag to show climatology will be output
       REAL(dp), INTENT(IN) :: ddobsini  ! Obs. ini time in YYYYMMDD.HHMMSS
       REAL(dp), INTENT(IN) :: ddobsend  ! Obs. end time in YYYYMMDD.HHMMSS
@@ -500 +501 @@
       ENDIF
       CALL obs_prof_alloc( profdata, kvars, kextr, iprof, iv3dt, &
-         &                 kstp, jpi, jpj, jpk )
+         &                 kstp, jpi, jpj, jpk, ldclim )
 
       ! * Read obs/positions, QC, all variable and assign to profdata
@@ -697 +698 @@
                                  &                inpfiles(jj)%padd(ij,ji,1,jvar)
                            ENDIF
+                           IF ( profdata%lclim ) THEN
+                               profdata%var(jvar)%vclm(ivart(jvar)) = fbrmdi
+                           ENDIF                          
                            ! Count number of profile var1 data as function of type
                            itypvar( profdata%ntyp(iprof) + 1, jvar ) = &