Commit 91ed66d4 authored by ph290's avatar ph290
Browse files

making sure nutrients were properly picked up from WOA

parent 21521a38
......@@ -485,7 +485,7 @@ END MODULE GRAPHICS_VARIABLES
read(1,'(i1)') imode
do i = 1,iline
read(1,'(16x,11f6.1)') smaj1, smin1, smaj2, smin2, smaj3, smin3, smaj4, smin4, smaj5, smin5
read(1,'(16x,10f6.1)') smaj1, smin1, smaj2, smin2, smaj3, smin3, smaj4, smin4, smaj5, smin5
!Note, reading in 1st 6 constituents of teh tides, but only using the first 5
! write(6,'(16x,6f6.1)') smaj1, smin1, smaj2, smin2, smaj3, smin3
if(i.eq.iline) then
......@@ -883,7 +883,7 @@ read(1,'(i1)') imode
do i = 2,iline+1
read(1,'(f8.3,f8.3,62x,f6.1)') lon,lat,alldepth
! read(1,'(f8.3,f8.3,38x,f6.1)') lon,lat,alldepth
!write(6,'(f8.3,f8.3,38x,f6.1)') lon,lat,alldepth
! write(6,'(f8.3,f8.3,38x,f6.1)') lon,lat,alldepth
if(i.eq.iline+1) newdepth = alldepth
enddo
close (1)
......@@ -911,8 +911,8 @@ endif
!open(1,file='s12_m2_s2_n2_h.dat',status='old')
!open(1,file='s12_m2_s2_n2_h_sec.dat',status='old')
!open(1,file='s12_m2_s2_n2_h_tim.dat',status='old')
vismax=0.1; Nz_bg=1.0e-5; Kz_bg=1.0e-5; first_temp=5.0; lambda=0.1; heat_shade=0.012; bed_din=7.0; din_rate=10.0
!vismax=0.1; Nz_bg=1.0e-5; Kz_bg=1.0e-5; first_temp=5.0; lambda=0.1; heat_shade=0.012; bed_din=woa_nutrient; din_rate=10.0
vismax=0.1; Nz_bg=1.0e-5; Kz_bg=1.0e-5; first_temp=5.0; lambda=0.1; heat_shade=0.012; bed_din=woa_nutrient; din_rate=10.0
par_atten=0.1; par_percent=0.45
uamp(1)=0.4; uamp(2:5)=0.0; uphase(1:5)=0.0; vphase(1:5)=0.0
vamp(1:5)=0.0
......@@ -1262,7 +1262,7 @@ vely_new(1:200)=0.0; temp_new(1:200)=first_temp; Nz(1:200)=vismax/2.0; Kz(1:200)
!
! ....and the biological arrays...
!
x_old(1:200)=0.1; ni_old(1:200)=0.1; ni_old(1:200)=0.1; s_old(1:200)=woa_nutrient; x_new(1:200)=0.1
x_old(1:200)=0.1; ni_old(1:200)=0.1; ni_old(1:200)=0.1; s_old(1:200)=bed_din; x_new(1:200)=0.1
ni_new(1:200)=0.1; x_new(1:200)=0.1
x_old(1:200)=x_new(1:200); ni_new(1:200)=0.1; s_new(1:200)=bed_din
tpn1=0.0; tpn2=0.0; tpg1=0.0; tpg2=0.0; total_flux=0.0
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment