- Timestamp:
- 2017-08-14T15:22:09+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r8280 r8436 43 43 USE sbc_oce, ONLY: lk_oasis 44 44 USE oce, ONLY: CO2Flux_out_cpl, DMS_out_cpl, chloro_out_cpl !! Coupling variable 45 #if defined key_foam_medusa 46 USE obs_const, ONLY: obfillflt ! Observation operator fill value 47 #endif 45 48 46 49 IMPLICIT NONE … … 329 332 IF(lwp) WRITE(numout,*) 'Or don t start from uncomplete restart...' 330 333 ENDIF 334 ! 335 # if defined key_foam_medusa 336 !! 2D fields of pCO2 and fCO2 for observation operator on first timestep 337 IF( iom_varid( numrtr, 'PCO2W', ldstop = .FALSE. ) > 0 ) THEN 338 IF(lwp) WRITE(numout,*) ' MEDUSA pCO2 present - reading in ...' 339 CALL iom_get( numrtr, jpdom_autoglo, 'PCO2W', f2_pco2w(:,:) ) 340 CALL iom_get( numrtr, jpdom_autoglo, 'FCO2W', f2_fco2w(:,:) ) 341 ELSE 342 IF(lwp) WRITE(numout,*) ' MEDUSA pCO2 absent - setting to fill ...' 343 f2_pco2w(:,:) = obfillflt 344 f2_fco2w(:,:) = obfillflt 345 ENDIF 346 # endif 347 # endif 348 # if defined key_foam_medusa 349 !! Fields for ocean colour assimilation on first timestep 350 IF( iom_varid( numrtr, 'pgrow_avg', ldstop = .FALSE. ) > 0 ) THEN 351 IF(lwp) WRITE(numout,*) ' MEDUSA pgrow_avg present - reading in ...' 352 CALL iom_get( numrtr, jpdom_autoglo, 'pgrow_avg', pgrow_avg(:,:) ) 353 CALL iom_get( numrtr, jpdom_autoglo, 'ploss_avg', ploss_avg(:,:) ) 354 CALL iom_get( numrtr, jpdom_autoglo, 'phyt_avg', phyt_avg(:,:) ) 355 CALL iom_get( numrtr, jpdom_autoglo, 'mld_max', mld_max(:,:) ) 356 ELSE 357 IF(lwp) WRITE(numout,*) ' MEDUSA pgrow_avg absent - setting to zero ...' 358 pgrow_avg(:,:) = 0.0 359 ploss_avg(:,:) = 0.0 360 phyt_avg(:,:) = 0.0 361 mld_max(:,:) = 0.0 362 ENDIF 331 363 # endif 332 364 … … 498 530 call trc_rst_dia_stat( f2_ccd_arg(:,:),'CCD_ARG') 499 531 !! 532 # endif 533 # if defined key_foam_medusa 534 !! Fields for assimilation and observation operator on first timestep 535 IF(lwp) WRITE(numout,*) ' MEDUSA OBS/ASM fields - writing out ...' 536 # if defined key_roam 537 CALL iom_rstput( kt, nitrst, numrtw, 'PCO2W', f2_pco2w(:,:) ) 538 CALL iom_rstput( kt, nitrst, numrtw, 'FCO2W', f2_fco2w(:,:) ) 539 # endif 540 CALL iom_rstput( kt, nitrst, numrtw, 'pgrow_avg', pgrow_avg(:,:) ) 541 CALL iom_rstput( kt, nitrst, numrtw, 'ploss_avg', ploss_avg(:,:) ) 542 CALL iom_rstput( kt, nitrst, numrtw, 'phyt_avg', phyt_avg(:,:) ) 543 CALL iom_rstput( kt, nitrst, numrtw, 'mld_max', mld_max(:,:) ) 500 544 # endif 501 545 !!
Note: See TracChangeset
for help on using the changeset viewer.