Changeset 10168
- Timestamp:
- 2018-10-03T16:38:42+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90
r9306 r10168 888 888 USE par_hadocc 889 889 USE had_bgc_const, ONLY: c2n_p 890 #elif defined key_medusa && defined key_foam_medusa890 #elif defined key_medusa 891 891 USE trc, ONLY : & ! MEDUSA variables 892 892 & trn … … 985 985 ! Chlorophyll from HadOCC 986 986 zprofvar(:,:,:,1) = HADOCC_CHL(:,:,:) 987 #elif defined key_medusa && defined key_foam_medusa987 #elif defined key_medusa 988 988 ! Add non-diatom and diatom chlorophyll from MEDUSA 989 989 zprofvar(:,:,:,1) = trn(:,:,:,jpchn) + trn(:,:,:,jpchd) … … 1023 1023 ! Chlorophyll from HadOCC 1024 1024 zprofvar(:,:,:,1) = HADOCC_CHL(:,:,:) 1025 #elif defined key_medusa && defined key_foam_medusa1025 #elif defined key_medusa 1026 1026 ! Add non-diatom and diatom chlorophyll from MEDUSA 1027 1027 zprofvar(:,:,:,1) = trn(:,:,:,jpchn) + trn(:,:,:,jpchd) … … 1039 1039 ! Dissolved inorganic nitrogen from HadOCC 1040 1040 zprofvar(:,:,:,1) = trn(:,:,:,jp_had_nut) 1041 #elif defined key_medusa && defined key_foam_medusa1041 #elif defined key_medusa 1042 1042 ! Dissolved inorganic nitrogen from MEDUSA 1043 1043 zprofvar(:,:,:,1) = trn(:,:,:,jpdin) … … 1054 1054 CALL ctl_stop( ' Trying to run psi4 observation operator', & 1055 1055 & ' but HadOCC does not simulate silicate' ) 1056 #elif defined key_medusa && defined key_foam_medusa1056 #elif defined key_medusa 1057 1057 ! Silicate from MEDUSA 1058 1058 zprofvar(:,:,:,1) = trn(:,:,:,jpsil) … … 1069 1069 CALL ctl_stop( ' Trying to run ppo4 observation operator', & 1070 1070 & ' but HadOCC does not simulate phosphate' ) 1071 #elif defined key_medusa && defined key_foam_medusa1071 #elif defined key_medusa 1072 1072 CALL ctl_stop( ' Trying to run ppo4 observation operator', & 1073 1073 & ' but MEDUSA does not simulate phosphate' ) … … 1084 1084 ! Dissolved inorganic carbon from HadOCC 1085 1085 zprofvar(:,:,:,1) = trn(:,:,:,jp_had_dic) 1086 #elif defined key_medusa && defined key_foam_medusa1086 #elif defined key_medusa 1087 1087 ! Dissolved inorganic carbon from MEDUSA 1088 1088 zprofvar(:,:,:,1) = trn(:,:,:,jpdic) … … 1099 1099 ! Alkalinity from HadOCC 1100 1100 zprofvar(:,:,:,1) = trn(:,:,:,jp_had_alk) 1101 #elif defined key_medusa && defined key_foam_medusa1101 #elif defined key_medusa 1102 1102 ! Alkalinity from MEDUSA 1103 1103 zprofvar(:,:,:,1) = trn(:,:,:,jpalk) … … 1114 1114 CALL ctl_stop( ' Trying to run pph observation operator', & 1115 1115 & ' but HadOCC has no pH diagnostic defined' ) 1116 #elif defined key_medusa && defined key_ foam_medusa1116 #elif defined key_medusa && defined key_roam 1117 1117 ! pH from MEDUSA 1118 1118 zprofvar(:,:,:,1) = f3_pH(:,:,:) … … 1129 1129 CALL ctl_stop( ' Trying to run po2 observation operator', & 1130 1130 & ' but HadOCC does not simulate oxygen' ) 1131 #elif defined key_medusa && defined key_foam_medusa1131 #elif defined key_medusa 1132 1132 ! Oxygen from MEDUSA 1133 1133 zprofvar(:,:,:,1) = trn(:,:,:,jpoxy) … … 1210 1210 ! Surface chlorophyll from HadOCC 1211 1211 zsurfvar(:,:) = HADOCC_CHL(:,:,1) 1212 #elif defined key_medusa && defined key_foam_medusa1212 #elif defined key_medusa 1213 1213 ! Add non-diatom and diatom surface chlorophyll from MEDUSA 1214 1214 zsurfvar(:,:) = trn(:,:,1,jpchn) + trn(:,:,1,jpchd) … … 1227 1227 CALL ctl_stop( ' Trying to run slchldia observation operator', & 1228 1228 & ' but HadOCC does not explicitly simulate diatoms' ) 1229 #elif defined key_medusa && defined key_foam_medusa1229 #elif defined key_medusa 1230 1230 ! Diatom surface chlorophyll from MEDUSA 1231 1231 zsurfvar(:,:) = trn(:,:,1,jpchd) … … 1243 1243 CALL ctl_stop( ' Trying to run slchlnon observation operator', & 1244 1244 & ' but HadOCC does not explicitly simulate non-diatoms' ) 1245 #elif defined key_medusa && defined key_foam_medusa1245 #elif defined key_medusa 1246 1246 ! Non-diatom surface chlorophyll from MEDUSA 1247 1247 zsurfvar(:,:) = trn(:,:,1,jpchn) … … 1260 1260 CALL ctl_stop( ' Trying to run slchldin observation operator', & 1261 1261 & ' but HadOCC does not explicitly simulate dinoflagellates' ) 1262 #elif defined key_medusa && defined key_foam_medusa1262 #elif defined key_medusa 1263 1263 CALL ctl_stop( ' Trying to run slchldin observation operator', & 1264 1264 & ' but MEDUSA does not explicitly simulate dinoflagellates' ) … … 1276 1276 CALL ctl_stop( ' Trying to run slchlmic observation operator', & 1277 1277 & ' but HadOCC does not explicitly simulate microphytoplankton' ) 1278 #elif defined key_medusa && defined key_foam_medusa1278 #elif defined key_medusa 1279 1279 CALL ctl_stop( ' Trying to run slchlmic observation operator', & 1280 1280 & ' but MEDUSA does not explicitly simulate microphytoplankton' ) … … 1292 1292 CALL ctl_stop( ' Trying to run slchlnan observation operator', & 1293 1293 & ' but HadOCC does not explicitly simulate nanophytoplankton' ) 1294 #elif defined key_medusa && defined key_foam_medusa1294 #elif defined key_medusa 1295 1295 CALL ctl_stop( ' Trying to run slchlnan observation operator', & 1296 1296 & ' but MEDUSA does not explicitly simulate nanophytoplankton' ) … … 1308 1308 CALL ctl_stop( ' Trying to run slchlpic observation operator', & 1309 1309 & ' but HadOCC does not explicitly simulate picophytoplankton' ) 1310 #elif defined key_medusa && defined key_foam_medusa1310 #elif defined key_medusa 1311 1311 CALL ctl_stop( ' Trying to run slchlpic observation operator', & 1312 1312 & ' but MEDUSA does not explicitly simulate picophytoplankton' ) … … 1324 1324 ! Surface chlorophyll from HadOCC 1325 1325 zsurfvar(:,:) = HADOCC_CHL(:,:,1) 1326 #elif defined key_medusa && defined key_foam_medusa1326 #elif defined key_medusa 1327 1327 ! Add non-diatom and diatom surface chlorophyll from MEDUSA 1328 1328 zsurfvar(:,:) = trn(:,:,1,jpchn) + trn(:,:,1,jpchd) … … 1340 1340 ! Surface phytoplankton nitrogen from HadOCC multiplied by C:N ratio 1341 1341 zsurfvar(:,:) = trn(:,:,1,jp_had_phy) * c2n_p 1342 #elif defined key_medusa && defined key_foam_medusa1342 #elif defined key_medusa 1343 1343 ! Add non-diatom and diatom surface phytoplankton nitrogen from MEDUSA 1344 1344 ! multiplied by C:N ratio for each … … 1358 1358 CALL ctl_stop( ' Trying to run slphydia observation operator', & 1359 1359 & ' but HadOCC does not explicitly simulate diatoms' ) 1360 #elif defined key_medusa && defined key_foam_medusa1360 #elif defined key_medusa 1361 1361 ! Diatom surface phytoplankton nitrogen from MEDUSA multiplied by C:N ratio 1362 1362 zsurfvar(:,:) = trn(:,:,1,jpphd) * xthetapd … … 1374 1374 CALL ctl_stop( ' Trying to run slphynon observation operator', & 1375 1375 & ' but HadOCC does not explicitly simulate non-diatoms' ) 1376 #elif defined key_medusa && defined key_foam_medusa1376 #elif defined key_medusa 1377 1377 ! Non-diatom surface phytoplankton nitrogen from MEDUSA multiplied by C:N ratio 1378 1378 zsurfvar(:,:) = trn(:,:,1,jpphn) * xthetapn … … 1408 1408 & ' as HADOCC_FCO2(:,:) == HADOCC_FILL_FLT' ) 1409 1409 ENDIF 1410 #elif defined key_medusa && defined key_ foam_medusa && defined key_roam1410 #elif defined key_medusa && defined key_roam 1411 1411 zsurfvar(:,:) = f2_fco2w(:,:) 1412 1412 #elif defined key_fabm … … 1448 1448 & ' as HADOCC_PCO2(:,:) == HADOCC_FILL_FLT' ) 1449 1449 ENDIF 1450 #elif defined key_medusa && defined key_ foam_medusa && defined key_roam1450 #elif defined key_medusa && defined key_roam 1451 1451 zsurfvar(:,:) = f2_pco2w(:,:) 1452 1452 #elif defined key_fabm
Note: See TracChangeset
for help on using the changeset viewer.