Changeset 13899 for NEMO/branches/2020/tickets_icb_1900/src/OCE/OBS
- Timestamp:
- 2020-11-27T17:26:33+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/tickets_icb_1900
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/tickets_icb_1900
- Property svn:externals
-
NEMO/branches/2020/tickets_icb_1900/src/OCE/OBS/find_obs_proc.h90
r10068 r13899 41 41 ! first and last indoor i- and j-indexes kldi, klei, kldj, klej 42 42 ! exclude any obs in the bottom-left overlap region 43 ! also any obs outside to whole region (defined by nlci and nlcj)43 ! also any obs outside to whole region (defined by jpi and jpj) 44 44 ! I am assuming that kobsp does not need to be the correct processor 45 45 ! number -
NEMO/branches/2020/tickets_icb_1900/src/OCE/OBS/mpp_map.F90
r10068 r13899 11 11 !!---------------------------------------------------------------------- 12 12 USE par_kind, ONLY : wp ! Precision variables 13 USE par_oce , ONLY : jpi, jpj 14 USE dom_oce , ONLY : mig, mjg, n ldi, nlei, nldj, nlej, nlci, nlcj, narea! Ocean space and time domain variables13 USE par_oce , ONLY : jpi, jpj, Nis0, Nie0, Njs0, Nje0 ! Ocean parameters 14 USE dom_oce , ONLY : mig, mjg, narea ! Ocean space and time domain variables 15 15 #if defined key_mpp_mpi 16 USE lib_mpp , ONLY : mpi_comm_oce ! MPP library16 USE lib_mpp , ONLY : mpi_comm_oce ! MPP library 17 17 #endif 18 18 USE in_out_manager ! I/O manager … … 65 65 66 66 ! ! Setup local grid points 67 imppmap(mig(1):mig( nlci),mjg(1):mjg(nlcj)) = narea67 imppmap(mig(1):mig(jpi),mjg(1):mjg(jpj)) = narea 68 68 69 69 ! Get global data -
NEMO/branches/2020/tickets_icb_1900/src/OCE/OBS/obs_grid.F90
r12933 r13899 129 129 IF ( cdgrid == 'T' ) THEN 130 130 CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 131 & 1, nlci, 1, nlcj,&131 & 1, jpi, 1, jpj, & 132 132 & nproc, jpnij, & 133 133 & glamt, gphit, tmask, & … … 136 136 ELSEIF ( cdgrid == 'U' ) THEN 137 137 CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 138 & 1, nlci, 1, nlcj,&138 & 1, jpi, 1, jpj, & 139 139 & nproc, jpnij, & 140 140 & glamu, gphiu, umask, & … … 143 143 ELSEIF ( cdgrid == 'V' ) THEN 144 144 CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 145 & 1, nlci, 1, nlcj,&145 & 1, jpi, 1, jpj, & 146 146 & nproc, jpnij, & 147 147 & glamv, gphiv, vmask, & … … 150 150 ELSEIF ( cdgrid == 'F' ) THEN 151 151 CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 152 & 1, nlci, 1, nlcj,&152 & 1, jpi, 1, jpj, & 153 153 & nproc, jpnij, & 154 154 & glamf, gphif, fmask, & … … 279 279 zmskg(:,:) = -1.e+10 280 280 ! Add various grids here. 281 DO jj = 1, nlcj282 DO ji = 1, nlci281 DO jj = 1, jpj 282 DO ji = 1, jpi 283 283 zlamg(mig(ji),mjg(jj)) = glamt(ji,jj) 284 284 zphig(mig(ji),mjg(jj)) = gphit(ji,jj) … … 819 819 820 820 CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 821 & 1, nlci, 1, nlcj,&821 & 1, jpi, 1, jpj, & 822 822 & nproc, jpnij, & 823 823 & glamt, gphit, tmask, & -
NEMO/branches/2020/tickets_icb_1900/src/OCE/OBS/obs_oper.F90
r12377 r13899 189 189 ! Initialize daily mean for first timestep of the day 190 190 IF ( idayend == 1 .OR. kt == 0 ) THEN 191 DO_3D _11_11(1, jpk )191 DO_3D( 1, 1, 1, 1, 1, jpk ) 192 192 prodatqc%vdmean(ji,jj,jk,1) = 0.0 193 193 prodatqc%vdmean(ji,jj,jk,2) = 0.0 … … 195 195 ENDIF 196 196 197 DO_3D _11_11(1, jpk )197 DO_3D( 1, 1, 1, 1, 1, jpk ) 198 198 ! Increment field 1 for computing daily mean 199 199 prodatqc%vdmean(ji,jj,jk,1) = prodatqc%vdmean(ji,jj,jk,1) & … … 209 209 IF (lwp) WRITE(numout,*) 'Calculating prodatqc%vdmean on time-step: ',kt 210 210 CALL FLUSH(numout) 211 DO_3D _11_11(1, jpk )211 DO_3D( 1, 1, 1, 1, 1, jpk ) 212 212 prodatqc%vdmean(ji,jj,jk,1) = prodatqc%vdmean(ji,jj,jk,1) & 213 213 & * zdaystp … … 750 750 ! Initialize night-time mean for first timestep of the day 751 751 IF ( idayend == 1 .OR. kt == 0 ) THEN 752 DO_2D _11_11752 DO_2D( 1, 1, 1, 1 ) 753 753 surfdataqc%vdmean(ji,jj) = 0.0 754 754 zmeanday(ji,jj) = 0.0 … … 761 761 imask_night(:,:) = INT( zouttmp(:,:) ) 762 762 763 DO_2D _11_11763 DO_2D( 1, 1, 1, 1 ) 764 764 ! Increment the temperature field for computing night mean and counter 765 765 surfdataqc%vdmean(ji,jj) = surfdataqc%vdmean(ji,jj) & … … 773 773 IF ( idayend == 0 ) THEN 774 774 IF (lwp) WRITE(numout,*) 'Calculating surfdataqc%vdmean on time-step: ',kt 775 DO_2D _11_11775 DO_2D( 1, 1, 1, 1 ) 776 776 ! Test if "no night" point 777 777 IF ( icount_night(ji,jj) > 0 ) THEN -
NEMO/branches/2020/tickets_icb_1900/src/OCE/OBS/obs_read_altbias.F90
r12377 r13899 125 125 ! Get the Alt bias data 126 126 127 CALL iom_get( numaltbias, jpdom_ data, 'altbias', z_altbias(:,:), 1)127 CALL iom_get( numaltbias, jpdom_global, 'altbias', z_altbias(:,:) ) 128 128 129 129 ! Close the file -
NEMO/branches/2020/tickets_icb_1900/src/OCE/OBS/obs_readmdt.F90
r12377 r13899 90 90 CALL iom_open( mdtname, nummdt ) ! Open the file 91 91 ! ! Get the MDT data 92 CALL iom_get ( nummdt, jpdom_ data, 'sossheig', z_mdt(:,:), 1)92 CALL iom_get ( nummdt, jpdom_global, 'sossheig', z_mdt(:,:) ) 93 93 CALL iom_close(nummdt) ! Close the file 94 94 … … 215 215 zeta2 = 0.0 216 216 217 DO_2D _11_11217 DO_2D( 1, 1, 1, 1 ) 218 218 zdxdy = e1e2t(ji,jj) * zpromsk(ji,jj) 219 219 zarea = zarea + zdxdy -
NEMO/branches/2020/tickets_icb_1900/src/OCE/OBS/obs_sstbias.F90
r12377 r13899 139 139 cl_bias_files(jtype) ) 140 140 ! Get the SST bias data 141 CALL iom_get( numsstbias, jpdom_ data, 'tn', z_sstbias_2d(:,:), 1 )141 CALL iom_get( numsstbias, jpdom_global, 'tn', z_sstbias_2d(:,:), 1 ) 142 142 z_sstbias(:,:,jtype) = z_sstbias_2d(:,:) 143 143 ! Close the file
Note: See TracChangeset
for help on using the changeset viewer.