2013年11月1日金曜日

Implementing Fe limitation in NPZD model

gotm-4.0.0/src/extras/bio/bio_npzd.F90

!HakaseA
   REALTYPE                  :: lfe=0.1
!HakaseZ

!HakaseA
!                      rzn,rdn,rpdu,rpdl,rzd,aa,g2
                      rzn,rdn,rpdu,rpdl,rzd,aa,g2,lfe
!HakaseZ

!HakaseA
!Co-limitation (N and Light)
!   fnp=rmax*min(par/iopt*exp(1.-par/iopt),n/(alpha+n))*(p+p0)
!Co-limitation (N, Light, and Fe)
   fnp=rmax*min(par/iopt*exp(1.-par/iopt),n/(alpha+n),lfe)*(p+p0)
!HakaseZ

Recompile the model (make command in 'src' directory)

Edit 'bio_npzd.nml' to add the new parameter (lfe)

!HakaseA
!  lfe=         Iron limitation term                              [ - ]
!HakaseZ

!HakaseA
  lfe=         0.1
!HakaseZ

Run the model

Test runs memo

Run 0: No iron limitation (where I left off. run with modified parameters and initial conditions)
Run 1: lfe=0.1 (Minimum value used in Denman & Pena 1999; DP99 hereafter)
Run 2: lfe=0.35 (Value used in Run 7 of DP99)
Run 3: lfe=0.2 (Simply selected)
Run 4: lfe=0.15 (Mean value of lfe in Run 1 and Run3)

0 件のコメント:

コメントを投稿