/* ** svn $Id: estuarybgc_var.h 2232 2012-01-03 18:55:20Z arango $ *************************************************** Hernan G. Arango *** ** Copyright (c) 2002-2012 The ROMS/TOMS Group ** ** Licensed under a MIT/X style license ** ** See License_ROMS.txt ** ************************************************************************ ** ** ** Assigns metadata indices for the Fennel et al. (2006) ecosystem ** ** model variables that are used in input and output NetCDF files. ** ** The metadata information is read from "varinfo.dat". ** ** ** ** This file is included in file "mod_ncparam.F", routine ** ** "initialize_ncparm". ** ** ** ************************************************************************ */ /* ** Model state biological tracers. */ CASE ('idTvar(iNO3_)') idTvar(iNO3_)=varid CASE ('idTvar(iNH4_)') idTvar(iNH4_)=varid CASE ('idTvar(iPhyt)') idTvar(iPhyt)=varid CASE ('idTvar(iZoop)') idTvar(iZoop)=varid CASE ('idTvar(iLDeN)') idTvar(iLDeN)=varid CASE ('idTvar(iSDeN)') idTvar(iSDeN)=varid CASE ('idTvar(iChlo)') idTvar(iChlo)=varid # ifdef SPECTRAL_LIGHT /* ** Photosynthetically Available Radiation (PAR).. */ CASE ('idPARo') idPARo=varid CASE ('idPARs') idPARs=varid CASE ('idSpKd') idSpKd=varid # ifdef CDOM_VARIABLE CASE ('idTvar(iCDMC)') load=.FALSE. varid=varid-1 DO i=1,Ndom varid=varid+1 idTvar(iCDMC(i))=varid DO ng=1,Ngrids Fscale(varid,ng)=scale Iinfo(1,varid,ng)=gtype END DO WRITE (Vname(1,varid),'(a,i1)') 'CDM_C', i WRITE (Vname(2,varid),'(a,i1)') & & 'color degradational matter, carbon group ', i WRITE (Vname(3,varid),'(a)') & & TRIM(ADJUSTL(Vinfo(3))) WRITE (Vname(4,varid),'(a,a)') & & TRIM(Vname(1,varid)), ', scalar, series' WRITE (Vname(5,varid),'(a)') & & TRIM(ADJUSTL(Vinfo(5))) END DO varid=varid+1 # endif # endif # ifdef SAV_BIOMASS CASE ('iddins') iddins=varid CASE ('iddinw') iddinw=varid CASE ('iddowc') iddowc=varid CASE ('idwsvl') idwsvl=varid CASE ('idsagb') idsagb=varid CASE ('idsbgb') idsbgb=varid CASE ('idsvpp') idsvpp=varid CASE ('idsvam') idsvam=varid CASE ('idsgar') idsgar=varid CASE ('idsvbr') idsvbr=varid CASE ('idsvrs') idsvrs=varid CASE ('idsvbg') idsvbg=varid CASE ('idsvag') idsvag=varid CASE ('idsbgr') idsbgr=varid CASE ('idsbgm') idsbgm=varid # endif # ifdef CARBON CASE ('idTvar(iTIC_)') idTvar(iTIC_)=varid CASE ('idTvar(iTAlk)') idTvar(iTAlk)=varid CASE ('idTvar(iLDeC)') idTvar(iLDeC)=varid CASE ('idTvar(iSDeC)') idTvar(iSDeC)=varid # endif # ifdef OXYGEN CASE ('idTvar(iOxyg)') idTvar(iOxyg)=varid # endif /* ** Adjoint sensitivity state biological tracers. */ #if defined AD_SENSITIVITY || defined IS4DVAR_SENSITIVITY || \ defined OPT_OBSERVATIONS || defined SENSITIVITY_4DVAR || \ defined SO_SEMI CASE ('idTads(iNO3_)') idTads(iNO3_)=varid CASE ('idTads(iNH4_)') idTads(iNH4_)=varid CASE ('idTads(iPhyt)') idTads(iPhyt)=varid CASE ('idTads(iZoop)') idTads(iZoop)=varid CASE ('idTads(iLDeN)') idTads(iLDeN)=varid CASE ('idTads(iSDeN)') idTads(iSDeN)=varid CASE ('idTads(iChlo)') idTads(iChlo)=varid # ifdef CARBON CASE ('idTads(iTIC_)') idTads(iTIC_)=varid CASE ('idTads(iTAlk)') idTads(iTAlk)=varid CASE ('idTads(iLDeC)') idTads(iLDeC)=varid CASE ('idTads(iSDeC)') idTads(iSDeC)=varid # endif # ifdef OXYGEN CASE ('idTads(iOxyg)') idTads(iOxyg)=varid # endif #endif /* ** Biological tracers open boundary conditions. */ CASE ('idTbry(iwest,iNO3_)') idTbry(iwest,iNO3_)=varid CASE ('idTbry(ieast,iNO3_)') idTbry(ieast,iNO3_)=varid CASE ('idTbry(isouth,iNO3_)') idTbry(isouth,iNO3_)=varid CASE ('idTbry(inorth,iNO3_)') idTbry(inorth,iNO3_)=varid CASE ('idTbry(iwest,iNH4_)') idTbry(iwest,iNH4_)=varid CASE ('idTbry(ieast,iNH4_)') idTbry(ieast,iNH4_)=varid CASE ('idTbry(isouth,iNH4_)') idTbry(isouth,iNH4_)=varid CASE ('idTbry(inorth,iNH4_)') idTbry(inorth,iNH4_)=varid CASE ('idTbry(iwest,iPhyt)') idTbry(iwest,iPhyt)=varid CASE ('idTbry(ieast,iPhyt)') idTbry(ieast,iPhyt)=varid CASE ('idTbry(isouth,iPhyt)') idTbry(isouth,iPhyt)=varid CASE ('idTbry(inorth,iPhyt)') idTbry(inorth,iPhyt)=varid CASE ('idTbry(iwest,iZoop)') idTbry(iwest,iZoop)=varid CASE ('idTbry(ieast,iZoop)') idTbry(ieast,iZoop)=varid CASE ('idTbry(isouth,iZoop)') idTbry(isouth,iZoop)=varid CASE ('idTbry(inorth,iZoop)') idTbry(inorth,iZoop)=varid CASE ('idTbry(iwest,iSDeN)') idTbry(iwest,iSDeN)=varid CASE ('idTbry(ieast,iSDeN)') idTbry(ieast,iSDeN)=varid CASE ('idTbry(isouth,iSDeN)') idTbry(isouth,iSDeN)=varid CASE ('idTbry(inorth,iSDeN)') idTbry(inorth,iSDeN)=varid CASE ('idTbry(iwest,iLDeN)') idTbry(iwest,iLDeN)=varid CASE ('idTbry(ieast,iLDeN)') idTbry(ieast,iLDeN)=varid CASE ('idTbry(isouth,iLDeN)') idTbry(isouth,iLDeN)=varid CASE ('idTbry(inorth,iLDeN)') idTbry(inorth,iLDeN)=varid CASE ('idTbry(iwest,iChlo)') idTbry(iwest,iChlo)=varid CASE ('idTbry(ieast,iChlo)') idTbry(ieast,iChlo)=varid CASE ('idTbry(isouth,iChlo)') idTbry(isouth,iChlo)=varid CASE ('idTbry(inorth,iChlo)') idTbry(inorth,iChlo)=varid #ifdef CARBON CASE ('idTbry(iwest,iSDeC)') idTbry(iwest,iSDeC)=varid CASE ('idTbry(ieast,iSDeC)') idTbry(ieast,iSDeC)=varid CASE ('idTbry(isouth,iSDeC)') idTbry(isouth,iSDeC)=varid CASE ('idTbry(inorth,iSDeC)') idTbry(inorth,iSDeC)=varid CASE ('idTbry(iwest,iLDeC)') idTbry(iwest,iLDeC)=varid CASE ('idTbry(ieast,iLDeC)') idTbry(ieast,iLDeC)=varid CASE ('idTbry(isouth,iLDeC)') idTbry(isouth,iLDeC)=varid CASE ('idTbry(inorth,iLDeC)') idTbry(inorth,iLDeC)=varid CASE ('idTbry(iwest,iTIC_)') idTbry(iwest,iTIC_)=varid CASE ('idTbry(ieast,iTIC_)') idTbry(ieast,iTIC_)=varid CASE ('idTbry(isouth,iTIC_)') idTbry(isouth,iTIC_)=varid CASE ('idTbry(inorth,iTIC_)') idTbry(inorth,iTIC_)=varid CASE ('idTbry(iwest,iTAlk)') idTbry(iwest,iTAlk)=varid CASE ('idTbry(ieast,iTAlk)') idTbry(ieast,iTAlk)=varid CASE ('idTbry(isouth,iTAlk)') idTbry(isouth,iTAlk)=varid CASE ('idTbry(inorth,iTAlk)') idTbry(inorth,iTAlk)=varid #endif #ifdef OXYGEN CASE ('idTbry(iwest,iOxyg)') idTbry(iwest,iOxyg)=varid CASE ('idTbry(ieast,iOxyg)') idTbry(ieast,iOxyg)=varid CASE ('idTbry(isouth,iOxyg)') idTbry(isouth,iOxyg)=varid CASE ('idTbry(inorth,iOxyg)') idTbry(inorth,iOxyg)=varid #endif /* ** Biological tracers point Source/Sinks (river runoff). */ CASE ('idRtrc(iNO3_)') idRtrc(iNO3_)=varid CASE ('idRtrc(iNH4_)') idRtrc(iNH4_)=varid CASE ('idRtrc(iPhyt)') idRtrc(iPhyt)=varid CASE ('idRtrc(iZoop)') idRtrc(iZoop)=varid CASE ('idRtrc(iLDeN)') idRtrc(iLDeN)=varid CASE ('idRtrc(iSDeN)') idRtrc(iSDeN)=varid CASE ('idRtrc(iChlo)') idRtrc(iChlo)=varid #ifdef CARBON CASE ('idRtrc(iTIC_)') idRtrc(iTIC_)=varid CASE ('idRtrc(iTAlk)') idRtrc(iTAlk)=varid CASE ('idRtrc(iLDeC)') idRtrc(iLDeC)=varid CASE ('idRtrc(iSDeC)') idRtrc(iSDeC)=varid #endif #ifdef OXYGEN CASE ('idRtrc(iOxyg)') idRtrc(iOxyg)=varid #endif #ifdef DIAGNOSTICS_BIO /* ** Biological tracers term diagnostics. */ # ifdef DENITRIFICATION CASE ('iDbio2(iDNIT)') iDbio2(iDNIT)=varid # endif # ifdef CARBON CASE ('iDbio2(iCOfx)') iDbio2(iCOfx)=varid CASE ('iDbio2(ipCO2)') iDbio2(ipCO2)=varid # endif # ifdef OXYGEN CASE ('iDbio2(iO2fx)') iDbio2(iO2fx)=varid # endif CASE ('iDbio3(iPPro)') iDbio3(iPPro)=varid CASE ('iDbio3(iNO3u)') iDbio3(iNO3u)=varid #endif