Changeset 7646 for trunk/NEMOGCM/NEMO/OPA_SRC/OBS
- Timestamp:
- 2017-02-06T10:25:03+01:00 (7 years ago)
- Location:
- trunk/NEMOGCM/NEMO/OPA_SRC/OBS
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/OBS/obs_conv_functions.h90
r2287 r7646 71 71 72 72 !! * Arguments 73 REAL(KIND=wp) :: pft ! in situ temperature in degrees celcius73 REAL(KIND=wp) :: pft ! in situ temperature in degrees Celsius 74 74 REAL(KIND=wp) :: pfs ! salinity in psu 75 75 REAL(KIND=wp) :: pfp ! pressure in bars -
trunk/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90
r6140 r7646 1251 1251 1252 1252 ! At the end of the day get interpolated means 1253 IF ( idayend == 0 .AND. ldnightav ) THEN 1254 1255 ALLOCATE( & 1256 & zsurfm(2,2,isurf) & 1257 & ) 1258 1259 CALL obs_int_comm_2d( 2, 2, isurf, kpi, kpj, igrdi, igrdj, & 1260 & surfdataqc%vdmean(:,:), zsurfm ) 1261 1253 IF (ldnightav ) THEN 1254 IF ( idayend == 0 ) THEN 1255 1256 ALLOCATE( & 1257 & zsurfm(2,2,isurf) & 1258 & ) 1259 1260 CALL obs_int_comm_2d( 2, 2, isurf, kpi, kpj, igrdi, igrdj, & 1261 & surfdataqc%vdmean(:,:), zsurfm ) 1262 1263 ENDIF 1262 1264 ENDIF 1263 1265 … … 1322 1324 1323 1325 ! At the end of the day also deallocate night-time mean array 1324 IF ( idayend == 0 .AND. ldnightav ) THEN 1325 DEALLOCATE( & 1326 & zsurfm & 1327 & ) 1326 IF ( ldnightav ) THEN 1327 IF ( idayend == 0 ) THEN 1328 DEALLOCATE( & 1329 & zsurfm & 1330 & ) 1331 ENDIF 1328 1332 ENDIF 1329 1333 -
trunk/NEMOGCM/NEMO/OPA_SRC/OBS/obs_prep.F90
r6140 r7646 13 13 !! obs_sor : Sort the observation arrays 14 14 !!--------------------------------------------------------------------- 15 !! * Modules used 16 USE par_kind, ONLY : & ! Precision variables 17 & wp 15 USE par_kind, ONLY : wp ! Precision variables 18 16 USE in_out_manager ! I/O manager 19 17 USE obs_profiles_def ! Definitions for storage arrays for profiles … … 24 22 USE obs_inter_sup ! Interpolation support 25 23 USE obs_oper ! Observation operators 26 USE lib_mpp, ONLY : & 27 & ctl_warn, ctl_stop 24 USE lib_mpp, ONLY : ctl_warn, ctl_stop 28 25 29 26 IMPLICIT NONE 30 31 !! * Routine accessibility32 27 PRIVATE 33 28 34 PUBLIC & 35 & obs_pre_prof, & ! First level check and screening of profile obs 36 & obs_pre_surf, & ! First level check and screening of surface obs 37 & calc_month_len ! Calculate the number of days in the months of a year 29 PUBLIC obs_pre_prof ! First level check and screening of profile obs 30 PUBLIC obs_pre_surf ! First level check and screening of surface obs 31 PUBLIC calc_month_len ! Calculate the number of days in the months of a year 38 32 39 33 !!---------------------------------------------------------------------- … … 63 57 !! ! 2015-02 (M. Martin) Combined routine for surface types. 64 58 !!---------------------------------------------------------------------- 65 !! * Modules used66 USE domstp ! Domain: set the time-step67 59 USE par_oce ! Ocean parameters 68 USE dom_oce, ONLY : & ! Geographical information 69 & glamt, & 70 & gphit, & 71 & tmask, & 72 & nproc 60 USE dom_oce, ONLY : glamt, gphit, tmask, nproc ! Geographical information 73 61 !! * Arguments 74 62 TYPE(obs_surf), INTENT(INOUT) :: surfdata ! Full set of surface data 75 63 TYPE(obs_surf), INTENT(INOUT) :: surfdataqc ! Subset of surface data not failing screening 76 64 LOGICAL, INTENT(IN) :: ld_nea ! Switch for rejecting observation near land 77 ! ! * Local declarations65 ! 78 66 INTEGER :: iyea0 ! Initial date 79 67 INTEGER :: imon0 ! - (year, month, day, hour, minute) … … 94 82 INTEGER :: inlasobsmpp ! - close to land 95 83 INTEGER :: igrdobsmpp ! - fail the grid search 96 LOGICAL, DIMENSION(:), ALLOCATABLE :: & 97 & llvalid ! SLA data selection 84 LOGICAL, DIMENSION(:), ALLOCATABLE :: llvalid ! SLA data selection 98 85 INTEGER :: jobs ! Obs. loop variable 99 86 INTEGER :: jstp ! Time loop variable 100 87 INTEGER :: inrc ! Time index variable 101 102 IF(lwp) WRITE(numout,*)'obs_pre_surf : Preparing the surface observations...' 103 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 88 !!---------------------------------------------------------------------- 89 90 IF(lwp) WRITE(numout,*) 'obs_pre_surf : Preparing the surface observations...' 91 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 104 92 105 93 ! Initial date initialization (year, month, day, hour, minute) … … 253 241 !! 254 242 !!---------------------------------------------------------------------- 255 !! * Modules used 256 USE domstp ! Domain: set the time-step 257 USE par_oce ! Ocean parameters 258 USE dom_oce, ONLY : & ! Geographical information 259 & gdept_1d, & 260 & nproc 243 USE par_oce ! Ocean parameters 244 USE dom_oce, ONLY : gdept_1d, nproc ! Geographical information 261 245 262 246 !! * Arguments … … 314 298 INTEGER :: jstp ! Time loop variable 315 299 INTEGER :: inrc ! Time index variable 300 !!---------------------------------------------------------------------- 316 301 317 302 IF(lwp) WRITE(numout,*)'obs_pre_prof: Preparing the profile data...'
Note: See TracChangeset
for help on using the changeset viewer.