Changeset 9192
 Timestamp:
 20180108T13:36:41+01:00 (3 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/UKMO/dev_r5518_obs_oper_update_bgc3d/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90
r9191 r9192 892 892 INTEGER :: jvar ! Variable number 893 893 INTEGER :: ji, jj ! Loop counters 894 REAL(wp) :: tiny 894 REAL(wp) :: tiny ! small number 895 895 REAL(wp), POINTER, DIMENSION(:,:,:) :: & 896 896 & zprofvar1, & ! Model values for 1st variable in a prof ob … … 907 907 & zgphi1, & ! Model latitudes for prof variable 1 908 908 & zgphi2 ! Model latitudes for prof variable 2 909 LOGICAL :: llog10 ! Perform log10 transform of variable 909 910 910 911 … … 1160 1161 !Defaults which might be changed 1161 1162 zsurfmask(:,:) = tmask(:,:,1) 1163 llog10 = .FALSE. 1162 1164 1163 1165 SELECT CASE ( TRIM(cobstypessurf(jtype)) ) … … 1205 1207 & ' but no biogeochemical model appears to have been defined' ) 1206 1208 #endif 1207 zsurfmask(:,:) = tmask(:,:,1) ! create a special mask to exclude certain things 1208 ! Take the log10 where we can, otherwise exclude 1209 tiny = 1.0e20 1210 WHERE(zsurfvar(:,:) > tiny .AND. zsurfvar(:,:) /= obfillflt ) 1211 zsurfvar(:,:) = LOG10(zsurfvar(:,:)) 1212 ELSEWHERE 1213 zsurfvar(:,:) = obfillflt 1214 zsurfmask(:,:) = 0 1215 END WHERE 1209 llog10 = .TRUE. 1216 1210 1217 1211 CASE('slchldia') … … 1229 1223 & ' but no biogeochemical model appears to have been defined' ) 1230 1224 #endif 1231 zsurfmask(:,:) = tmask(:,:,1) ! create a special mask to exclude certain things 1232 ! Take the log10 where we can, otherwise exclude 1233 tiny = 1.0e20 1234 WHERE(zsurfvar(:,:) > tiny .AND. zsurfvar(:,:) /= obfillflt ) 1235 zsurfvar(:,:) = LOG10(zsurfvar(:,:)) 1236 ELSEWHERE 1237 zsurfvar(:,:) = obfillflt 1238 zsurfmask(:,:) = 0 1239 END WHERE 1225 llog10 = .TRUE. 1240 1226 1241 1227 CASE('slchlnon') … … 1254 1240 & ' but no biogeochemical model appears to have been defined' ) 1255 1241 #endif 1256 zsurfmask(:,:) = tmask(:,:,1) ! create a special mask to exclude certain things 1257 ! Take the log10 where we can, otherwise exclude 1258 tiny = 1.0e20 1259 WHERE(zsurfvar(:,:) > tiny .AND. zsurfvar(:,:) /= obfillflt ) 1260 zsurfvar(:,:) = LOG10(zsurfvar(:,:)) 1261 ELSEWHERE 1262 zsurfvar(:,:) = obfillflt 1263 zsurfmask(:,:) = 0 1264 END WHERE 1242 llog10 = .TRUE. 1265 1243 1266 1244 CASE('slchldin') … … 1278 1256 & ' but no biogeochemical model appears to have been defined' ) 1279 1257 #endif 1280 zsurfmask(:,:) = tmask(:,:,1) ! create a special mask to exclude certain things 1281 ! Take the log10 where we can, otherwise exclude 1282 tiny = 1.0e20 1283 WHERE(zsurfvar(:,:) > tiny .AND. zsurfvar(:,:) /= obfillflt ) 1284 zsurfvar(:,:) = LOG10(zsurfvar(:,:)) 1285 ELSEWHERE 1286 zsurfvar(:,:) = obfillflt 1287 zsurfmask(:,:) = 0 1288 END WHERE 1258 llog10 = .TRUE. 1289 1259 1290 1260 CASE('slchlmic') … … 1302 1272 & ' but no biogeochemical model appears to have been defined' ) 1303 1273 #endif 1304 zsurfmask(:,:) = tmask(:,:,1) ! create a special mask to exclude certain things 1305 ! Take the log10 where we can, otherwise exclude 1306 tiny = 1.0e20 1307 WHERE(zsurfvar(:,:) > tiny .AND. zsurfvar(:,:) /= obfillflt ) 1308 zsurfvar(:,:) = LOG10(zsurfvar(:,:)) 1309 ELSEWHERE 1310 zsurfvar(:,:) = obfillflt 1311 zsurfmask(:,:) = 0 1312 END WHERE 1274 llog10 = .TRUE. 1313 1275 1314 1276 CASE('slchlnan') … … 1326 1288 & ' but no biogeochemical model appears to have been defined' ) 1327 1289 #endif 1328 zsurfmask(:,:) = tmask(:,:,1) ! create a special mask to exclude certain things 1329 ! Take the log10 where we can, otherwise exclude 1330 tiny = 1.0e20 1331 WHERE(zsurfvar(:,:) > tiny .AND. zsurfvar(:,:) /= obfillflt ) 1332 zsurfvar(:,:) = LOG10(zsurfvar(:,:)) 1333 ELSEWHERE 1334 zsurfvar(:,:) = obfillflt 1335 zsurfmask(:,:) = 0 1336 END WHERE 1290 llog10 = .TRUE. 1337 1291 1338 1292 CASE('slchlpic') … … 1350 1304 & ' but no biogeochemical model appears to have been defined' ) 1351 1305 #endif 1352 zsurfmask(:,:) = tmask(:,:,1) ! create a special mask to exclude certain things 1353 ! Take the log10 where we can, otherwise exclude 1354 tiny = 1.0e20 1355 WHERE(zsurfvar(:,:) > tiny .AND. zsurfvar(:,:) /= obfillflt ) 1356 zsurfvar(:,:) = LOG10(zsurfvar(:,:)) 1357 ELSEWHERE 1358 zsurfvar(:,:) = obfillflt 1359 zsurfmask(:,:) = 0 1360 END WHERE 1306 llog10 = .TRUE. 1361 1307 1362 1308 CASE('schltot') … … 1375 1321 & ' but no biogeochemical model appears to have been defined' ) 1376 1322 #endif 1377 zsurfmask(:,:) = tmask(:,:,1) ! create a special mask to exclude certain things1378 1323 1379 1324 CASE('sspm') … … 1453 1398 1454 1399 END SELECT 1400 1401 IF ( llog10 ) THEN 1402 ! Take the log10 where we can, otherwise exclude 1403 tiny = 1.0e20 1404 WHERE(zsurfvar(:,:) > tiny .AND. zsurfvar(:,:) /= obfillflt ) 1405 zsurfvar(:,:) = LOG10(zsurfvar(:,:)) 1406 ELSEWHERE 1407 zsurfvar(:,:) = obfillflt 1408 zsurfmask(:,:) = 0 1409 END WHERE 1410 ENDIF 1455 1411 1456 1412 CALL obs_surf_opt( surfdataqc(jtype), kstp, jpi, jpj, &
Note: See TracChangeset
for help on using the changeset viewer.