- Timestamp:
- 2015-09-23T15:26:27+02:00 (9 years ago)
- Location:
- branches/UKMO/dev_r5518_coupling_GSI7_GSI8_landice/NEMOGCM/NEMO/OPA_SRC/SBC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_coupling_GSI7_GSI8_landice/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90
r5708 r5756 157 157 REAL(wp), PUBLIC :: rn_greenland_calving_fraction ! Set fraction of total freshwater flux for iceberg calving - remainder goes to iceshelf melting. 158 158 REAL(wp), PUBLIC :: rn_antarctica_calving_fraction ! Set fraction of total freshwater flux for iceberg calving - remainder goes to iceshelf melting. 159 REAL(wp), PUBLIC :: rn_iceshelf_fluxes_tolerance ! Fractional thresholdfor detecting differences in icesheet masses.159 REAL(wp), PUBLIC :: rn_iceshelf_fluxes_tolerance ! Absolute tolerance for detecting differences in icesheet masses. 160 160 161 161 !! * Substitutions -
branches/UKMO/dev_r5518_coupling_GSI7_GSI8_landice/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r5708 r5756 1280 1280 zgreenland_icesheet_mass_in = zgreenland_icesheet_mass_in / zmask_sum 1281 1281 greenland_icesheet_timelapsed = greenland_icesheet_timelapsed + rdt 1282 IF( ABS( zgreenland_icesheet_mass_in - greenland_icesheet_mass ) /( greenland_icesheet_mass + zepsilon )> zepsilon ) THEN1282 IF( ABS( zgreenland_icesheet_mass_in - greenland_icesheet_mass ) > zepsilon ) THEN 1283 1283 zgreenland_icesheet_mass_b = greenland_icesheet_mass 1284 greenland_icesheet_mass = zgreenland_icesheet_mass_in 1284 1285 ! Only update the mass if it has increased 1286 IF ( (zgreenland_icesheet_mass_in - greenland_icesheet_mass) > 0.0 ) THEN 1287 greenland_icesheet_mass = zgreenland_icesheet_mass_in 1288 ENDIF 1289 1285 1290 IF( zgreenland_icesheet_mass_b /= 0.0 ) & 1286 1291 & greenland_icesheet_mass_rate_of_change = ( greenland_icesheet_mass - zgreenland_icesheet_mass_b ) / greenland_icesheet_timelapsed 1287 1292 greenland_icesheet_timelapsed = 0.0_wp 1288 1293 ENDIF 1289 IF(lwp) WRITE(numout,*) 'Greenland icesheet mass (kg) is ', greenland_icesheet_mass 1294 IF(lwp) WRITE(numout,*) 'Greenland icesheet mass (kg) read in is ', zgreenland_icesheet_mass_in 1295 IF(lwp) WRITE(numout,*) 'Greenland icesheet mass (kg) used is ', greenland_icesheet_mass 1290 1296 IF(lwp) WRITE(numout,*) 'Greenland icesheet mass rate of change (kg/s) is ', greenland_icesheet_mass_rate_of_change 1291 1297 IF(lwp) WRITE(numout,*) 'Greenland icesheet seconds lapsed since last change is ', greenland_icesheet_timelapsed … … 1302 1308 zantarctica_icesheet_mass_in = zantarctica_icesheet_mass_in / zmask_sum 1303 1309 antarctica_icesheet_timelapsed = antarctica_icesheet_timelapsed + rdt 1304 IF( ABS( zantarctica_icesheet_mass_in - antarctica_icesheet_mass ) /( antarctica_icesheet_mass + zepsilon )> zepsilon ) THEN1310 IF( ABS( zantarctica_icesheet_mass_in - antarctica_icesheet_mass ) > zepsilon ) THEN 1305 1311 zantarctica_icesheet_mass_b = antarctica_icesheet_mass 1306 antarctica_icesheet_mass = zantarctica_icesheet_mass_in 1312 1313 ! Only update the mass if it has increased 1314 IF ( (zantarctica_icesheet_mass_in - antarctica_icesheet_mass) > 0.0 ) THEN 1315 antarctica_icesheet_mass = zantarctica_icesheet_mass_in 1316 END IF 1317 1307 1318 IF( zantarctica_icesheet_mass_b /= 0.0 ) & 1308 1319 & antarctica_icesheet_mass_rate_of_change = ( antarctica_icesheet_mass - zantarctica_icesheet_mass_b ) / antarctica_icesheet_timelapsed 1309 1320 antarctica_icesheet_timelapsed = 0.0_wp 1310 1321 ENDIF 1311 IF(lwp) WRITE(numout,*) 'Antarctica icesheet mass (kg) is ', antarctica_icesheet_mass 1322 IF(lwp) WRITE(numout,*) 'Antarctica icesheet mass (kg) read in is ', zantarctica_icesheet_mass_in 1323 IF(lwp) WRITE(numout,*) 'Antarctica icesheet mass (kg) used is ', antarctica_icesheet_mass 1312 1324 IF(lwp) WRITE(numout,*) 'Antarctica icesheet mass rate of change (kg/s) is ', antarctica_icesheet_mass_rate_of_change 1313 1325 IF(lwp) WRITE(numout,*) 'Antarctica icesheet seconds lapsed since last change is ', antarctica_icesheet_timelapsed
Note: See TracChangeset
for help on using the changeset viewer.