Ignore:
Timestamp:
2019-11-27T16:43:05+01:00 (8 months ago)
Author:
dford
Message:

Get the nitrogen balancing working with 3D chlorophyll increments.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc_3dnitbal/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

    r9292 r11990  
    370370         mld_max(:,:)   = 0.0 
    371371      ENDIF 
     372      IF( iom_varid( numrtr, 'pgrow_avg_3d', ldstop = .FALSE. ) > 0 ) THEN 
     373         IF(lwp) WRITE(numout,*) ' MEDUSA pgrow_avg present - reading in ...' 
     374         CALL iom_get( numrtr, jpdom_autoglo, 'pgrow_avg_3d',  pgrow_avg_3d(:,:,:)  ) 
     375         CALL iom_get( numrtr, jpdom_autoglo, 'ploss_avg_3d',  ploss_avg_3d(:,:,:)  ) 
     376         CALL iom_get( numrtr, jpdom_autoglo, 'phyt_avg_3d',   phyt_avg_3d(:,:,:)   ) 
     377      ELSE 
     378         IF(lwp) WRITE(numout,*) ' MEDUSA pgrow_avg_3d absent - setting to zero ...' 
     379         pgrow_avg_3d(:,:,:) = 0.0 
     380         ploss_avg_3d(:,:,:) = 0.0 
     381         phyt_avg_3d(:,:,:)  = 0.0 
     382      ENDIF 
    372383# endif 
    373384 
     
    555566      CALL iom_rstput( kt, nitrst, numrtw, 'phyt_avg',  phyt_avg(:,:)  ) 
    556567      CALL iom_rstput( kt, nitrst, numrtw, 'mld_max',   mld_max(:,:)   ) 
     568      CALL iom_rstput( kt, nitrst, numrtw, 'pgrow_avg_3d', pgrow_avg_3d(:,:,:) ) 
     569      CALL iom_rstput( kt, nitrst, numrtw, 'ploss_avg_3d', ploss_avg_3d(:,:,:) ) 
     570      CALL iom_rstput( kt, nitrst, numrtw, 'phyt_avg_3d',  phyt_avg_3d(:,:,:)  ) 
    557571# endif 
    558572!! 
Note: See TracChangeset for help on using the changeset viewer.