͜STL2STM2BookC_StringPrint_SetupTPrintStartupSound_TabLISTBook_PrefsAct List Scene Drawing_Model Draw_Specs diagram_struct_arrayDependentsScriptPad_ListManuscriptPlayPlay_PrefsModelSubscript_Set_ListSubscript_SetArraySimulation_ModelRun_Specsentity_struct_arrayentity_structint16_arraydouble_arraytoken_type_arraytoken_type run_handle_array!Draw_Index"Poster_Info#queue$Units%diagram_struct&Point_array'movie_attach(Picture_Tab)Pad*Model_Index_array+Picture_Data,Picture-AliasRecord.Pathname/graph_type0curve_type_array1curve_type2Model_Index3Panel4Sector_Specs5SET6Sense_Setup7select_struct_array8allow_struct_array9Picture_array:Section_List;Font_Map<Font_Map_Record_Array=Font_Map_Record>Window_List?Model_Window@Pad_WindowA@@c / / 0 / /p /xxHHدG(HHK'`>>  ,, / 0 0      / 0    JJ  drmd   caNc`@@c_ / / # /wwhPP@8?2   VVc / # c / #X** #  # &1 Dim_Name_1 P #|Shoots Roots_RhyWrack EpiphytesPhytoplankton FoodSourceb'jj@?=Time'2c 7??????????????convert_wspd1 !BEHM: Wind speed is probably in units of knots. To convert from knots to m/hr: 1 knot = 1 nautical mile/hr = 1.85 km/hr = 1850 m/hr" #&$%c ????????????+*^???Alpha 1+Ua_Ub^2+2*Ua_Ub*COS(Angle)NB4  N     N       !" #&$%c @3?@7@7@3@3@3????@3Hours MOD(TIME,24)   , !" #&$%s??@"@v????????Days0}tgUJ !" #&$%f  ????????? Count_hoursIF Hours = 0 THEN 1*DT ELSE 00          !" #&$%c ?Q?Q???Q?Q?Q?Q????Q?x?Q Alpha_PhBio0.43?   ? Q! was 0.43" #&$%c2 @aNJ?@aNJ??@c1>@^ @e8U?pR@aNJ?@aNJ?@aNJ?@aNJ?@j\$L@aNJ??@aNJ?Angle3 (ABS(Current_direction-Wind_direction)/8)*(2*PI) 0            !" #&$%c ?bM񖩟?bM񖩟??bM񖩟?bM񖩟?bM񖩟?bM񖩟?bM񖩟?bM񖩟?bM񖩟?bM񖩟????bM񖩟? ?bM񖩟ic_TSS0.0083   ?bM 񖩟!" #&$%4ci#@h^TM!@h?@FS??@w |@g@h^TM!@h?@FS@h^TM!@h^TM!{???@h?@FS?-@h^TM!Solar_radiationj0.7*(916.73*(h_angle*SIN(Latitude*0.0174533)*SDCLN+SIN(h_angle) *COS(Latitude*0.0174533)*COS(Decline)))bg%##  ?ff ffff ( @Դ =p b  * ?ΧG "@_     g  b    * ?ΧG "@_            !|Two versions of this were written in the model: 677.5-371.5*((2*PI*(DAY+10))/365) does not work next is the QEXT parameter." #&$%c@L@L?@L0}4@,#?@L0}4@L @L@L@L@L???@L? @L SF_watercell_size*water_depth 32^X~ixl     !" #&$%f C>vn>M(g??>vn>M(g>s9(>uv^???>vn >uv^egested_matter cons_egestIC I!" #&$%c o??????j  River_NACL0 !" #&$%c ??`bM񖩟?`bM񖩟?`bM񖩟?`bM񖩟*??????  ic_wrack06 !+{kg/m2, adjust this value for IC changes}" #&$%c  @@@@@@@@@@Aj@? @ cell_size10000"+67| q{1<}ts0 '!{M^2} =100m * 100m" #&$%c ? ? ??"j?E12*A?`YH?0? ? ? ? ???? ?? Current_velocitysqrt(EW_vector^2+NS_vector^2)7N$    7     !" #&$%c9J @eZ@Qςz??@4[x?<>@eZ@Qςz@Qςz@QςzS???@aUk?-@QςzPAR_Epi:PAR_Surface*exp(-k_H2O*max(0,(water_depth-leaf_length)))%4 %            !TBEHM: Calculates light for epiphytes: calculated for top of leaf in water column " #&$%c a?At D8Imc?????wrack_X_N@S:00 !" #&$%>c N-@!o?e3??@#1Adפ?g/L1@!o?e3@!o@!o???@;~wt?-@!o Cloud_coverIF Rand < few_clouds THEN 3*(Rand/few_clouds) ELSE IF Rand > some_clouds THEN 6.6+3*((Rand-some_clouds)/(100-some_clouds)) ELSE 3.3+3*((Rand-few_clouds)/(some_clouds-few_clouds))t%--  t      t   -       : @ff ffff     d      -  # @ ff ffff   t   t     !" #&$%cS ? d"Kmj?W,???cēcW? d"Kmj?W,?W,?W, Kn? d"Kmj??W,Loc_SAT -TAN(Latitude*PI/180)*TAN(max_Sun_angle) {tan (lat in radians) x tan declination}4            ! Trust it" #&$%!cA@"7@"7?@0r7d ,Fj>@,Fj<0r7d @"7@"7@"7@"7/???@"7?@"7EW_wind_vectorB(Wind_speed*convert_wspd)*SIN(((180+Wind_direction)*2*PI)/360)""@          h     !0BEHM: convert_wspd converts from knots to m/hr" #&$%c {@ h(/@ `^ ??@/@ =ڛD@ h(/@ `^ @ h(/@ h(/???@ `^ ?-@ h(/ Photo_period7.639437*h_angleb  @ h b  !" #&$%!cA @uPr@uPr?@;#.Zd@)(#g1qp@uPr@uPr@uPr@uPr/???@uPr?@uPrNS_wind_vectorB(Wind_speed*convert_wspd)*Cos(((Wind_direction+180)*2*PI)/360)n@          h     !0BEHM: convert_wspd converts from knots to m/hr" #&$%"c]K?J*?J*?CE-/?.2p3 P*?J*?J*?J*?J*k????J*??J* E_vector^if E_Current_Vector<0 then 0 else (E_Current_Vector*convert_currents)+(0.035*EW_wind_vector)a"D        a  ( ? Q    !lWind drift factor is 3.5%. McDonald 1992 oilspill model. BEHM: convert_currents converts from m/s to m/hr" #&$%f>">! v?m33]?KNOGAj*P?KNO>">! v>"d7>">! v>"cons_mortalityConsumers*(C_mort_rt*DT)++ +     !" #&$%c?zG{?zG{???zG{?zG{?zG{?zG{????zG{?-?zG{ k_water0.04   ?z G{!Short: .00007 is base light extinction coeff for water w/o particulates Behm: 0.04 m^-1 is PAR extinction coefficient for water (Buzzelli et al. 1999) " #&$%"c]n?携˥q?携˥q??-F k6a?e=kD6T?携˥q?携˥q?携˥q?携˥q1????携˥q??携˥q N_vector^if N_Current_Vector<0 then 0 else (N_Current_Vector*convert_currents)+(0.035*NS_wind_vector)a7D        a  ( ? Q    !2BEHM: convert_currents converts from m/s to m/hr" #&$%c? ?ۍ05????? ?ۍ05?ۍ05?ۍ05???? ?-?ۍ05k_H2Ok_water+(k_Pht*Pht_conc)l3   l   !.0007 is coef of absorption of the combined suspended solids .00007 is base light extinction coeff for water w/o particulates Behm: Based on Buzzelli et al. 1999" #&$%g-??????$Light!" #&$%cF ?$ʹz?$ʹz??^҃V? ,C?+׼!?$ʹz?$ʹz?$ʹz?$ʹz????$ʹz??$ʹzN_wrack_vectorGIF N_vector<0 THEN 0 ELSE N_vector+(Wind_Wrack_factor*NS_wind_vector)(t4            !" #&$%c11?+jg?+jg?@\zG@\ }-@\zG@[bñ ?+jg?+jg?+jg?+jgC????+jg? ?+jgSettling_velocity2MIN((settling_rt),(settling_rt)/Eddy_viscosity)&h & &    !D((2.6-1)*9.8/18*0.001)*0.0039^2 = 1.325 cm per sec (11.448 meters per day) is the maximum settling velocity of silts in non tubulent water.(2.6 kg/l is density of minerals, 1kg/l is density of fluid (water).9.8 is gravity; 0.001 poise (kgm^-1s^-1) is dynamic viscosity of clear water (normally does change with Temperature and Clay concentration, this is not in the model) m0.0039 is the diameter of silt particles (mm). Reworking the formular with excluding the viscosity gives this particular equation for silt. The eddy viscosity is related to the shearstress in the water" #&$%*f?@????Wrack_OM_decompIF on_map=3 THEN wrack ELSE IF on_map = 2 THEN wrack*decomp_temp_Air*Shredding_by_cons ELSE wrack*decomp_temp_H2O*Shredding_by_cons 65&p6d      6        6 5  &    6 p  &  !AThe decomp temp coeff rises from near 0 to 1.0, signmoid, with temp increasing from 0 to 40C. Because this is above sed organic material, aerobic decomp is assumed. Decomp increases with the moisture proportion, remains at max while in water column. kg*dimensionless(0-1)*dimensionless(0-1)*some function of shredders" #&$%c  ?zG{@Ym?zG{?zG{??H?zG{?b7?zG{?zG{?zG{?zG{?zG{/????zG{? ?zG{Eddy_viscosity Shear_stressBXX@$@4@>@D@I@N@Q@T@V@Y?zG{@1@@p@E (@J@S@Oy=c@R*0@Tu@V@W_o @Y B!0Larry talk to Larry Sandford (student of Grand)" #&$%c T ?CK$?CK$??ݑ1;?͐c@?pU?CK$?CK$?CK$?CK$?????^?CK$? ?CK$ DW_Wave_L(9.8*(Wave_period^2))/(2*PI)c0 ( @#         !" #&$%c??????y!River_P0 !" #&$%cD @&@gEב@&@gEב? rʣF;?j6CÃZH@&@gEב@&@gEב@&@gEב@&@gEב???@&@gEב?"@&@gEבE_wrack_vectorEIF E_vector<0 THEN 0 ELSE E_vector+Wind_Wrack_factor*EW_wind_vector(s4            !" #&$%c {?@r@r@r@r????@ ?# OMtotBio@S0Y !This variable repesents the aggregate amount of biomass in the four adjacent cells that could provide mobile consumers with a food source." #&$%c 3?????????????x$?Alpha_NPhBio1 !what is this ? RB .5" #&$%c ' @g,k @S'G??@hW@P`Bf@g,k @S'G@S'G@S'G???@c?Wg?-%@S'G PAR_Surface(Solar_radiation*(1-(0.071*Cloud_cover))3$   ( ?- V     !" #&$%c????????????&?Shredding_by_cons1 !" #&$%!cL!>-z_>Nb34??%U?o5z!H*RBEl>hsÆ>-z_>Nb34>wSc@>wSc@???>-z_?'>wSc@Cons_ingest_EpiMif Epiphytes<1 then 0 else Min(Epiphytes,(EP_pref/pref_tot)*Con_pot_ingest)-sU@            -   !" #&$% cQ@*OIt@*OIt?@wV@0^@ab@@*OIt@*OIt@*OIt@*OIt???@*OIt?(@*OItVector_Abs_SumRABS(E_wrack_vector)+ABS(N_wrack_vector)+ ABS(S_wrack_vector)+ABS(W_wrack_vector)""ts< "           "    !" #&$%c )@@? 6C-@@? 6C-@@@@2@?y)@DIN_sed_conc4.2   @ !3kilograms per M3 the same concentration everywhere" #&$%OcT >?pxl9R$ɺ??KaWtMD?. 6eWtMD?pxl9R$ɺd\d\???ͅ_?*d\TideU0.6096*(-0.3+0.14 *COS(( 15.041 *TIME -75.1 )*2*PI/360)+ 0.065 *COS(( 30.082 *TIME -304.3 )*2*PI/360)+ 0.041 *COS(( 29.528 *TIME -298.9 )*2*PI/360)+ 1.468 *COS(( 28.984 *TIME -270.7 )*2*PI/360)+ 0.1 *COS(( 57.968 *TIME -16.6 )*2*PI/360)+ 0.164 *COS(( 86.952 *TIME -3.7 )*2*PI/360)+ 0.32 *COS(( 28.439 *TIME -227.7 )*2*PI/360)+ 0.043 *COS(( 27.895 *TIME -199.7 )*2*PI/360)+ 0.047 *COS(( 13.943 *TIME -51.2 )*2*PI/360)+ 0.047 *COS(( 14.958 *TIME -70.7 )*2*PI/360)+ 0.239 *COS(( 30 *TIME -310.6 )*2*PI/360)+ 0.014 *COS(( 29.958 *TIME -283.1 )*2*PI/360)+ 0.062 *COS(( 28.512 *TIME -234.5 )*2*PI/360))>>  ? - ?33 3333  + ?瑷 Q ( @. E   - @Rf ffff      h      + ? =p ( @> E   - @s      h      + ? E ( @=+  I   - @rf ffff      h      + ?| hr ( @< lC   - @p3 3333      h      + ? ( @L lC   - @0      h      + ? E ( @U hr   - @      h      + ?z G{ ( @߼ >#]?UC?y- DIN_concDIN_Maps*DIN_scen*River_N,EM ,     !Ddin maps * ratio of current total river input/max total river input" #&$%c ?PbM񖩟?PbM񖩟???PbM񖩟?PbM񖩟?PbM񖩟?PbM񖩟????PbM񖩟? .?PbM񖩟 River_CF0.001   ?PbM 񖩟!" #&$%c0??F????y/ PO4_Maps0 !GThese are the hydrology maps created using 1 minute timesteps (kg/m3)." #&$%c=?PbM񖩟?PbM񖩟???PbM񖩟?PbM񖩟?PbM񖩟?PbM񖩟????PbM񖩟?0?PbM񖩟ic_Pht0.0010  * ?PbM 񖩟!" #&$%c ?@r@r@r@r????@ ?1 OMtotBio@W0V !This variable repesents the aggregate amount of biomass in the four adjacent cells that could provide mobile consumers with a food source." #&$%c @4@4??@4@4@4@4???@4?2@4half_sat_H2O20p !" #&$%c</ ????q????-3PAR_Sw=PAR_Surface*exp(-k_H2O*max(0,(water_depth-leaf_length_Sw)))%;4 %    ;        !rBEHM: Calculates light for seaweed: calculated for top of leaf (assuming leaf length is 0.5 m) in water column " #&$%c %>uMUi>uMUi??`bM񖩟?`bM񖩟?`bM񖩟?`bM񖩟>uMUi>uMUi>uMUi>uMUiN????>uMUi?4>uMUi ic_consumer&(ic_mac_PhBio/sept_to_jan_conv)*.018p+ 8 p  * ?z G{  !OFor every 100 grams of eelgrass there is 1 gram of consumers. Guessed {kg/m2}" #&$%f '?@@@@???5 Cons_in_X(ConstoE@W+ConstoN@S+ConstoS@N+ConstoW@EPFLO+ P F  L  O  !" #&$%cL?Q?Q???Q?Q?Q?Q?Q?QC????Q?x6?Q ic_NPHB.015   ?Q !DInitial Non-photosynthic Biomass in Kg =1/3 PHB Short 1980 .03 .02" #&$%c j>u?5(v>uMUh??v?`bM񖩟GT`?`bM񖩟>u?5(v>uMUh>u?5(v>u?5(v????>uMUh?7>u?5(v Consdens!(River_Con+Consumers)/cell_size+ QSXZ  +    ! {kg/m2}" #&$%c W?Q?Q???Q?Q?Q?Q?Q?Q.????Q?x8?Qic_mac_PhBio.06 4<P   ?Q !/Initial Photosynthic Biomass in Kg/m2 .09 .06" #&$%c ?>Ëh>Ëh>Ëh>Ëh?????9 Consdens@N0Q !" #&$%c ?>Ëh>Ëh>Ëh>Ëh?????: Consdens@S0X !" #&$%c??????????????-;?leaf_length_Sw0.53   ? !" #&$%cB ??PM3?V#K6?????j> Salt_X_NCIF sf_wt_X_N >0 THEN Salt_conc*sf_wt_X_N*convert_currents ELSE 0a4        a    !" #&$%c ??j? Salt_X_N@S0 !" #&$%g ??????A$TSS!" #&$%c / :pj:pj?>u47:t?? /%/45 b:pj:pj:pj:pj?Ehe:pj?B:pj Shear_stress0 .5*fric_coef*Alpha*SQRT((orbital_velocity)^2) Tq4  ? T    q      !{Grant and Madsen}" #&$%c ?>Ëh>Ëh>Ëh>Ëh?????C Consdens@E0Z !" #&$%c ?>Ëh>Ëh>Ëh>Ëh?????D Consdens@W0S !" #&$%c0'>ս݁>`D?>J]5=uZ@n )>ս݁>`D>`hi7>`hi7???>ս݁?E>`hi7Cons_ingest_SW1Min(Seaweeds,(SW_pref/pref_tot)*Con_pot_ingest)|-dU  |    -   !1/h" #&$%c ?????????F ConstoN@S05 !" #&$%f ?@ A-m? =pGhj? =p???G Cons_out_X ConstoE+ConstoN+ConstoS+ConstoWK\JH+ K \  J  H  !" #&$% cR?? A-m? =pGhj? =p????HConstoWSif (Density_X_W+Food_X_W) > 1 then Consumers else (Density_X_W+Food_X_W)*ConsumersSV+G< S V      +   S V  +  !" #&$%f I>vn>M(g?*M\?g6RdwGRBEy^>&>vn>M(g>s9(>uv^ >vnI>uv^ cons_egest Cons_ingest*(C_egest_eff)U +  U   ! units kg/h" #&$% cR?? A-m? =pGhj? =p????JConstoSSif (Density_X_S+Food_X_S) > 1 then Consumers else (Density_X_S+Food_X_S)*ConsumersXY+G< X Y      +   X Y  +  !" #&$% cR?? A-m? =pGhj? =p????KConstoESif (Density_X_E+Food_X_E) > 1 then Consumers else (Density_X_E+Food_X_E)*ConsumersZ[+G< Z [      +   Z [  +  !" #&$%c ?????????L ConstoS@N05 !" #&$%c FR ??D҉^?[?-ͨn?M Epi_Nut_limGmin( DIN_conc/(DIN_conc+Epi_Sat_N), PO4_conc/(PO4_conc+Epi_Sat_P) ) -0 - -          !oMich-Menton, assumes only one limiting nutrient. half sat coeffs are broad approximations (2.5E-5g/l=25 g/l) " #&$%cU ??f6=LN"?NUa_UbV If orbital_velocity > Current_velocity THEN Current_velocity/orbital_velocity ELSE 0qT, q      q    !" #&$%c ?????????O ConstoW@E05 !" #&$%c ?????????P ConstoE@W05 !" #&$%$c e??ΐ3ZZ?xE?????Q Density_X_NfIF Consdens < Consdens@N OR Consdens <1 THEN 0 ELSE ((Travel_time)*(Consdens-Consdens@N))/Consdens 798\L 7 9  7        8 7 9   7  !" #&$%$cf?????R Food_X_NgIF OMtotBio > OMtotBio@N OR OMtotBio <1 THEN 0 ELSE ((Travel_time)*(OMtotBio@N-OMtotBio))/OMtotBio >n8\L > n  >        8 n >   >  !" #&$%$c d??ΐ3ZZ?xE?????S Density_X_WeIF Consdens < Consdens@W OR Consdens <1 THEN 0 ELSE ((Travel_time)*(Consdens-Consdens@W))/Consdens7D8HL 7 D  7        8 7 D   7  !" #&$%c  ?tzG{?????????????g??T? fric_coefUa_UbNB?zG{?zG{?Q?zG{??Q?Q?zG{? =p ???7K?úYJ?&J?8YJ?MjOw?u!.?|[X?x@?vR<64?toiE N!" #&$%f dU >%Nk>3pq?*M\?og6RdwH*RBEy^? &>%Nk>3pq>L܁X>NN" >%NkU>NN" Cons_ingeste(Cons_ingest_Epi+Cons_ingest_NPhBio+Cons_ingest_PhBio+Cons_ingest_Pht+Cons_ingest_SW+Cons_ingest_Wr) 'LKdE)I+, ' L  K  d  E  )  ! units kg/h" #&$%$ce?????V Food_X_WfIF OMtotBio > OMtotBio@W OR OMtotBio <1 THEN 0 ELSE ((Travel_time)*(OMtotBio@W-OMtotBio))/OMtotBio>18HL > 1  >        8 1 >   >  !" #&$%c c@4@4??@4@4@4@4@4@4???@4?xW@4L_to_R_limit20{~ !" #&$%$c d??ΐ3ZZ?xE?????X Density_X_SeIF Consdens < Consdens@S OR Consdens <1 THEN 0 ELSE ((Travel_time)*(Consdens-Consdens@S))/Consdens7:8JL 7 :  7        8 7 :   7  !" #&$%$ce?????Y Food_X_SfIF OMtotBio > OMtotBio@S OR OMtotBio <1 THEN 0 ELSE ((Travel_time)*(OMtotBio@S-OMtotBio))/OMtotBio>#8JL > #  >        8 # >   >  !" #&$%$c d??ΐ3ZZ?xE?N????Z Density_X_EeIF Consdens < Consdens@E OR Consdens <1 THEN 0 ELSE ((Travel_time)*(Consdens-Consdens@E))/Consdens7C8KL 7 C  7        8 7 C   7  !ODispension of consumers over the area while travelling across cell boundaries." #&$%$ce?????[ Food_X_EfIF OMtotBio > OMtotBio@E OR OMtotBio <1 THEN 0 ELSE ((Travel_time)*(OMtotBio@E-OMtotBio))/OMtotBio>+8KL > +  >        8 + >   >  !" #&$% cRB?? A-m? =pGhj? =p????\ConstoNSif (Density_X_N+Food_X_N) > 1 then Consumers else (Density_X_N+Food_X_N)*ConsumersQR+G< Q R      +   Q R  +  !" #&$%c ] ?68tJ?68tJ?$Dڌ?_;+@ z+$ ?A3?68tJ?68tJ?68tJ?68tJ?`5 ~=?68tJ?]?68tJg_d! (water_depth*9.7)/Wind_speed^2_$  * @#ff ffff      !" #&$%c^ @%y"mI@%y"mI@DPq>+@cϷ@zyO'B3?5c%@%y"mI@%y"mI@%y"mI@%y"mI@1JV@%y"mI?^@%y"mIg_f (Fetch*9.7)/Wind_speed^2;`$ ; * @#ff ffff      !" #&$%c_ ?&yq?&yq?W?_?0P?gd>Ը=?&yq?&yq?&yq?&yq?{HI?&yq?_?&yqH1 .53*(g_d^.75)]a$  ? \( ] ^ ?   !" #&$%c`?Ň?Ň?nn?CȮD?%?~qHy ?Ň?Ň?Ň?Ň?e?Ň?`?ŇH2 .00565*(SQRT(g_f))^b  ?w$t SڗM ^    !" #&$%"c$a?~&?~&?`͘7@?v6b?ؾ,C?݄@ld?~&?~&?~&?~&?ZN d?~&?a?~&Ha% (EXP(H1)-EXP(-H1))/EXP(H1)+EXP(-H1)_bD _   _     _    _     !" #&$%*c2b?cf3?cf3?I?&پ? |?{G?cf3?cf3?cf3?cf3?U٩?cf3?b?cf3Habb3 (EXP(H2/Ha)-EXP(-H2/Ha))/(EXP(H2/Ha)+EXP(-H2/Ha))`ad ` a    `  a     ` a    `  a      !" #&$%Tc C?CK$?CK$??ݑ)?͐c@ ?pU?CK$?CK$?CK$?CK$?????7N?CK$?c?CK$Wave_LIF DW_Wave_L<=0 THEN 0 ELSE IF water_depth > 2* DW_Wave_L THEN DW_Wave_L ELSE DW_Wave_L*SQRT((EXP((2*PI*water_depth)/DW_Wave_L)-EXP(-(2*PI*water_depth)/DW_Wave_L))/(EXP((2*PI*water_depth)/DW_Wave_L)+EXP(-(2*PI*water_depth)/DW_Wave_L))) q CC     C          C                                               !" #&$%c5d>ս݁>`D??3$DjÌ>/22Ah R>ս݁>`D>`hi7>`hi7????>ս݁?d>`hi7Cons_ingest_Pht6Min(Phytoplankton,(PH_pref/pref_tot)*Con_pot_ingest)0-U  0    -   ! 1/h" #&$%g x??????e $Eelgrass!" #&$%!fah@ E@7Z%?@r_%px=0@r_%px=0@ E@7Z%@)Pp@7Z%????@ E h@7Z% Sed_depobIF Settling_velocity > water_depth THEN SUS_SED ELSE (Settling_velocity/water_depth)/DT * SUS_SED33@      3        3  !" #&$%f%i ????? i Sed_susp&IF Pot_Eros > 0 THEN Pot_Eros ELSE 0@3$ @    @   !" #&$%f m?@ =\@P iǦ=B7x}j Epi_Gross_PPMAX(0,Epi_prod_lim*Epiphytes)     !1=alpha .03 = plankton prim. productivity rate, 1/h What is the max Epiphyte production? Look up in Nickols 1994 In the spatial model this became negative, not sure why." #&$%c0 ?PbM񖩟?PO Q??E{r-Y?N#@!?PbM񖩟?PO Q?PO Q?PO Q ?????PbM񖩟?l?PO Q Pht_conc1IF SF_water<1 THEN 0 ELSE Phytoplankton/SF_water 0,        0   ! {g/l=kg/m^3}" #&$%cg@@@@?@@@@@@@@@@@@@@@@???@@?m@@ F_North500; !" #&$%c O?@@@@@@@@????@ ?n OMtotBio@N0R !This variable repesents the aggregate amount of biomass in the four adjacent cells that could provide mobile consumers with a food source." #&$%f ?xtj~?xws_??sIR?r٤#@q a?xtj~?xws_?xxz?xws_v?????xtj~o?xws_Pht_mortality Phytoplankton*Pht_mortality_rt00  0   !wwas: ALGAE^t_Al_Mort, unsure of its origin now: is specific mort rate, 1/h units in kg/d .005*DT = t_Phyto_Mort " #&$%$cG??{!AZ?????MD?pdecomp_temp_H2OHIF Air_temp<=0 THEN 0 ELSE (Air_temp)^4/(half_sat_H2O^4+Air_temp^4)2L          2        !" #&$%?cI.;HV .;HV .??Q)8 O~?w"5s\BL;HV .;HV .;HV .;HV .?????~;HV .?q;HV .orbital_velocityIF Wave_L<=0 THEN 0 ELSE ((Wave_height*9.8*Wave_period)/(2*Wave_L))*(1/((EXP((2*PI*water_depth)/Wave_L)+EXP(-(2*PI*water_depth)/Wave_L))/2)) cNB.. c    .   * @#     c         c          c         !" #&$%g ??????r $Consumers!" #&$%f s>-z_>Nb34??3$DjÌ>/22Ah R>-z_>Nb34>u ]>wSc@>-z_s>wSc@ Epi_consumedCons_ingest_Epi' '!" #&$%c @Y@4@4??@>@$@4@4@4@4???@4?-t@4 few_cloudsDays@B@@R@@[`@b@@f@k`@o@r@@t@v@y@{`@}@@@@@d@@@@4@4@$@$@>@>@4@$@$@$@4@4@$@$@>@>@4@$@$@$@4 !" #&$%f~>-z_>Nb34?>>vE>6>-z_>Nb34>u ]>wSc@4?ɨcA>-z_xu>wSc@PhBio_consumedCons_ingest_PhBioK< K!5 connects to consumers with Cons_ingest_PhBio {kg/d}" #&$%c?*swܚs?;LB??@TnNI@P?G?*swܚs?;LB?;LB?;LB#@ j?{'|?xv?;LBMac_NPPNPhBio_NPP+PhBio_Net_Prodz   z  !$SHORT: Net plant biomass production" #&$%f.??@*ILݕ]@^xwPhBio_litterfall/litterfall_season*Shoots*max_litterfall_rate}<,< } <    !^0.05 is max spec rate of litterfall, 1/d litterfall_season ranges from 0,1 dimensionless .05" #&$%f )?z???@?z??sdF?zH@?hs?d9R,xz?z PhBio_Resp*PHBio_resprate*Mac_Temp_Resp_Lim* Shoots o<z<  o  <  !ITemp limitation dimensionless Units kg/d PhBio specific resp rate is 1/d" #&$%c ??F?F?F?F?|1&y? | TSS_X_N@S0 !smaller than ic of SIS in cell" #&$%c)?z???Zx˧???}litterfall_seasonDayswH@2@@B@@K`@R@@V@[`@_@b@@d@f@i@k`@m@o@q@r@@sd@t@u@v@w@y@z<@{`@|@}@~@@@@@@@@d@@@@@>@H?zG{?zG{?zG{??Q?Q?Q?瑷Q?333333?(\)?\(?\(\?(\)?333333?Q?zG{?zG{?zG{??Q?Q?Q?瑷Q?333333?(\)?\(?\(\?(\)?333333?Q !" #&$%#cD??n:~X???@b_eI??n:~X?n:~X?n:~Xg????@2.#C??x~?n:~XLeaf_Limit_CoeffEif RootsRhy=0 then 0 else MAX(0,(1-Shoots/(RootsRhy*L_to_R_limit)))1<W`H 1       < 1 W     !h{SHORT: Similar to "max leaves". Computes limit of leaf growth due to root mass} can never exceed one!" #&$%c n@ @.@"# =p@"Q??@"# =p@"Q@"# =p@"# =p???@"Q?@"# =p day_lengthDays p?@@@@@@@ @"@$@&@(@*@,@.@0@1@2@3@4@5@6@7@8@9@:@;@<@=@>@?@@@@@A@A@B@B@C@C@D@D@E@E@F@F@G@G@H@H@I@I@J@J@K@K@L@L@M@M@N@N@O@O@P@P@@P@P@Q@Q@@Q@Q@R@R@@R@R@S@S@@S@S@T@T@@T@T@U@U@@U@U@V@V@@V@V@W@W@@W@W@X@X@@X@X@Y@Y@@Y@Y@Z@Z@@Z@Z@[@[@@[@[@\@\@@\@\@]@]@@]@]@^@^@@^@^@_@_@@_@_@`@` @`@@``@`@`@`@`@a@a @a@@a`@a@a@a@a@b@b @b@@b`@b@b@b@b@c@c @c@@c`@c@c@c@c@d@d @d@@d`@d@d@d@d@e@e @e@@e`@e@e@e@e@f@f @f@@f`@f@f@f@f@g@g @g@@g`@g@g@g@g@h@h @h@@h`@h@h@h@h@i@i @i@@i`@i@i@i@i@j@j @j@@j`@j@j@j@j@k@k @k@@k`@k@k@k@k@l@l @l@@l`@l@l@l@l@m@m @m@@m`@m@m@m@m@n@n @n@@n`@n@n@n@n@o@o @o@@o`@o@o@o@o@p@p@p @p0@p@@pP@p`@pp@p@p@p@p@p@p@p@p@q@q@q @q0@q@@qP@q`@qp@q@q@q@q@q@q@q@q@r@r@r @r0@r@@rP@r`@rp@r@r@r@r@r@r@r@r@s@s@s @s0@s@@sP@s`@sp@s@s@s@s@s@s@s@s@t@t@t @t0@t@@tP@t`@tp@t@t@t@t@t@t@t@t@u@u@u @u0@u@@uP@u`@up@u@u@u@u@u@u@u@u@v@v@v @v0@v@@vP@v`@vp@v@v@v@v@v@v p@"Q@"# =p@"(\)@".zG@"333333@"=p =@"B\(@"L@"W =p@"\(\@"ffffff@"p =q@"zG{@"Q@"zG@"Q@"\)@"Q@"\(@"R@"Gz@"ݣ =q@"G{@# =p @#@#(\)@#=p =@#L@#\(\@#kQ@#@#\(@# =p@#33333@#zH@#(\@#Q@$@$zG@$(\)@$=p =@$QR@$ffffff@$zG{@$\(@$ =p@$Q@$@$fffff@$G{@%\(@%# =p@%=p =@%QR@%kQ@%@%zG@%zG@%\(@%(\@%ݣ =q@& =p @&Q@&8Q@&QR@&ffffff@&@&@&zG@&zH@&(\@&\(@'\(@'(\)@'=p =@'W =p@'p =q@'Q@'Q@'Q@'R@'fffff@(@(@(.zG@(GzH@(aGz@(zG{@(\(@(\)@(\(@( =p@(ݣ =q@) =p @)# =p@)8Q@)QR@)ffffff@)@)@)zG@)zH@)Gz@)\(@*\(@*# =p@*=p =@*QR@*kQ@*@*zG@*zG@*\(@*(\@*ݣ =q@+Q@+@+333333@+GzH@+\(\@+p =q@+Q@+@+zG@+\(@+ =p@+Q@,@,zG@,# =p@,8Q@,L@,\(\@,p =q@,@,zG@, =p@,33333@,\(@, =p@,fffff@,\(@-Q@-\(@-Q@-.zG@-=p =@-GzH@-W =p@-aGz@-kQ@-zG{@-Q@-\(@-@- =p@-\)@-33333@-p =@-\(@-@-R@- =p@-(\@-Gz@-fffff@-Q@-ݣ =q@-ݣ =q@-\(@-\(@-\(@-G{@-G{@-G{@-G{@-\(@-\(@-\(@-ݣ =q@-Q@-Q@-fffff@-Gz@-(\@- =p@-@-zH@-p =@-Q@-zG@-\)@-Q@-zG@-=p @-@-u\(@-ffffff@-\(\@-L@-B\(@-333333@-(\)@-@- =p @,G{@,Q@,(\@,@,p =@,zG@,@,=p @,zG{@,ffffff@,W =p@,B\(@,.zG@,Q@, =p @+\(@+Gz@+@+Q@+ =p@+\(@+zG{@+ffffff@+QR@+=p =@+(\)@+\(@*G{@*fffff@*R@*Q@* =p@*=p @*u\(@*aGz@*GzH@*333333@*@*@)Q@)R@)p =@) =p@)\(@)u\(@)\(\@)GzH@).zG@)zG@)@(fffff@(@(Q@(Q@(Q@(kQ@(W =p@(=p =@(# =p@( =p @'\(@'(\@'\(@'zG@'zG@'zG{@'aGz@'L@'333333@'@'Q@&Q@&R@&p =@& =p@&=p @&u\(@&\(\@&B\(@&.zG@&zG@&@%fffff@%R@%Q@% =p@%=p @%u\(@%\(\@%GzH@%333333@%@%Q@$ݣ =q@$(\@$\(@$zG@$@$Q@$p =q@$\(\@$GzH@$333333@$Q@$ =p @#\(@#fffff@#R@#p =@#zG@#@#=p @#u\(@#ffffff@#W =p@#B\(@#333333@## =p@#zG@#Q@"\(@"fffff@" =p@"@"p =@"33333@" =p@"@"\(@"@"u\(@"kQ@"aGz@"\(\@"QR@"GzH@"B\(@"8Q@"333333@".zG@"# =p@"Q@"@"zG@"zG@"\(@"\(@" =p @" =p @"Q@"Q@"Q@"Q@" =p @" =p @" =p @"\(@"\(@"zG@"@"Q !" #&$%g??????$Phytoplankton!" #&$%c ??F?F?F?F?|1&y?  TSS_X_S@N0 !smaller than ic of SIS in cell" #&$%gj??????$Salt!" #&$%cA ??(q ?(q ???? TSS_X_SBIF sf_wt_X_S > 0 THEN sf_wt_X_S*TSS_conc*convert_currents ELSE 0aF4        a    !" #&$%cAR ?>dbz)>=]l+???? TSS_X_NBIF sf_wt_X_N >0 THEN sf_wt_X_N*TSS_conc*convert_currents ELSE 0aF4        a    !" #&$%c6?bM񖩟?xNJ t?333333?bM񖩟<8|;?333333?bM񖩟?xNJ t?xNJ t?xNJ t?bM񖩟? ?xNJ t TSS_conc7IF SF_water<=0 THEN 0 ELSE MAX((SUS_SED/SF_water),0 ) 38       3    ! {g/l = kg/m^3}" #&$%g y?????? $Nutrients!" #&$%c ?"?Fq6g??@v?Fq6g?Fq6g?Fq6g?????Fq6g Seed_prodDaysH@B@@R@@[`@b@@f@k`@o@r@@t@v@y@{`@}@@@@@d@@@??(\)?Q?ffffff?(\)? =p ?zG{?zG{?zG{?Q??(\)?Q?ffffff?(\)? =p ?zG{?zG{?zG{?Q !" #&$%c ?tzG{?tzG{??tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?????tzG{??tzG{ wrack_OM_PC0.005    ?tz G{!guess {Phos/carbon ratio}" #&$%cA ?>aO뼯???? TSS_X_EBIF sf_wt_X_E > 0 THEN sf_wt_X_E*TSS_conc*convert_currents ELSE 0aF4        a    !" #&$%cA ??68Sx 0 THEN sf_wt_X_W*TSS_conc*convert_currents ELSE 0aF4        a    !" #&$%f ?j[*2W?in@&??@ٛ/B?j[*2W?in@&?isQk?in@&@Q%>F|?j[*2Wx?in@& NPhBio_mort RootsRhy*NPhBio_mort_rate11  1   !kg*1/d" #&$%c?7m%M?7NPhBio_resp_rate*Mac_Temp_Resp_Lim*RootsRhy^Alpha_NPhBio o1$1  o  1 $   !P1/d*dimensionless(0-1)*(kg)^dimensionless Kr*Tlim*(NB)^a3 where a3 is unknown" #&$%c?|1&y?|1&y???\(?\(?|1&y?|1&y?|1&y?|1&y"?\(?|1&y?x?|1&yNPhBio_resp_rate0.007   ?| 1&y!# 0.17 {1/d} 0.17/24 = 0.007 {1/hr}" #&$%c  ?aQ6? ??0?WL??aQ6? ? ? X?*H`?>?x?  Mac_Prod_lim{(Mac_water_lim^switch_water_lim)*}(Mac_Light_lim^switch_light_lim)*(Mac_Nut_lim^switch_nut_lim)*(Mac_Temp_lim^switch_temp_lim) /?, /           !Y{Roel: I don't expect a water defficit} {SHORT: The exponents incorporates the switches}" #&$%c E ?D҉^?[?-ͨt?Pht_nut_limitFmin( DIN_conc/(DIN_conc+Pht_Sat_N), PO4_conc/(PO4_conc+Pht_Sat_P) )-!#0 - - !     #    !uMich-Menton, assumes only one limiting nutrient. half sat coeffs are broad approximations (2.5E-5g/l=25 g/l) {g/l}" #&$%c?tzG{?tzG{??tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?????tzG{??tzG{ Cons_PC0.005   ?tz G{!guessed {Phos/carbon ratio}" #&$%f?@ګ@8o CW^`Y????Wrack_producedPhBio_litterfallw6 w!" #&$%cB ??IMк^?IMк^?????j Salt_X_SCIF sf_wt_X_S > 0 THEN Salt_conc*sf_wt_X_S*convert_currents ELSE 0a4        a    !" #&$%g ?????? $Epiphytes!" #&$%"c X?As^[KmDÎ#`???? wrack_X_SYif Vector_Abs_Sum=0 then 0 else (wrack*(S_wrack_vector/Vector_Abs_Sum))/sqrt(cell_size)(6 D (      6  (       !" #&$%"c X?C P???? wrack_X_WYif Vector_Abs_Sum=0 then 0 else (wrack*(W_wrack_vector/Vector_Abs_Sum))/sqrt(cell_size)(6" D (      6 " (       !" #&$%c )??????  River_TSS0 !" #&$%f'> 6\>ľ8S??$"G1>3cPҤF!>3cPҤ> 6\>ľ8S>ľ(M>ľS> 6\>ľSCons_respiration(Consumers*Cons_activity*(C_resp_rt*DT)+M+  + M      !" #&$% f DX?A->rD.fh(<???? Wrack_X_outEIF on_map = 2 THEN 0 ELSE (wrack_X_E+wrack_X_N+wrack_X_S+wrack_X_W) st6<        s t      !" #&$%c ??j Salt_X_W@E0 !" #&$%c ??j Salt_X_E@W0 !" #&$%$f ?7>Ëh??6C-?6C->Ëh>Ëh???>Ëh>Ëh Epi_seeding@If (Random(0,1) < 0.1 AND Shoots>1) then Epi_Seeding_rt ELSE 0<L    < ?  <         !" #&$%g ?????? $Detritus!" #&$%c ??>ܮʵ??>ܮʵ??????y?PO4_sed_conc1 !" #&$%c ?A2B϶lC;?????wrack_X_S@N:00  !" #&$%c?M(A^>ˮʵ?#Pv>>ˮʵC?y PO4_concPO4_Maps*River_P*PO4_scen/!EM / !    !DPO4 maps * ratio of current total river input/max total river input" #&$%c ?@o$ꔅ?5ٛ=? sf_wt_X_E0l  !MAX(0, EW_vector*sqrt(cell_size)*water_depth) Surface water flow out of cell to East cell. Uses Manning's equation. Uses the Manning's roughness coeff of the donor cell. 1.49 * (HeadDiff/dist between cells)^.5 * cross-sectional area * hydraulic radius(=water depth)^.666 / Manning's roughness {ABS function because depth goes slightly negative (1e-15)} (dimensionless) * (m/m) * (m^2) * (m^.666) / (time/(m^.333)) = m^3/time units" #&$%g??????$Global!" #&$%f'??f?f?f?f????  TSS_X_in(TSS_X_E@W+TSS_X_N@S+TSS_X_S@N+TSS_X_W@E|3  |      !" #&$% f P@@"Wbg(@"|&ZAnO1wx8 Pht_Gross_PPQif SF_water=0 then 0 else Pht_Prod_lim*Pht_Gross_PP_rt*Phytoplankton^Alpha_Pht 00<          0    !91=alpha .03 = plankton prim. productivity rate, 1/h {kg}" #&$%f >ս݁>`D?3$DjÌ>/22A/>ս݁>`D>`fk9>`hi7 >ս݁>`hi7 Pht_consumedCons_ingest_Pht d0 d! {1/timestep}" #&$%g ?????? $seaweed!" #&$%c ?????????????? C_egest_eff.75I   ? ! Assimilation efficientcy (0-1)" #&$%c K>6C->6C-??>6C->6C->6C->6C-???>6C-?>6C- C_resp_rt1.25e-5   >6 C-!respiration rate 0.0003 /24" #&$%cP)b2+Q/@&@;x?]SAi0ŋ'G>P)b2+Q/@&P)b2+P)b2+Q/@&?P)b2+Pht_NPP Pht_Gross_PP-Pht_resp     ! {kg/timestep}" #&$%f &??@@VF3Է7???x PI_transfer'RootsRhy*Translocation*PI_trans_rate1q1< 1     !8COOH transfer to shoots transfer rate = 0.7/24 = 0.03" #&$%c @Y@@"@v?PrecipitationDays@$@4@>@D@I@N@Q@T@V@Y !" #&$%c Y>ʵ>ʵ??>ʵ>ʵ>ʵ>ʵ5???>ʵ?>ʵ C_mort_rt1e-6   > ʵ!6.001 mortality rate, 1/d 0.001/24=mortality rate, 1/h" #&$%c ? =p ? =p ??@2@2? =p ? =p ? =p ? =p '@.? =p ?x? =p  Mac_PP_rate0.09?   ? = p !(18=spec rate, 1/d 18/24 = 0.75 = 1/hr " #&$%c (?ƾ?Imk?Vg? Pht_Prod_lim)Pht_Temp_lim*Pht_Light_lim*Pht_nut_limit      !" #&$%c [? bM񖩟? bM񖩟??? bM񖩟? bM񖩟? bM񖩟? bM񖩟'???? bM񖩟?? bM񖩟Ingestion_rt1.25e-4-   ? bM 񖩟!(Ingestion rt consumer specific 0.003/24" #&$%c@@??@@@@???@?@ SW_pref7E !" #&$%c@@??@@@@???@?@ PH_pref7d !" #&$%cc c%0 THEN Pht_resp_rt*Pht_Temp_lim*Phytoplankton else 0  004        0    !%1.1=alpha coeff t_Pl_Resp, 1/d {kg}" #&$%c@@??@@@@???@?@ EP_pref4' !" #&$%c@@??@@@@???@?@ SH_pref2K !" #&$%c 9?zG{?zG{?*s 6?o?*s 6<5{Q>?zG{?zG{?zG{?zG{&>d5?zG{?j?zG{ Salt_conc:IF SF_water<1 THEN 0 ELSE MAX((SALT_SURF_WT/SF_water),0 ) 2>kl8        2    !'{Behm}: this is the same as salinity." #&$%f #@ x?w/ovm@ x=Ej Salt_X_out$Salt_X_E+Salt_X_N+Salt_X_S+Salt_X_Wl>k2 l >    k  !" #&$%c?????????????? RR_pref1L !" #&$%c5 @6@6??@6@6@6@6???@6?@6 pref_tot6EP_pref+PH_pref+RR_pref+SH_pref+SW_pref+Wr_pref ( E)d'KL,          (  !" #&$%c??????????????Wind_Wrack_factor.5""   ? !The additional influence wind has on speed and direction of wrack movement, besides the influence on the general water movement. " #&$%c?p)^ ?p)^ ???p)^ ?p)^ ?p)^ ?p)^ ????p)^ ??p)^  Sink_rt4.1e-3G   ?p) ^ !" #&$%f>-z_>Nb34?>>fE>6>-z_>Nb34>u ]>wSc@???>-z_x>wSc@NPhBio_consumedCons_ingest_NPhBioL1 L!" #&$%c*?MN?'F<1??@ӵS%??MN?'F<1?'F<1?'F<1????MN?x?'F<1Q10_2EXP(0.069*H20_Temp)U ( ? lD U   ! Fred - please give ref for this" #&$%c?>??@X CT???TranslocationDays@B@@R@@[`@b@@f@k`@o@r@@t@v@y@{`@}@@@@@d@@@??? =p??Q??? =p??Q !" #&$%c???????????Y???????on_map1 !Zdetermines if there is a non-null habitat in that cell 0=0ff map 1= On map, no forced inputs, no boundary cell 2= On map, no forced inputs, boundary cells 3 and above are Rivers and tidal forcing inputs (Boundary cells) 3 = Adams Point 4 = Squamscott River 5 = Lamprey River 6 = Crommet Creek 7 = Winnicut River 8 = Little Bay" #&$%c #?:.a?:E;?Ps? ^&oXBl5q? ?:.a?:E;?:.a?:.a?:E;? ?:.a Shear_res$Resistance_coeff*(1+veg_structure)' '     !" #&$%c 1??K?LWd?D@^h? Epi_prod_lim2Epi_Temp_lim*Epi_Light_lim*Epi_Nut_lim*Epi_pp_rtMj    M    !" #&$%f&?c%ϓ>?>ygt~?@mvÄ?F, NC ol*?F, N?>ϓ>?>ygt~?>ϓ>?>ϓ>?????>ygt~??>ϓ>veg_structureB(veg_dens_coeff/cell_size)*(Shoots+RootsRhy+Seaweeds+Epiphytes)  <1|,    < 1  |     !the density coefficient (stems/kg) multiplied by the density (kg/m^2), added for different types of standing plant parts units are #stems (or other structural attributes) /m^2" #&$%f 7?A9%D1}vhC???? Wrack_X_in8wrack_X_E@W:0+wrack_X_N@S:0+wrack_X_S@N:0+wrack_X_W@E:0uv6 u     v  !" #&$%c)?C???\(\??leaf_grow_seasonDays`H@2@@B@@K`@R@@V@[`@_@b@@d@f@i@k`@m@o@q@r@@sd@t@u@v@w@y@z<@{`@|@}@~@@@@@@@@d@@@@@>@H?zG{???ffffff?GzH?\(?\(\?\(\?GzH??Q??zG{???ffffff?GzH?\(?\(\?\(\?GzH??Q? !" #&$%f +,?????j Salt_X_in,Salt_X_E@W+Salt_X_N@S+Salt_X_S@N+Salt_X_W@E? 2  ?      !" #&$%c  E-y)l???1&1dRrpE-y)lE-yE-y&?`g)l?E-ymax_Sun_angle 0.4093*SIN((Days-82)/58.1) 4 ( ?1 -  R  / @M     !'{declination of sun; radians} Trust it" #&$%c?tzG{?tzG{??tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?????tzG{??tzG{Pht_PC0.005   ?tz G{!guessed {Phos/carbon ratio}" #&$%+c >?#???lu?=e]@cF??#???b?Lcl.?#?? Epi_Temp_lim?min(1, 0.22*exp((0.07-0.00001*exp(0.28*H20_Temp))*H20_Temp))Uh  , ?( \) ( ? Q - > Ëh ( ? Q U      U      !cmin(1, 0.22*exp((0.07-0.00001*exp(0.28*H20_Temp))*H20_Temp)) (H20_Temp*Epi_Temp_coeff)/Max_Temp" #&$%s @@g?@եU1K@g&BC@g&@@g@g@g%?dznM@?@g Epiphytescell_size*ic_epiphytes P'>jjs   P  !&Carbon biomass of epiphytes {kg/m^2}" #&$%c ?F?F?F?F?|1&y?  TSS_X_E@W0 !smaller than ic of SIS in cell" #&$%c ?F?F?F?F?|1&y?  TSS_X_W@E0 !smaller than ic of SIS in cell" #&$%c @Y@N#Àb1@N??@Q~==@D@N#Àb1@N@N#Àb1@N#Àb1???@N?-@N#Àb1 some_cloudsDays@B@@R@@[`@b@@f@k`@o@r@@t@v@y@{`@}@@@@@d@@@@N@Q@I@I@N@I@D@D@I@I@N@Q@I@I@N@I@D@D@I@I@N !" #&$%f< ?'縴?'`j=i?@2?̤@XuB?̤@Xu?'縴?'`j=i?'e~Ѻ)?'`j=i[?????'縴?'`j=i Epi_mort=if Shoots=0 then Epiphytes else Epiphytes*(Shoot_Mort_Ratio)<, <          !\Based on Buzzelli et al. 1999 = total shoot mortality*(ratio of epi biomass/shoot biomass)" #&$%c T>݄l>݄l??>݄l>݄l>݄l>݄ls???>݄l?->݄lk_Pht 0.0000138   >݄ l!tBehm: 0.0138 m^2 mg Pht^-1 is PAR extinction coefficient for ChlA (Buzzelli et al. 1999) =1.38e^-5 m^2 gm Pht^-1 " #&$%c  ???+?I!Ul???+?I!Ul?I!Ul?I!Ul????+??I!Ul old_cons_act H20_TempUh?UUUUUU@UUUUUU@ @UUUUUU@@@UUUUUT@!UUUUUU@#@%@'UUUUV@*h?? =p ?Gz?GzH?Q?=p =?zG?Q???zG??zG{ U!" #&$%c ? ? stress_effect@ IF Shear_stress > Shear_res THEN Shear_stress-Shear_res ELSE 0B@, B     B     !" #&$%cl@@@@?@@@@@@@@@@@@@@@@???@@?@@ F_South500; !" #&$%cu@@?@@@@@@@@???@?@F_East800;  !" #&$%c?M}H?M}H???M}H?M}H?M}H?M}H????M}H??M}HPht_Gross_PP_rt0.0009   ?M} H!" #&$%cv@I@I?@I@I@I@I@I@I@I@I???@I?@IF_West50; 2!" #&$%c>Ëh>Ëh??>Ëh>Ëh>Ëh>Ëh.???>Ëh?>ËhEpi_Seeding_rt0.00001   > Ëh!/BEHM 8/24: changed from 0.000000001 to 0.00001" #&$%c ??L??L? F21?fv?t^d??e??L??L??L??L?+??L???LT1 .833*(g_d^.375)]$  ? - ] ^ ?   !" #&$%"c&?pa#?pa#?d?Ἦ>?U ?! ,;?pa#?pa#?pa#?pa#?ַ?pa#??pa#T1a' (EXP(T1)-EXP(-T1))/(EXP(T1)+EXP(-T1))D                  !" #&$%2cJ!?˿Ճ\[?˿Ճ\[?=-~?U6?Jm?j"X?˿Ճ\[?˿Ճ\[?˿Ճ\[?˿Ճ\[?cdD1G?˿Ճ\[??˿Ճ\[T1abKIF T1a=0 THEN 0 ELSE (EXP(T2/T1a)-EXP(-T2/T1a))/(EXP(T2/T1a)+EXP(-T2/T1a))!!     !                           !" #&$%c ?z@?z@?;YTT?.2?9ȷK?J=el?5ٛ=? sf_wt_X_W0k !MAX(0, -EW_vector*sqrt(cell_size)*water_depth) Surface water flow out of cell to West cell. Uses Manning's equation. Uses the Manning's roughness coeff of the donor cell. 1.49 * (HeadDiff/dist between cells)^.5 * cross-sectional area * hydraulic radius(=water depth)^.666 / Manning's roughness {ABS function because depth goes slightly negative (1e-15)} (dimensionless) * (m/m) * (m^2) * (m^.666) / (time/(m^.333)) = m^3/time units " #&$%c ?pbM񖩟?pbM񖩟???pbM񖩟?pbM񖩟?pbM񖩟?pbM񖩟????pbM񖩟??pbM񖩟 Epi_resp_rt0.004A   ?pbM 񖩟!" #&$%c m4@D333333"<؆7j E&.nA^@7yQ4sՆw3333333"<؆"<؆"<؆@C~"."333333?"<؆ Air_tempDays5po h? @ @@ @ @@@ @ @@@@;;@ @ @@" Ԭ@$@&xֶx@( @*I$I"@,;8@.QN@0 @ @@1 K@2 Ԭ@3 \\@4@5N@6xֶx@7-,@8 @9Q@:I$I @;/+@<;6@=eFTeA@>QL@?\W@@ @ @@@@A K@ANN@B Ԭ@B-@C \\@C٦m@D@Djj@EN@E@Fxֶx@F/@G--@Gƽp @H@H;;@IQ@IΕ#@JI$I)@J*2/@K/5@KW5sW;@L;A@L @ G@MeFTeM@MKS@NQY@NsW5s_@O\e@O'bv'k@P @ @@PKmm@P@P肹{@Q K@QL!!@QNN@Q{{@R Ԭ@RL\@R-@R//@S \\@SMůŧ@S٦m@S=@T@TN==@Tjj@T~@UN@UN@U@UKK@Vxֶy@VO_@V/@V @W-- @WPZZ@Wƽp@Wܫ @@X@XQ@X;;@Xhh!@YQ$@YQ!'@YΕ*@Y!-@ZI$I0@ZRv'bv3@Z*26@Z-9@[/<@[S*2*?@[W5sWB@[8CE@\;H@\S=K@\ @ N@\8C8Q@]eFTeT@]TI$W@]KZ@]N]@^Q`@^UFTeFc@^sW5sf@^Zi@_\l@_U_o@_'bv'r@_TeFTu@` @ @@`+W5sW=@`Kmm@`k8C?@`@`DZ;A@`肹{@`=C@a K@a, @ E@aL!!@al8C8G@aNN@aeFTeI@a{{@aI$K@b Ԭ@b,KM@bL\@blNO@b-@bQQ@b//@bFTeFS@c \\@c-sW5sU@cMůŅ@cmZW@c٦m@c\Y@c=@c_[@d@d.'bv']@dN==@dnTeFT_@djj@dha@d~@djc@eN@e.me@eN@eopg@e@e5sW5i@eKK@ebv'bk@fxֶx@f/xֶm@fO_@fo{o@f/@f~q@fؘ@fhs@g-,@g0C8Cu@gPZY@gpppw@gƽp@gůy@gܫ @Ҝ@gԬ{@h@h0ֶx}@hQ @hq$I$@h;;@hQQ@hhh@h~~@iQ@i1Ƿ@iQ!@iqů@iΕ@iZ@i! @i2*2@jI$I @j2__@jRv'bv@jrԬԶ@j*2@j@j-@jj@k/@k3;@kS*2*@ks@ @@kW5sW@kmm@k8C@k@l;@l3{@lS=@lsK@l @ @l!!@l8C8"@lNN@meFTe$@m4{{@mTI$&@mtԬ@mK(@m\ͩ@mN*@m-@nQ,@n5//@nUFTeF.@nu\\@nsW5s0@nůű@nZ2@nm@o\4@o5=@oU_6@ov@o'bv'8@o==@oTeFT:@ojj@p @ @@pKK@p+W5sW @p;bv'ba@pKmm@p[xֶx@pk8C$@p{xֶe@p@p_@pDZ;(@pȹ{i@p肹{@p/@p=,@p~m@q K@qؕ@q, @ 0@q<hq@qL!!@q\-,@ql8C84@q|C8Cu@qNN@qZY@qeFTe8@qppy@q{{@qp@qI$<@qů}@r Ԭ@r @@r,K@@r<Ԭ@rL\@r\@rlND@r|ֶxַ@r-@r@rQH@r$I$@r//@r;; @rFTeFL@rQQ@s \\@shh@s-sW5sP@s=~~@sMů@s]Q@smZT@s}@s٦m@s!@s\X@sů@s=@s@s_\@sZ@t@t!@t.'bv'`@t>2*2@tN==@t^I$I#@tnTeFTd@t~__@tjj@tv'bv'@thh@tԬԩ@t~@tޣ*2+@tjl@t@uN@u-/@u.mp@u>j@uN@u^/3@uopt@u;@u@u*2*7@u5sW5x@u@ @@uKK@uW5sW;@ubv'b|@umm@vxֶx@v8C?@v/xֶ@v?@vO_@v_;C@vo{@v{@v/@v=G@v~@vK@v  h333333"&333333!333333 ffffff!"ffffff(ffffff033333(ffffff*$%333333 .1333333033333*"&333333! ffffff! ffffff333333$'ffffff)*3ffffff4%333333333333!"ffffff333333 ffffff333333"""ffffff@ffffff@ffffff333333!$"ffffffffffff333333 ffffffffffff333333333333@@"@333333@ffffff@ffffff@@@@ffffff@@ffffff@ ffffff333333 ffffff ffffff@"@)@$@@ffffff@@ ffffff@ffffff@ @@(ffffff@(ffffff?333333@!@333333@ ffffff?333333333333@@ffffff@%333333@ @333333@*@&333333@'ffffff@3ffffff@1333333@)@&333333@333333@%333333@%333333@)@'ffffff@ffffff@ffffff@$@1@033333?333333@@&333333@$@2fffff@3ffffff@7L@8ffffff@4@3ffffff@ @ffffff@'ffffff@&333333@(ffffff@'ffffff@*@'ffffff@.@,@1@2L@2fffff@2fffff@1@033333@4@4@1@*@ @2fffff@2L@@ffffff@,@1@3ffffff@1333333@)@ @3ffffff@1@1@0@.@4@2L@/333333@/333333@/333333@033333@033333@*@033333@6333333@4@0@,@.@2fffff@7L@6@3ffffff@2L@.@)@0@.@/333333@.@0@2L@/333333@,@,@4@1333333@%333333@/333333@,@.@1@3ffffff@7L@8ffffff@7L@5@1@1333333@/333333@'ffffff@%333333@*@(ffffff@(ffffff@333333@ffffff@,@/333333@+@'ffffff@%333333@&333333@.@+@*@033333@033333@1@.@$@@,@)@.@*@@333333@ffffff@ffffff@@$@ @@!@@)@033333@*@333333@"@&333333@@ @%333333@333333@@ ffffff@ @ffffff333333@@ ffffff@333333@ffffff@ffffff@"@ @333333@@ffffff@333333@@ffffff?333333@@ ffffff?333333@@ffffff333333@@*@333333?333333?@@333333"333333@?333333?333333@&333333@ffffff$@333333 ffffff ffffff333333ffffff@@ ffffffffffff!!333333333333ffffff@ffffff@ffffffffffff ffffff!$(ffffff333333!ffffffffffff %333333/3333334%333333 !" #&$%c ?????? River_Con07 !" #&$%c ?????????????? Alpha_Epi1.1   ? !" #&$%c !?!? =@-SQ>?@@@@a????@?@ water_depth3]hqc| ~03 !busing new grid output from hydrology model Equation was: Max(0, (Tide+Datum_Diff)-ic_sed_elev)" #&$%c  >6C->6C-??>6C->6C->6C->6C-???>6C-?>6C- Epi_Sat_N2.5E-5M   >6 C-!" #&$%c  >Ëh>Ëh??>Ëh>Ëh>Ëh>Ëh???>Ëh?>Ëh Epi_Sat_P2.5E-6M   > Ëh!" #&$%c@UW3o @jDͣ??@X,gZ?'j)@UW3o @jDͣ@UW3o @UW3o ???@3\pS>?-@UW3o RandRANDOM(0,100)  d   !" #&$%c??????yRiver_N0- !" #&$%c ?PbM񖩟?PbM񖩟???PbM񖩟?PbM񖩟?PbM񖩟?PbM񖩟????PbM񖩟??PbM񖩟 Seed_wgt0.001H   ?PbM 񖩟!guessed in kg per individual" #&$%9c n[(?#?????#???@?07Ik?#?x? Mac_Temp_limoif switch_temp=1 then 1-(40-Air_temp)/40 else min(1, 0.22*exp((0.07-0.00001*exp(0.28*H20_Temp))*H20_Temp))wU(( w      (   (   (   , ?( \) ( ? Q - > Ëh ( ? Q U      U      !ASHORT: Air or Water temp limited depending on switch PTS Thesis " #&$%&c0@""@"$'??@-I5@"zUk@""@"$'@""@""A@*,j@"$'?@"" L_of_Day1 24*((PI/2)-ARCTAN(Loc_SA/SQRT(1-Loc_SA^2)))/PI T                    !Btotal day length in radians 24 x fraction of day dawn til noon/pi" #&$%4c '#?T?-QNB????c8?T?-QNB?-QNB?-QNB????T?-?-QNB Reflection0.5*((((SIN(Zenith-Angle_of_refraction))^2) /((SIN(Zenith+Angle_of_refraction))^2))+ (((TAN(Zenith-Angle_of_refraction))^2)/((TAN(Zenith+Angle_of_refraction))^2))) f##  ? f       f        f       f          !" #&$%c ?tzG{?tzG{???tzG{?tzG{?tzG{?tzG{ ????tzG{??tzG{ sw_pp_rt0.005c   ?tz G{! was 0.00125" #&$%c @Cఄ+@C??A̬Gh@࣎@Cఄ+@C@Cఄ+@Cఄ+@N:~@C?x@Cఄ+ mac_tot_biom Shoots+RootsRhy<1}  < 1  !" #&$%c o@@?@@@@@@@@????@@?@ ic_sed_elev3 !IF there is water BELOW MSL then use the bathymetry data to determine elevation of sediments ELSE use the land elevation data indicating elevation above MSL plus the datum below mean sea level ALL VALUE ARE POSITIVE" #&$%c |?@^Q!X@ll)+K?5ٛ=? sf_wt_X_S0 !MAX(0, -NS_vector*sqrt(cell_size)*water_depth) Surface water flow out of cell to East cell. Uses Manning's equation. Uses the Manning's roughness coeff of the donor cell. 1.49 * (HeadDiff/dist between cells)^.5 * cross-sectional area * hydraulic radius(=water depth)^.666 / Manning's roughness {ABS function because depth goes slightly negative (1e-15)} (dimensionless) * (m/m) * (m^2) * (m^.666) / (time/(m^.333)) = m^3/time units " #&$%cT??PbM񖩟?PbM񖩟?PbM񖩟?PbM񖩟?????PbM񖩟? SW_X_N@S0 !" #&$%co??PbM񖩟?PbM񖩟?PbM񖩟?PbM񖩟?????PbM񖩟? SW_X_S@N0 !" #&$%c >6C->6C-??>6C->6C->6C->6C- ???>6C-?>6C- sw_resp_rt1.25e-5x   >6 C-!!Sea weed respiration rate (kg/h)" #&$%c  ?j  Salt_X_S@N0 !" #&$%c9 ??*;=i???? SW_X_E:IF sf_wt_X_E>0 THEN seaweed_rate*sf_wt_X_E*SW_conc ELSE 0~e4        ~    !" #&$%Oc>@R@@R@?@^֙=@uÉҡ@f~U?qoӑ@R@@R@@R@@R@????@R@? @R@Current_directionIF NS_vector=0 AND EW_vector=0 THEN 0 ELSE (IF (450-(ARCTAN(NS_vector/EW_vector)*(360/(2*PI)))) >360 THEN (450-(ARCTAN(NS_vector/EW_vector)*(360/(2*PI))))-360 ELSE (450-(ARCTAN(NS_vector/EW_vector)*(360/(2*PI)))))7>> 7       >   7     h       h  0   7     h       h  >   7     h       !" #&$%c ?`bM񖩟?`bM񖩟???`bM񖩟?`bM񖩟?`bM񖩟?`bM񖩟????`bM񖩟? ?`bM񖩟 sw_mort_rt0.002   ?`bM 񖩟!Behm: was 2e-3" #&$%c ?@RIP_@kG'O?5ٛ=? sf_wt_X_N0> !MAX(0, NS_vector*sqrt(cell_size)*water_depth) Surface water flow out of cell to East cell. Uses Manning's equation. Uses the Manning's roughness coeff of the donor cell. 1.49 * (HeadDiff/dist between cells)^.5 * cross-sectional area * hydraulic radius(=water depth)^.666 / Manning's roughness {ABS function because depth goes slightly negative (1e-15)} (dimensionless) * (m/m) * (m^2) * (m^.666) / (time/(m^.333)) = m^3/time units " #&$%c9 ??qU ?qU????SW_X_W:IF sf_wt_X_W>0 THEN seaweed_rate*sf_wt_X_W*SW_conc ELSE 0~e4        ~    !" #&$%c9 ??PK?ęQ????SW_X_S:IF sf_wt_X_S>0 THEN seaweed_rate*sf_wt_X_S*SW_conc ELSE 0~e4        ~    !" #&$%c??BY y|?BY y|?????E_Current_Vector0 !" #&$%c ?N?C_JJ?ݝӲ K??? V46C->6C-??>6C->6C->6C->6C-???>6C-?>6C- sw_Sat_N2.5e-5 E   >6 C-!" #&$% c !?L칭4?L칭4?:0g?S?9Q=?XPi'p?L칭4?L칭4?L칭4?L칭4?MG(?L칭4??L칭4 Wave_height" (.283*Wind_speed^2*Ha*Habb)/9.7abq< ( ? 1&     a  b  / @#ff ffff  !" #&$%c; ??ɬK?v}????SW_X_N 0 THEN seaweed_rate*sf_wt_X_N*SW_conc ELSE 0~e4        ~    !" #&$%&c 4?9y"?9y"?n^K?soE??Aŷ?EXo8?9y"?9y"?9y"?9y"?6?9y"??9y" Wave_period5IF T1a=0 THEN 0 ELSE (7.54*Wind_speed/9.8)*T1ab*T1a qT       ( @( \)   / @#      !" #&$%cm@v@o@o@q@i@t@K@o@o@o@o@@o?@oWind_directionDays; h?@@@@@@@ @"@$@&@(@*@,@.@0@1@2@3@4@5@6@7@8@9@:@;@<@=@>@?@@@@@A@A@B@B@C@C@D@D@E@E@F@F@G@G@H@H@I@I@J@J@K@K@L@L@M@M@N@N@O@O@P@P@@P@P@Q@Q@@Q@Q@R@R@@R@R@S@S@@S@S@T@T@@T@T@U@U@@U@U@V@V@@V@V@W@W@@W@W@X@X@@X@X@Y@Y@@Y@Y@Z@Z@@Z@Z@[@[@@[@[@\@\@@\@\@]@]@@]@]@^@^@@^@^@_@_@@_@_@`@` @`@@``@`@`@`@`@a@a @a@@a`@a@a@a@a@b@b @b@@b`@b@b@b@b@c@c @c@@c`@c@c@c@c@d@d @d@@d`@d@d@d@d@e@e @e@@e`@e@e@e@e@f@f @f@@f`@f@f@f@f@g@g @g@@g`@g@g@g@g@h@h @h@@h`@h@h@h@h@i@i @i@@i`@i@i@i@i@j@j @j@@j`@j@j@j@j@k@k @k@@k`@k@k@k@k@l@l @l@@l`@l@l@l@l@m@m @m@@m`@m@m@m@m@n@n @n@@n`@n@n@n@n@o@o @o@@o`@o@o@o@o@p@p@p @p0@p@@pP@p`@pp@p@p@p@p@p@p@p@p@q@q@q @q0@q@@qP@q`@qp@q@q@q@q@q@q@q@q@r@r@r @r0@r@@rP@r`@rp@r@r@r@r@r@r@r@r@s@s@s @s0@s@@sP@s`@sp@s@s@s@s@s@s@s@s@t@t@t @t0@t@@tP@t`@tp@t@t@t@t@t@t@t@t@u@u@u @u0@u@@uP@u`@up@u@u@u@u@u@u@u@u@v@v@v @v0@v@@vP@v`@vp@v@v@v@v@v@v h@o@l@j@m@i@n @q0@q@k@m@p@i @o`@t@q@n@f@@i@b`@g@p@g@r@q@o@@b@d`@p`@n @i@p0@k@@n @p@g@@q@r@l @p@h@k@r@q@n @f@rp@q@h@U@e`@g @c@qP@k@m`@sP@q0@p @n @l@rP@f @g @j@` @o@rp@i@i@qp@j@T@g`@e@@a@o@j`@g@@qp@p@j@g @p@i @g @f@pP@h@b@r0@h@n@q@i @Y@a@P@@b`@g@p@c @f@h@e@f@@h@a`@Y@@a@@c@@a`@o @e@c@]@\@d @c@c@]@f@p@r@p@d@S@V@a@@a`@g @]@[@a@`@i @`@^@\@a@V@a@U@@_@e`@c@@Y@e@k@i`@o@`@@`@m@pp@g@d@d@^@h@q@j@d@Q@k@@j@j@d @q@p@o@@_@@e@a@d@@k@q@pp@c@l`@p@r @o@h@n@k@h@a@c@o @p@q@q0@c@\@^@qp@h`@j`@d@@V@g@p@c`@g@i`@g@f @e@f@]@W@b@e`@f@@d`@q@k@g@f@h@n@k@@j@i@h@c @e@k@b`@N@i@c@e@@q@i@n@qP@l @k`@r0@l@m@@g@h@b @K@`@@d@@k@p`@i@@k@@h@h@g@p@q@pP@l@q@k@@i@^@`@g@c@a@p@j@a`@k@a@h@p@b@f @o@r@o@m@@g@qP@s@e@@f@h@c@n@j@@Y@o@m@d`@j @j`@o@sp@k@c @j@p`@i@p@k`@k @k@@k@o@g@p@j@g@l@q@f@m`@p@n`@k@g@j`@j@h@k`@i@q@h @l @r@l@q@k@m@@tP@tP@j@e@m@o@k@n @n@d@l`@l @n@@s@l`@m@@b@`@d@g @s@q@i@l @g@b@qP@q@n @o @i@@t  !" #&$%c @4@#Q@#Q@A@)fffff@!Q@4!5`@ wJ8v@#Q@#Q@#Q@#Q@zG@#Q?@#Q Wind_speedDays ]^?@--@ZZ@8C8D@ZZ@pp@pp@ NN@"ZY@$TeFTd@&po@({z@*pp@, $I$@. ů@0FTeFS@1ZY@2__@3TeFTe@4jk@5pq@6bv'bw@7{}@8h@9pp@:Ԭ@; $I$@< ~~@= ů@> 2*2@? Ԭ@@sW5sY@@Z\@A\_@A_b@B'bv'e@BTeFTh@Chk@Cjn@Dmq@Dpt@E5sW5w@Ebv'bz@Fxֶ}@F{@G~@Gh@HC8C@Hpp@Iů@IԬ@Jֶx@J$I$@K QQ@K~~@L @Lůů@M Z@M2*2@N __@NԬԡ@O @Oɾj@Půņ@PEZ]@Pm@P\_@Q=@QE_a@Q@Q'bv'c@R==@RFTeFTe@Rjj@Rhg@S~@SFji@SN@Smk@T@TGpm@T@T5sW5o@UKK@UGbv'bq@Uxֶx@Uxֶs@V_@VG{u@V/@V~w@Wؿ@WHhy@W-,@WC8C{@XZY@XHpp}@Xýp@Xů@Y @@YHԬ@YÇ@Yֶxְ@Z @ZI$I$@Z;;@ZQQ@[ hh@[I~~@[Q @[@\ ! @\Iů@\ŕ@\Z@] !@]J2*2@]I$I@]__@^ v'bv@^JԬ@^ɣ*2@^@_ -@_Jj@_/@_;@`Q@`%ZO@`E@`em@`!@`\S@`ů@`=@a@a%_W@aFZ@af@a!@a'bv'[@a2*2@a==@bI$I@b&TeFT_@bF__@bfjj@bv'bv"@bϪhc@bԬ@b~@c*2&@c&jg@cFլ@cfN@c-*@cmk@cj@c@d/.@d'po@dG;@dg@d*2*2@d5sW5s@d@ @@dKK@eW5sW6@e'bv'bw@eGmm@egxֶx@e8C:@exֶ{@e@e_@f;>@f'{@fG{@fg/@f=B@f~@fK@f@g @ F@g(h@gH!!@gh-- @g8C8J@gC8C@gNN@gZZ @heFTeN@h(pp@hH{{@hhp@hI$R@hů@hԬ@h @@iKV@i(Ԭ@iH\ͯ@ih@iÓNZ@iֶx@i-@i@j Q^@j)$I$@jI//@ji;;!@jFTeFb@jQQ@j\\@jhh%@k sW5sf@k)~~@kIůŏ@kiQ)@kZj@k@km@k!-@l \n@l)ů@lI=@li1@lŜ_r@lZ@l@l!5@m 'bv'v@m*2*2@mJ==@mjI$I9@mTeFTz@m__@mjj@mv'bv=@n h~@n*Ԭ@nJ~@nj*2A@nɨj@n@n怓N@n-E@o m@o*j@oJ@oj/I@op@o;@o @o*2*M@p@pZg@p%_@p5Ǥ@pE;G@pUm@pe{@pu!'@p{@p\g@p/@pů@p=G@p͋=@p~@p'@qK@q_g@q&@q6Z@qF @ G@qV@qfh@qv!'@q!!@q'bv'g@q--@q2*2@q8C8G@q==@qC8C@qI$I'@rNN@rTeFTg@r&ZZ@r6__@rFeFTeG@rVjj@rfpp@rvv'bv'@r{{@rhg@rϽp@rԬԤ@rI$G@r~@rů@r*2'@sԬ@sjg@s& @@s6դ@sFKG@sVN@sfԬ@sv-'@s\͂@smg@sχ@sj@sNG@s@sֶxֽ@s/'@t-@tpg@t'@t7;@tGQG@tW@tg$I$@tw*2*'@t//@t5sW5g@t;;@t@ @@tFTeFG@tKK@tQQ@tW5sW'@u\\@ubv'bg@u'hh@u7mm@uGsW5sG@uWxֶx@ug~~@uw8C'@uůł@uxֶg@uQ@u@uZG@u״_@uǽ@uֱ;'@vm@v{g@v'!@v7{@vG\G@vW/@vgů@vw='@v=@v~g@v@vK@v_G@v؍@vZ@v @ '@w@whg@w(!@w8!!@wH'bv'G@wX-,@wh2*2@wx8C8'@w==@wC8Cg@wI$I@wNN@wTeFTG@wZY@w__@weFTe'@xjj@xppg@x(v'bv@x8{{@xHhG@xXp@xhԬԽ@xxI$'@x~@xůg@x*2@xԬ@xjG@x @ҍ@xս@xK'@yN@yԬg@y(-@y8\ͤ@yHmG@yX@yhj@yxN'@yÖ@yֶxg@y/@y-@ypG@y @y;@yQ'@z @z$I$g@z)*2*@z9//@zI5sW5G@zY;:@zi@ @@zyFTeF'@zKK@zQQg@zW5sW@z\\@zbv'bG@zhg@zmm@zsW5s'@{ xֶx@{~~g@{)8C@{9ůŤ@{IxֶG@{YQ@{i@{yZ'@{Ŵ_@{g@{;@{զm@{{G@{!@{{@{\'@| /@|ůg@|)=@|9=@|I~G@|Y@|iK@|y_'@|؂@|Zg@| @ @|@|hG@|!@|!!@|'bv''@} -,@}2*2g@}*8C8@}:==@}JC8CG@}ZI$H@}jNN@}zTeFT'@}ZY@}__g@}eFTe@}jj@}ppG@}v'bu@}{{@}h'@~ p@~Ԭg@~*I$@~:~@~JůG@~Z*2@~jԬ@~zj'@~ɫ @҂@~g@~K@~N@~ԬG@~-@~\ͽ@~m'@ @jg@*N@:@JֶxG@Z/@k-@{p'@ @;g@Q@@$I$G@*2)@//@5sW5'@ůd@ Z4@*2@_@%Ԭ@-t@5jD@=;@E @Ҋ@Mm@Uճ@]{T@eK$@m!@uN@}{@Ԭd@\4@-@/@\@Яůt@mD@=@@=@;j@~T@N$@@@K@ֶxd@ _4@/@@&-@.Zt@6pD@> @ @F @N@V;@^hT@fQ$@n!@v@~!!@$I$d@'bv'4@*2*@-,@//@2*2t@5sW5D@8C8@;:@==@@ @@C8CT@FTeF$@I$H@KK@NN@QQd@TeFT4@W5sW@ZY@&\\@.__t@6bv'bD@>eFTe@Fhg@Njj@Vmm@^ppT@fsW5s$@nv'bu@vxֶx@~{{@~~d@ڪh4@8C@p@ů@Ԭt@xֶD@I$@Q@~@@ůT@Z$@*2@_@߬Ԭ@d@j4@;@ @@&m@.t@6{D@>K@F!@NN@V{@^ԬT@f\$@n-@v/@~\@ůd@m4@=@@ϋ=@jt@~D@N@@떕@K@ֶxT@_$@/@؀@-@Zd@p4@ @ @ @'@/;t@7hD@?Q@G!@O@W!!@_$I$T@g'bv'$@o*2)@w-,@//@2*2d@5sW54@8C8@;:@==@@ @t@C8CD@FTeF@I$H@KK@NN@QQT@TeFT$@W5sV@ZY@\\@__d@bv'b4@eFTe@hg@'jj@/mmt@7ppD@?sW5s@Gv'bu@Oxֶx@W{{@_~~T@gh$@o8C@wp@ů@Ԭd@xֶ4@I$@Q@~@t@ůD@Z@*2@촜_@׬Ԭ@T@j$@;@֫ @Ҁ@ئm@d@{4@K@!@'N@/{t@7ԬD@?\@G-@O/@W\ͳ@_ůT@gm$@o=@w@=@jd@~4@N@ĕ@@Kt@ֶxD@_@/@ث@#Q@!Q@&uU.r@"(I@%Ck&@)&fffff@"@&{+@@)fffff@'33333@{ȴ];@*wY@ȴ];@-:&@%hÿa @wم@33333@33333@r)@ȴ];@ U`@ȴ];@'9@!"-Qq@#D9@4!5`@0[Ȱ,@)UU.r@# @ @&33333@DZ@#ffffff@ @@@)^ t@33333@On@s33333@ J,@@U>@'o@ Q@"qQ@ @@%hÿa @#I@ J,@"Y@,Ck&@!D9@@&wY@&33333@!wlH'@wم@wم@DZ@j J,@Hr)@" ̭@!U`@0l@.@@ J,@!1Q@#@@U>@" @%Ck&@@DZ@'o@.$Ck&@.7@&]nY@@*@' @<}_@(wwY@#D9@&ÿa @$Y@#I@%fffff@'@'fffff@&:&@<}_@&"Ξ?@ wlH'@%U.r@!hI@@b" ̭@!@'o@&wY@{ȴ];@ J,@ UU`@$Ffffff@)*:&@"I@U>@#ffffff@#Ȯ |@(@,7@)&fffff@7wم@"{Z@33333@@fffff@">@(wY@(17@!N>@"33333@&fffff@ J,@"WwlH'@$33333@@ @.nY@$nY@'o@!@@wم@ȴ];@#D9@">@!wlH'@'o@ȴ];@<}_@"B"-Qq@"&fffff@#Ȯ |@ fffff@"" ̭@ (I@ȴ];@!I@'l@(Ck&@':&@]<}_@'o@33333@wم@&Ffffff@"WwlH'@U>@!}@@ Ȯ |@ȴ];@#N>@+^ t@%N^ t@@@ D9@j J,@U>@&uU.r@*7wY@%+@@On@@"*{Z@#Ȯ |@$:&@$U.r@7wم@#5U`@@UU>@DZ@ I@U>@#Ȯ |@b" ̭@'o@@r)@!@ {Z@@@nOn@ Y@ȴ];@33333@33333@ fffff@%""Ξ?@#Ȯ |@'o@fffff@ffffff@%@j J,@ @@UU>@<}_@j J,@ ɻ{Z@" ̭@ {Z@wم@.On@.On@!Ȯ |@On@@@<}_@UU>@ @On@$*:&@@@ I@ =@0fffff@(17@ qQ@&17@U>@"" ̭@ J,@DDZ@33333@* J,@#@<}_@On@L@DZ@* J,@ȴ];@@%ÿa @$ÿa @!@$U.r@ fffff@+DCk&@fffff@U>@U>@33333@On@33333@<}_@On@!D9@On@!@@@<}_@""-Qq@$Ffffff@;ȴ];@.On@ UU`@"" ̭@*dCk&@" ̭@%nY@&333333@"D9@wم@" ̭@U>@j J,@%UU.r@&fffff@%j:&@!U`@U>@Y@j J,@b" ̭@)Ck&@"uU`@Hr)@* J,@!`@ J,@DZ@&wwY@&fffff@%^ t@'o@DZ@(33333@<}_@!33333@* J,@wم@"" ̭@%Ck&@ wJ8v@{ȴ];@)s33333@$@#@wم@{ȴ];@"" ̭@,@!"-Qq@2:9m@(33333@* J,@)Y@"" ̭@!}@#{Z@s33333@{ȴ];@)@#j{Z@@/dCk&@("Ξ?@'UU.r@"=@!Ȯ |@*+@@(wY@.On@@ @%nY@ @"33333@r)@ n>@$Ck&@-Hÿa @ I@!5U`@"" ̭@!Ȯ |@'o@&7@7wم@ J,@ďh3@<}_@[@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk&@+DCk& !was graph of 6*days" #&$%c %>Ëh>Ëh??>Ëh>Ëh>Ëh>Ëh???>Ëh?>Ëh sw_Sat_P2.5e-6E   > Ëh!" #&$%c??sͼ?-V?sͼ?-V?????N_Current_Vector0 !" #&$%c0?*6C-?*6C-???*6C-?*6C-?*6C-?*6C-????*6C-??*6C-Pht_mortality_rt2e-4o   ?*6 C-!" #&$%c ?Ul΄;???({ @M&$L-l΄;l΄;l΄;????-l΄;DeclineSDCLNgbfH????H!.Hǫh/{|?|?h/{?ǫ?!.H g!ASIN function of the argument" #&$%c 6?"Y׮?"Y׮???"Y׮?"Y׮?"Y׮?"Y׮????"Y׮??"Y׮ Pht_resp_rt1.4e-4   ?"Y ׮!" #&$%c?Ëh?Ëh???`bM񖩟?`bM񖩟?Ëh?Ëh?Ëh?Ëh$?????`bM񖩟?Ëh?x ?ËhPHBio_resprate0.00008z   ? Ëh!%0.002 {1/d} 0.002/24 = 0.00008{1/hr}" #&$%c <>6C->6C-??>6C->6C->6C->6C-???>6C-?!>6C- Pht_Sat_N2.5e-5   >6 C-!" #&$%cF @&@gEב@&@gEב???j6CÃZH@&@gEב@&@gEב@&@gEב@&@gEב???@&@gEב?"@&@gEבW_wrack_vectorGIF W_vector<0 THEN 0 ELSE W_vector+(Wind_Wrack_factor*EW_wind_vector)(4            !" #&$%c @>Ëh>Ëh??>Ëh>Ëh>Ëh>Ëh???>Ëh?#>Ëh Pht_Sat_P2.5e-6   > Ëh!" #&$%cV@a@a??@a@a@a@a???@a?-%@aIk_SW140 !" #&$%c F?+jg?+jg???+jg?+jg?+jg?+jg????+jg? &?+jg settling_rt0.0468   ?+ jg!in meters per hour" #&$%cH?6C-?6C-???6C-?6C-?6C-?6C-????6C-? '?6C-Resistance_coeff 0.0001   ?6 C-!" #&$%c?????????????(? Wr_pref1) !" #&$%c-??@n )????)Cons_ingest_Wr.Min(wrack,(Wr_pref/pref_tot)*Con_pot_ingest)6(-Ur  6 (   -   !*kg*1/h*kg*(Ingestion_rt*DT)" #&$%s +?3O1?333333@UQ/@4H)=@4?3O1?333333?3O1?3O1/>M?333333?+?3O1 Consumersic_consumer *cell_size 4 HJK7\-U5IG  4   !0Carbon biomass of an aggregated consumer (kg). " #&$%c1v ??????x,PhBio_Lttr_ratio2IF Shoots=0 THEN 0 ELSE (PhBio_litterfall/Shoots)<w, <      w <  !" #&$%c5{>D! 3>540??>D! 3>540>R 4<>R 4<???>D! 3?->R 4<Con_pot_ingest6MIN(OMtotBio,(Ingestion_rt*Cons_activity)*Consumers)>M+ E)d'KL  >  M  +   !" #&$%c !? "?՟@???? "?՟@?՟@?՟@E??U??-/?՟@Mac_Light_lim"PAR_Shoots/(PAR_Shoots+Ik_Shoot)      !FRoel's formulation greatly reduces NNP min(Light_to_Plants_NET*4,1)" #&$%s 0@>@= @V3/@BlRؔADz5T@>@= @= @= l?|se@>?0@= Phytoplanktonwater_depth*cell_size*ic_Pht 0 d>loo    0  !mCarbon biomass of algae, including phytoplankton and filamentous algae (kg). water_depth*cell_size*ic_Pht" #&$%s @"@!NC?@ ws2o@Aj@"@!NC@!NC@!NCE?p㚍@"?x1@!NC RootsRhycell_size*ic_mac_NPhBio >L{O~B     !FCarbon biomass of the non-photosynthetic portion of macrophytes (kg)." #&$%s 2@@@g%@rzI0@g%=(@@@@g@?j2@ SALT_SURF_WTSF_water*ic_SALT :^X^   :  !hsalt (kg) dissolved in surface water SF_water{m^3}*ic_SALT{ppt,gms/l adjust this value for IC changes}" #&$%s3@n@f@Ц@rlRؔ=0@rlRؔ=0@n@f@Ц@f@Ц@f@Ц@7Z%@n? 3@f@ЦSUS_SEDSF_water*ic_TSS hiFh     !Total suspended Seds Range of values g/l: source: Roel, Lousiania data mean 0.061 max 1.623 min 0.008 { g/l same as kg/m^3}" #&$%c @9@9??@9@9@9@9???@9?4@9half_sat_Air255 !" #&$%$cG??????5decomp_temp_AirHIF Air_temp<=0 THEN 0 ELSE (Air_temp)^4/(half_sat_Air^4+Air_temp^4)4L          4        !" #&$%s 6A1q C[W?6wrack ic_wrack *cell_size )G>tsH rHG     !XOrganic carbon mass that is incapable of being suspended but floats on top of the water" #&$%c ?携˥q?携˥q???携˥q?携˥q?携˥q?携˥q????携˥q?7?携˥q NS_vectorN_vector+S_vectorn    n  !" #&$%c [?TzG{?TzG{???TzG{?TzG{?TzG{?TzG{ ????TzG{?8?TzG{ Travel_time1.25e-3QRSXYZ[V   ?Tz G{!!cell width per timestep 0.03/24=" #&$%s.?AMUvCނÂ????9Seeds0DIHID !" #&$%c?zG{?zG{??zG{?zG{?zG{?zG{?zG{?zG{?zG{?zG{r?????6C-?zG{?j:?zG{ ic_SALT0.022   ?z G{!sppt,gms/l adjust this value for IC changes. Habitat dependant or from GIS data file for water quality measurments" #&$%2c;!@I@I@@@I@@I@I@I@I@I@T@I?;@IFetchIF Wind_direction < 45 THEN F_North ELSE IF Wind_direction <135 THEN F_East ELSE IF Wind_direction < 225 THEN F_South ELSE If Wind_direction < 315 THEN F_West ELSE F_North m^!!  -    m !       !       !   ;     !  m!" #&$%s, @?F@>@?A5@"U@@?F@>@?F@?FF?_UlvZ@>?x<@?FShoots-cell_size*(ic_mac_PhBio/sept_to_jan_conv) 8pzK>,O?~w? Bzuw  8 p   !GCarbon biomass of the photosynthetic portion of macrophytes (kg/cell)." #&$%c ?zG{?zG{???zG{?zG{?zG{?zG{?zG{?zG{?????zG{?zG{?x=?zG{Mac_PhBio_NC.01   ?z G{!guessed {N/C ratio}" #&$%c6 @SH`[W@S ?AvzW.@f@SH`[W OMtotBio7Phytoplankton+Shoots+RootsRhy+Epiphytes+Seaweeds+wrack 0<1|6 RY[V-, 0 <  1    |  6  !6Total biomass of all organic carbon food sources, kg." #&$%"f L?gi?=Kye??@!G{@?gi?=Kye?v?=Kye??????"Mm*x??=Kye PhBio_GPPMif(Shoots<0.001) then 0.0 else Mac_PP_rate*Mac_Prod_lim*Shoots^Alpha_PhBio<qzB<D < < ?PbM 񖩟         <    !{SHORT: Photosynthesis. NOTE: does not represent actual leaf growth! Actual leaf growth can be computed as the difference between leaf growth and root growth}" #&$%c  ?PbM񖩟@? @ Pot_Erosstress_effecti?zG{?zG{?Q?zG{??Q?Q?zG{? =p ??2Y׮??uMUi?DPM?H?ξFV?Ku%F?LX%[[?MSn?Mŵ ?Mŵ ?N% !" #&$%c(M ?GLk?G&X\???GLk?G&X\?G&X\?G&X\????GLk?A?G&X\Epi_Resp_Temp_lim)Epi_resp_rt*(0.16*exp(0.054*H20_Temp))U4  ( ?z G{ ( ?Ǵ S U      !" #&$%f .?@!náu?9xO@9????xB NPhBio_GPPPHMAC_to_Roots*PhBio_GPP`?qz<1  ` ?  !:{SHORT: Roots now grow from photosynthesis in leaves} 1.4" #&$%s ??4<????4<??4-z_>Nb34?> jA>^E7F]Z)>-z_>Nb34>wSc@>wSc@ ????>-z_?K>wSc@Cons_ingest_PhBio/Min(Shoots,(SH_pref/pref_tot)*Con_pot_ingest)<-Uu  <    -   ! units 1/h" #&$%c0>-z_>Nb34?>5= A@nA>-z_>Nb34>wSc@>wSc@ ????>-z_?L>wSc@Cons_ingest_NPhBio1Min(RootsRhy,(RR_pref/pref_tot)*Con_pot_ingest)1-U  1    -   ! units kg/h" #&$%c  ??S'?9??.q????S'?9?9?9?????zG{?S'?M?9Cons_activity H20_TempU-h@@@@"@'UUUUUV@,@0UUUUUV@2@5@7UUUUUU@9@;h?p =q?ˣ =p?GzH?Q?=p >? =p????zG{?p =q?Gz? U!" #&$%c D?@"@v????N Seed_germ0I ! was: If Days =100 THEN 1 ELSE 0" #&$%c'a ?333333?bf ??@$?333333?bf ?bf ?bf ,?????333333?xO?bf R_to_L(if Shoots=0 then 0 else RootsRhy/Shoots<1, <      1 <  !-SHORT: Ratio of root biomass to leaf biomass" #&$%c %r?3*0U2b?3*0U2b??zG|?zG|?zG|?zG|?3*0U2b?3*0U2b?3*0U2b?3*0U2bH??????3*0U2b?P?3*0U2b ic_epiphytes&0.1*(ic_mac_PhBio/sept_to_jan_conv)8p  ? 8 p   !IFor every 100 grams of eelgrass there is 10 grams of epiphytes {kg/m^2}" #&$%(fS?????RSW_reproductionTIf sw_repro_season<1 then 0 else If (Random(0,1) < 0.1) then sw_seeding_rt ELSE 0JS|\ J          < ?    S   !Behm: random seeding which can only occur during the reproductive season Reproductive season from: Niemeck, R.A., and A.C. Mathieson. 1976. An Ecological Study of Fucus spiralis L. J. exp. mar. Biol. Ecol. 24: 33-48." #&$%c~?zG{?zG{???zG{?zG{?zG{?zG{????zG{?S?zG{sw_seeding_rt0.01R   ?z G{!" #&$%c  n@<|?333333?GzH???333333?GzH?GzH?GzH-????333333?U?GzH H20_TempDaysoMAy\ p?@@@@@@@ @"@$@&@(@*@,@.@0@1@2@3@4@5@6@7@8@9@:@;@<@=@>@?@@@@@A@A@B@B@C@C@D@D@E@E@F@F@G@G@H@H@I@I@J@J@K@K@L@L@M@M@N@N@O@O@P@P@@P@P@Q@Q@@Q@Q@R@R@@R@R@S@S@@S@S@T@T@@T@T@U@U@@U@U@V@V@@V@V@W@W@@W@W@X@X@@X@X@Y@Y@@Y@Y@Z@Z@@Z@Z@[@[@@[@[@\@\@@\@\@]@]@@]@]@^@^@@^@^@_@_@@_@_@`@` @`@@``@`@`@`@`@a@a @a@@a`@a@a@a@a@b@b @b@@b`@b@b@b@b@c@c @c@@c`@c@c@c@c@d@d @d@@d`@d@d@d@d@e@e @e@@e`@e@e@e@e@f@f @f@@f`@f@f@f@f@g@g @g@@g`@g@g@g@g@h@h @h@@h`@h@h@h@h@i@i @i@@i`@i@i@i@i@j@j @j@@j`@j@j@j@j@k@k @k@@k`@k@k@k@k@l@l @l@@l`@l@l@l@l@m@m @m@@m`@m@m@m@m@n@n @n@@n`@n@n@n@n@o@o @o@@o`@o@o@o@o@p@p@p @p0@p@@pP@p`@pp@p@p@p@p@p@p@p@p@q@q@q @q0@q@@qP@q`@qp@q@q@q@q@q@q@q@q@r@r@r @r0@r@@rP@r`@rp@r@r@r@r@r@r@r@r@s@s@s @s0@s@@sP@s`@sp@s@s@s@s@s@s@s@s@t@t@t @t0@t@@tP@t`@tp@t@t@t@t@t@t@t@t@u@u@u @u0@u@@uP@u`@up@u@u@u@u@u@u@u@u@v@v@v @v0@v@@vP@v`@vp@v@v@v@v@v@v p?333333?GzH?\(\?? =p?zG?zG{?Gz?GzH?zG??=p =?zG?p =q?333333?\(?p =q??ffffff?333333??Gz?\(?ˣ =p??Q??Q?zG{?zG{?zG{?zG{?Q??Q??ˣ =p?\(?Gz??333333?ffffff??p =q?\(?333333?p =q?zG?=p =??zG?GzH?Gz?zG{?zG? =p??\(\?GzH?333333@\(@@ =p@\(@@@Q@QR@ Q@ Q@ @ GzH@ =p@Q@ =p@\(\@ =p @\(@zG{@=p =@\(@Q@Q@GzH@zG@Gz@zG@zG{@QR@Q@\(@ =p@zG@\(@ 333333@ =p@!@!=p @!G{@"p =q@"fffff@#W =p@#@$B\(@$p =@%333333@%\)@&# =p@&Q@'zG@'\(@( =p @(Q@)@)zG{@)\(@*p =q@*Q@+ffffff@+Gz@,\(\@, =p@-QR@-@.GzH@.\(@/=p =@/Q@0@0W =p@0zG@0\(@1 @1GzH@1Q@1@1G{@25\(@2p =q@2ǷQ@2fffff@3Q@3W =p@3ΑR@3=p @4@48Q@4p =q@4fffff@4(\@5R@5GzH@5zG{@5zG@5Gz@6zG@6GzH@6xQ@6\)@6@7zH@78Q@7ffffff@7ΑR@7@7Q@8 =p@8@@8kQ@8zG@8G{@8 =p@9 =p @9.zG@9TzG@9xQ@9@9p =@9Q@9p =@:(\@::G{@:Y@:u\(@:ΑR@:ǷQ@:Q@:Q@:\(@; @;# =p@;8Q@;L@;^Q@;p =q@;@;ΑR@;Q@;zG@;G{@;Q@;\(@;@;Gz@;\)@;ݣ =q@;\(@;Q@;p =@<@<@<@<@;p =@;Q@;\(@;ݣ =q@;\)@;Gz@;@;\(@;Q@;G{@;zG@;Q@;ΑR@;@;p =q@;^Q@;L@;8Q@;# =p@; @:\(@:Q@:Q@:ǷQ@:ΑR@:u\(@:Y@::G{@:(\@9p =@9Q@9p =@9@9xQ@9TzG@9.zG@9 =p @8 =p@8G{@8zG@8kQ@8@@8 =p@7Q@7@7ΑR@7ffffff@78Q@7zH@6@6\)@6xQ@6GzH@6zG@5Gz@5zG@5zG{@5GzH@5R@4(\@4fffff@4p =q@48Q@4@3=p @3ΑR@3W =p@3Q@2fffff@2ǷQ@2p =q@25\(@1G{@1@1Q@1GzH@1 @0\(@0zG@0W =p@0@/Q@/=p =@.\(@.GzH@-@-QR@, =p@,\(\@+Gz@+ffffff@*Q@*p =q@)\(@)zG{@)@(Q@( =p @'\(@'zG@&Q@&# =p@%\)@%333333@$p =@$B\(@#@#W =p@"fffff@"p =q@!G{@!=p @!@ =p@ 333333@\(@zG@ =p@\(@Q@QR@zG{@zG@Gz@zG@GzH@Q@Q@\(@=p =@zG{@\(@ =p @\(\@ =p@Q@ =p@ GzH@ @ Q@ Q@QR@Q@@@\(@ =p@@\(?333333 !.graph based on: 14-14*cos(2*pi*(Days-31)/365)" #&$%cw?tzG{?tzG{??tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?????tzG{?V?tzG{Epi_PC0.005    ?tz G{!guessed {Phos/carbon ratio}" #&$%cx?zG{?zG{??zG{?zG{?zG{?zG{?zG{?zG{?zG{?zG{?????zG{?zG{?W?zG{Epi_NC.01   ?z G{!guessed {N/C RATIO}" #&$%c2 ?zG{?zG{???zG{?zG{?zG{?zG{????zG{?jX?zG{ Salinity3if SF_water=0 then 0 else (SALT_SURF_WT/SF_water)  2,       2   !(BEHM) Converts salt mass (kg) to salinity (kg of salt per 1000 kgs water). Assuming a meter cubed of water weighs 1000 kg. kg/((1000 kg/m^3)*(m^3/1000))" #&$%c v?hI?hI???hI?hI?hI?hI;????hI?Y?hIseed_loss_rt5.8e-5D   ?h I!<0.0014 assumes 100% loss of a seeds after 2 years 0.0014/24" #&$%c;???sͼ?-V?????[S_Current_Vector0n !" #&$%c ??qOp????qOp?qOp?qOp?????\?qOpSw_Resp_Temp_lim0.16*exp(0.054*H20_Temp)Ux,  ?z G{ ( ?Ǵ S U     !" #&$%f ' ?????j^salt_downflow(if SF_water=0 then SALT_SURF_WT else 0  22$     2   !" #&$%c&f???Ra?/\)G?????x`PHMAC_to_Roots'(1-Leaf_Limit_Coeff)*leaf_grow_season~B  ~    !HSHORT: Computes net siphoning (as a percentage) of leaf growth to roots" #&$%c ????????i?????a?convert_currents1n>kl !jBEHM: Hydrology output from Safak's model is in m/s. We need m/hr. There are 3600 seconds in 1 hour." #&$%c& @Y ?.QS?&c??@tx(7?lPZ?.QS?&c?.QS?.QS????&c?-b?.QSh_angle'-TAN(Latitude*0.0174533)*TAN(Decline)fH????H@ "܌`A@Z1'@˰n?+ I^?"܌`A?+ J?˰n?"܌`A0  * ?ΧG "@_         !" #&$%f ?8c SW_Gross_PP SW_Prod_lim*sw_pp_rt*Seaweedsi|h| i   |  !91=alpha .03 = plankton prim. productivity rate, 1/d {kg}" #&$%f ,>ս݁>`D??3$DjÌ>/22?Tm>ս݁>`D>`fk9>`hi7 >ս݁d>`hi7 SW_consumedCons_ingest_SWE| E! {1/timestep}" #&$%f-??~ r`?~ e SW_X_outSW_X_E+SW_X_W+SW_X_S+SW_X_N |        !" #&$%,cZ @\?-4?Vg???߶&?۵Fޱd?-4?Vg?Vg?Vg ????-4?-f?VgZenith[(SIN(Latitude*0.0174533)*SIN(Decline)+COS(Latitude*0.0174533)*COS(Decline)*SIN(h_angle))bH????H@ !n@YE @To?,]N?!.H?CA?To? l  * ?ΧG "@_         * ?ΧG "@_        b     !ACOS function" #&$%Lc;p???ugȿqp_Op???p?-gSDCLN0.00678+0.39762*COS(0.0172142*(Days-172)) +0.00613*SIN(0.0172142*(Days-172))-0.00661*COS(0.034428*(Days-172))-0.00159*SIN(0.034428*(Days-172));;  ?{X dE# + ?r ( ( ? hrL        + ?y XdE$ ( ? hrL        - ?{ ]9 ( ? "$        - ?Z |[ ( ? "$        !" #&$%cbyb]ztF?3`Q}8Q}8byb]ztFbybyb]ztF?hbySW_NPP SW_Gross_PP-SW_respcx  c x  ! {kg/d}" #&$%c ? ?'e26+3?i SW_Prod_lim@if SF_water >0 then SW_Light_lim*SW_Temp_lim*SW_Nut_lim else 0 yEc4      y  E    !" #&$%cB ?@-q?By5@-q?????jk Salt_X_WCIF sf_wt_X_W > 0 THEN sf_wt_X_W*Salt_conc*convert_currents ELSE 0a4        a    !" #&$%cB] ??߃?????jl Salt_X_ECIF sf_wt_X_E > 0 THEN Salt_conc*sf_wt_X_E*convert_currents ELSE 0a4        a    !" #&$%c<???BY y|?????mW_Current_Vector0 !" #&$%"c]@?携˥q?携˥q???e=kD6T?携˥q?携˥q?携˥q?携˥q1????携˥q?n?携˥q S_vector^if S_Current_Vector<0 then 0 else (S_Current_Vector*convert_currents)+(0.035*NS_wind_vector)[a7D [      [ a  ( ? Q    !2BEHM: convert_currents converts from m/s to m/hr" #&$%*cL?ٕPҚh?^8??5}?'?pNtAۘ?'?ٕPҚh?^8?^8?^8?????ٕPҚh?xo?^8Mac_Temp_Resp_LimMif switch_temp=1 then 1-(40-Air_temp)/40 else 0.0107*exp(0.16*(H20_Temp))wUzd w      (   (     # ? ' ( ?z G{ U     !if switch_air_temp_lim=1 then 1-(40-Air_temp)/40 else .0107*exp(.047*(H2O_Temp)) {SHORT: Separate limitation function for respiration using water temp. Switch can go back and forth as needed.}" #&$%ca@4@4??@4@4@4@4???@4?xp@4sept_to_jan_conv20 4<P !" #&$%c2 ?{$B?sdb]ztF>by??Q}8?3`?Q}8>b]ztF>by>dP=>by%>b]ztFx>bySW_respBif SF_water >0 then sw_resp_rt*Sw_Resp_Temp_lim*Seaweeds else 0 \|h|4      \  |    !&1.1=alpha coeff t_Pl_Resp, 1/d {kg}" #&$%+c >?#???lu?=e]@cF??#????Lcl.?#?y? SW_Temp_lim?min(1, 0.22*exp((0.07-0.00001*exp(0.28*H20_Temp))*H20_Temp))Uih  , ?( \) ( ? Q - > Ëh ( ? Q U      U      !" #&$%c?e3S?8Ė??@qBH{q?e3S?8Ė?8Ė?8Ė?????<͵C?xz?8ĖPhBio_Net_Prod PhBio_GPP-NPhBio_GPP-PhBio_Resp?Bzv ? B  z  !SHORT: Net leaf production" #&$%c 7 ׌hS??d{Pq*׌hS׌h׌h?????ušS?x{׌h Root_depth8(1-L_to_R_limit)*RootsRhy/cell_size*Nph_dpth_to_densW1 $  W  1      !(1-L to R limit) = proportion of root biomass to total non photo biomass MAC_NOPH_BIOMAS/cell_size = nonphoto density, kg/m^2 2.5 = root depth/root density m/(kg/m^2) of plant type UNITS m" #&$%s!@@SJ@@@BlRؔ@S21!@BlRؔ@@SJ@SJ@SJ?.db@?|@SJ Seaweeds"water_depth*cell_size*ic_seaweed E>~cxcRexd      !Carbon biomass of seaweeds ic_Pht{mg/l}*.001{convert to g/l = kg/m^3} Behm: took out conversion, not necessary for seaweeds. Was: water_depth*cell_size*ic_seaweed{mg/l}*.001{convert to g/l = kg/m^3}" #&$%c  ?p`DP?oq???޲!?F\?p`DP?oq?p`DP?p`DP?????19(?oq?x}?p`DP mac_tot_densmac_tot_biom/cell_size      !kg/m2" #&$%c5 ?#*0U2a?"I:K??_e.?%n _+?_e.?#*0U2a?"I:K?"I:K?"I:K$?????#*0U2a?~?"I:KSW_conc6IF (water_depth > 0.01) then Seaweeds/SF_water else 0|  4  > ?z G{    |     !%BEHM: Before, SW_conc calculated like this: IF (on_map < 2) AND (water_depth > 0.01) then Seaweeds/SF_water else 0 WHY THE CHANGE? Took on_map considerations out of concentration calculation and put it into seaweed export so that seaweed export was calculated similarly to wrack export." #&$%cl???????????U???????veg_dens_coeff.1   ? !Vhabitat dependent variable, stems (or structural attribute) per kg of plant material" #&$%c??????????m??????x?switch_nut_lim1 !n{SHORT: 1=nutrient limitations are used in computing Prod Lim 0=no nutrient limitations are used}" #&$%f?nO0.75 then 1 else 0 ,  > ?       !SHORT: seaweed is moved only under high water velocity conditions. Need help calibrating this. {SHORT: generalizing algae model} Behm: Was modeled as graph where seaweed moved only when current velocity > 0.250" #&$%c @Dfffff@Dfffff??@Dfffff@Dfffff@Dfffff@Dfffff@Dfffff@Dfffff@C@@Dfffff?-@Dfffff Latitude 41.3 {degrees}bf   @Df ffff!Ninigret Pond, RI" #&$%"c]S?J*?J*???.2p3 P*?J*?J*?J*?J*k????J*??J* W_vector^if W_Current_Vector<0 then 0 else (W_Current_Vector*convert_currents)+(0.035*EW_wind_vector)ma"D m      m a  ( ? Q    !lWind drift factor is 3.5%. McDonald 1992 oilspill model. BEHM: convert_currents converts from m/s to m/hr" #&$%c??????????e??????x?switch_light_lim1 !fSHORT: 1=light limitations are used in computing Prod Lim 0=no light limitations are used" #&$%c????????i??????x?switch_temp_lim1 !jSHORT 1=temperature limitations are used in computing Prod Lim 0=no temp limitations are used " #&$%6c m%@` PSDD@Jۨf???@` PSDD@Jۨf@Jۨf@JۨfU????@Zp!߱?-@Jۨf PAR_Shootsnif Shoots=0 then 0 else if (1.15*Epiphytes)>Shoots then 0 else PAR_Epi*(1-(0.75*SQRT(1.15*Epiphytes/Shoots)))</%% <    %  ( ?ff ffff   <    %    ( ? ( ?ff ffff   <       !VSource: Buzzelli et al. 1999 (BEHM) PAR_Epi*(1-(0.75*SQRT(Epiphytes/(Shoots*2.9))))" #&$%c ?-F:?5S{???(~[|2M?-F:?5S{?-F:?-F:????5S{?-?-F:Percent_LightPAR_Shoots/PAR_Surface%   %  !" #&$%c)?Ÿ{?Ÿ{?`bM񖩟??Ÿ{?Ÿ{?Ÿ{?Ÿ{????Ÿ{?x?Ÿ{Shoot_Mort_Ratio*(PhBio_litterfall+PhBio_Mortality)/Shootsw< w   <  !" #&$%f>?^N۲o?]-V??@?^N۲o?]-V?^ s?^N۲o@?????]-Vx?^N۲oPhBio_MortalityPhBio_Mort_Rate*Shoots <<   <  !A{SHORT: Leaf mortality as a function of leaf weight. PhD Thesis}" #&$%cA?Ÿ{?Ÿ{???W縴?W縴?Ÿ{?Ÿ{?Ÿ{?Ÿ{?????Ÿ{?x?Ÿ{PhBio_Mort_Rate(0.005/24)*0.29  ( ?Ÿ {!0.005*leaf_weight*10 Base 0.005*leaf_weight*10*.1 {SHORT: Computation of leaf mortality based on leaf weight. Constant in this simulation since leaf weight is constant.}" #&$%cF@@??@@@@@@????@?x@Nph_dpth_to_dens2.5{   @ !" #&$%c HMYRN-G& ??@@MYRN-G& MYRMYR????N-G& ?xMYR NPhBio_NPPNPhBio_GPP-NPhBio_RespBv  B   !SHORT: Net root production" #&$%c@L@L??@L@L@L@L???@L?-@LIk_Epi57.5   @L !" #&$%c ?tzG{?tzG{???tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?????tzG{?tzG{?x?tzG{Mac_PhBio_PC0.005   ?tz G{! {Phos/carbon ratio}" #&$%c ?zG{?zG{???zG{?zG{?zG{?zG{?zG{?zG{?????zG{?zG{?x?zG{mac_nonph_NC.01   ?z G{!guessed. {N/C RATIO}" #&$%c &d?#*0U2b?#*0U2b??zG|?zG|?zG|?zG|?#*0U2b?#*0U2b?#*0U2b?#*0U2bP??????#*0U2b??#*0U2b ic_seaweed'(ic_mac_PhBio/sept_to_jan_conv)*0.058p| 8 p  * ?  !Qthere are 5 grams of seaweeds for every 100 grams of eelgrass {mg/l, general #}" #&$%ce?tzG{?tzG{??tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?tzG{?????tzG{??tzG{SW_PC0.005   ?tz G{!guessed {Phos/carbon ratio}" #&$%cf?zG{?zG{??zG{?zG{?zG{?zG{?zG{?zG{?zG{?zG{?????zG{?zG{??zG{SW_NC.01   ?z G{!guessed {N/C RATIO}" #&$%c @L@L??@L@L@L@L???@L?-@L Ik_Shoot57.5/   @L !" #&$%c5@h^TM!@h?@FS??@@h^TM!@h?@FS@h^TM!@h^TM!????@DI$I%@h?@FS?-@h^TM! sunlightSolar_radiation !" #&$%c  ?????????????y? PO4_scen1 !" #&$%c6 ?6C.?~^????6C.?~^??~^??~^?????6C.??~^? Epi_conc7IF (water_depth > 0.01) then Epiphytes/SF_water else 0 4  > ?z G{         !" #&$%c /?M}H?M}H???Q?Q?M}H?M}H?M}H?M}H??????M}H?x?M}Hic_mac_NPhBio0Perc_NPhBio*(ic_mac_PhBio/sept_to_jan_conv) 8p1  8 p   ! {kg/m2} " #&$%c ?333333?333333???333333?333333?333333?333333x????333333?x?333333 Perc_NPhBio0.30   ?33 3333!y30% of total Eelgrass Biomass is stored in the roots and rhizomes (estimate can be substantiated with the mesocosm data." #&$%c \ ?؜y?؜y??؜y?؜y??؜y?؜y?؜y?؜y?qq?؜y?x?؜y Mac_Nut_lim]min( DIN_sed_conc/ (DIN_sed_conc+Mac_Sat_N), PO4_sed_conc/ (PO4_sed_conc+Mac_Sat_P) ) )0 ) )          !Mich-Menton, assumes only one limiting nutrient. half sat coeffs are broad approximations (2.5E-5g/l=25 g/l) {g/l SHORT: Need to add Surface H2O concentrations and limitations} " #&$%c >￵Ëh>￵Ëh??>￵Ëh>￵Ëh>￵Ëh>￵Ëh???>￵Ëh?x>￵Ëh Mac_Sat_N5.0E-6   >￵ Ëh!" #&$%c >Ëh>Ëh??>Ëh>Ëh>Ëh>Ëh???>Ëh?x>Ëh Mac_Sat_P2.5E-6   > Ëh!" #&$%c?a 8C?a 8C???a 8C?a 8C?a 8C?a 8C6????a 8C?x?a 8Cmax_litterfall_rate0.00208w   ?a  8C!70.05 is max spec rate of litterfall, 1/d 0.00208, 1/hr" #&$%c W? d? 1N??? d? 1N? 1N? 1N????ըX?-? 1NEpi_Light_limPAR_Epi/(PAR_Epi+Ik_Epi)      !" #&$%c ?ˇ ?DšY???ˇ ?DšY?DšY?DšY????D^oy?-?DšYPht_Light_limPAR_Pht/(PAR_Pht+Ik_Pht)      !" #&$%c ?ˇ ?DšY???ˇ ?DšY?DšY?DšY????D^oy?-?DšY SW_Light_limPAR_Sw/(PAR_Sw+Ik_SW)3%i 3 3 %   !" #&$%c@a@a??@a@a@a@a???@a?-@aIk_Pht140 !" #&$%c ?????????????? Alpha_Pht1 !" #&$%c?Q?Q???Q?Q?Q?Q????Q?x?QPI_trans_rate0.03   ?Q !" #&$%cF: ?$ʹz?$ʹz???+׼!?$ʹz?$ʹz?$ʹz?$ʹz????$ʹz??$ʹzS_wrack_vectorGIF S_vector<0 THEN 0 ELSE S_vector+(Wind_Wrack_factor*NS_wind_vector)n(4 n      n     !" #&$%c ?J*?J*???J*?J*?J*?J*????J*??J* EW_vectorE_vector+W_vector      !" #&$%c ?????????{???????-? leaf_length1 !|PTS Thesis{SHORT: Leaf length as computed in Thesis} leaves*Shoots/(leaf_weight*leaf_density) 5*Shoots/(.29*500*cell_size)" #&$%c ?????????????y? DIN_scen1- !" #&$%f ?????j River_SaltRiver_NACL*SF_water 2     !Input of salt from Little Bay" #&$%f ?????  River_SedRiver_TSS*SF_water*River_CF .3    .  !gm/l*hr = kg/m3*hr *m3 = kg/hr" #&$%  @"    &&&&&&&&&&&&&&&a D@k$'() Graph 2*:+ !&&DEؕ@'R'() Graph 1*2+ !&&&&&&&&&&&&P@1 @>f8F KZ[\#'( $Phosphate)*+ !&PD@DJt<R'($DIN)*+ !&DN+@CbCPAl>'() Graph 3*D+ !&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&C$#o h K9D' << ()*+ !&&&&&&&&&&&&&&&Cp' )iZXH'`@`()*+ !&C$#o i[9D' <<()*+ !&C}(n )k\H'mm@()*+ !&&&&&&&&&&&&&&D$@Sؑj7'() Graph 6*s+ !&&&&&&&&&&&&&&&&&&&&D)@؇~ mG'() Graph 11*+ !&&&D@إ'() Graph 9*+ !&D@؝  |'() Graph 10*+ !&&&&&&D*T@ >Hj'() Graph 13*+ !&&DB@m 5z`'() Graph 14*+ !&&&&&&&&&&&&&&&&&&&Cv L`_ 'g()*+ !&&&&&&&&&D@ .'() Graph 17*+ !&&&D @%  > '() Graph 15*+ !&D@%  '() Graph 16*+ !&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&P 42 4<<,DDE'($TSS)*+ !&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&CP JRE8'@@()*+ !&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&C&n ]dh'pp()*+ !&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&FD'#,YH'44@()*+ !&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&C}(n .H'm@m()*+ !&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&P0 ~ Zd#'( seaweed$)*+ !&&Cv  Xc_5'g()*+ !&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&P;  ;g oE'( Consumers$)*+ !&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&P  n  \d'( eelgrass$)*+ !&P   @ oHK['(Phytoplankton$)*+ !&&PD D<Tq<D\'($Salt)*+ !&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&Czhn _P}'gXX()*+ !&&&&&&&&&&&&&&&&&&&&&C/ o Z_!'g()*+ !&&Czhn  p_P}'XgX()*+ !&&C/ o Z_!'g()*+ !drmd 1 '  nb   1&aL@ i*b@ʧ`ZZ((/s'()*+ p&c w'()*+ p&c4X؝ؕ ا.As |}'()*+ p&a;ب4@Te@%N'()*+ p&c|@f'()*+ p&cؕpw'()*+ p&cjؔ؝ ]Xw'()*+ p&c*.@/##=.vE'()*+ p&c*؟?؟ ?79 '()*+ p&c!nؕؖؖ@ *2 GHQx'()*+ p&coaPAb J'()*+ p&cؚؔ y Bxq'()*+ p&f @C ~6' ()'n +n*+ p&a7 ؚy@0"Dt['()*+ p&cg @&WCEZ1f'()*+ p&aI-/@FD'()*+ p&a8;@B@0 G,'()*+ p&a)@@0 H'()*+ p&fGR؟I2ؒ2 I,a5'GG()'QQ *+ p&cf ؔ))J0s -6'()*+ p&c$2ؔ؞@ K1?P'()*+ p&c6Lp'()*+ p&czlA@32767.5ZMm^'()*+ p&c ؕ ؊ ،N'()*+ p&aG3 h,,mmO%%B'()*+ p&a IKؼI@ P2O'()*+ p&aT @)@0 Q$Y'()*+ p&c؟M؝ R  N u'()*+ p&cebitat orSX kl'()*+ p&c^؟ P TQ'()*+ p&f3U2ؿ؏2 Ul^DLOy'()'4\4\*+ p&aK- `i!!V@@2'()*+ p&c/8؜ؕWvI]d'()*+ p&fXؘؖ\v XAc \  rj'XX()'oLL*+ p&c}@(YpZ'()*+ p&a6Yd?I@ Z4}'()*+ p&cDؔ@;9[wQ'()*+ p&al {A#@ `\DX-'()*+ p&c] ]9Q'()*+ p&cG^ ^:T8>'()*+ p&c9_ _,W'()*+ p&cCC` `66aP'()*+ p&cZa aMx'()*+ p&cshb bf[u'()*+ p&cWؔ؝pceu'()*+ p&cid؝< d> }'()*+ p&acx' Z     e!!'()*+ p&a1Ej؄5@ f>fF'()*+ p&d^@I?g}Ug'()*+ p&fh؏؞& h?C #O<'G()'*+ p&f/iؖQk i :z'//r()'xx*+ p&c8F7ؔ؟E_ j+ ES@f'()*+ p&av<Sg@0 kb{'()*+ p&aIS`@0 lDe'()*+ p&fj2ؘ2 m(  /M' ()'{ { *+ p&cؕؖؖn. %S+'()*+ p&c @]*@ Jo'()*+ p&a} bUUp++e'()*+ p&aBI 66q&&'()*+ p&a_7؉@0"r~'()*+ p&fM~s2؏2 s%6! tP'>M>M()'p   *+ p&a?-J@y_t:O'()*+ p&a%cxض8 آآئئu>>'()*+ p&c >|x@kBv # Y'()*+ p&a?x+@zw&D'()*+ p&gne n nSxQ v'( $Eelgrass)*+ p&c 1zx@(|/y  Y'()*+ p&cW}ؔؿ zJMdy'()*+ p&cG؜`{$^'()*+ p&aL{-؆) ؠp|;M]Z'()*+ p&d'x }0'()*+ p&fu2؞ؖ2 ~\ %8'd ()' -*+ p&c; Qt%Q8'()*+ p&cpl- ؕ@\s;`R} '()*+ p&cvؔ؟ '()*+ p&aek2F| ߮R||؀؀00`Ap'()*+ p&cfjSP|'()*+ p&cUF؟ بVH9l'()*+ p&c) ؕDPTor'()*+ p&cؘ؎ ؛r '()*+ p&ajMjk@gؑؑ##;;8g'()*+ p&aHx 4kNJJ%% '()*+ p&fw2 2 \ R>}'d()' K w w*+ p&dj1 @a(s:'()*+ p&a-/@,'()*+ p&c؟;؞I c @ 0'()*+ p&f z2ؖ2 D = b7'<()' - -_*+ p&a` yx@Np *'()*+ p&c|@}Q '()*+ p&c%~ؙؔMBg *2ehe'()*+ p&a{xث A@:mHq'()*+ p&cHؕ 0j WXAB'()*+ p&d 7P  '()*+ p&cO f  Bm'()*+ p&a-/ *@'()*+ p&cIؔpSip V'()*+ p&cTؘ  LP/ /  4a '()*+ p&c-?<p? K  j'()*+ p&c-B ؕ :'()*+ p&a+jk'@ W| J6'()*+ p&dh }0y'()*+ p&a /-"G@'()*+ p&a{-E@  J2'()*+ p&c&4@@O1'()*+ p&a?'hiww__N'()*+ p&ab y-#@ Lp'()*+ p&&ac-<@/@7^'()*+ p&&ciY \7{'()*+ p&d7i jш.@x'()*+ p&&aUSG9q@0 `'()*+ p&c{ؕ ؊ ،] kS'()*+ p&&cEؕ 'R V 29'()*+ p&&&&&&&fr2ؖ2 +1 r*'3()'! G! *+ p&f؟2 ؓ2 M '<()'   6*+ p&&&cP ؛ؕ W2] '()*+ p&c؟,0  ^:Br'()*+ p&c~ؕ q'()*+ p&f؟2ؖ2 `  5uvb'h ()'  *+ p&&cMn؟ؘ ؛>aZ('()*+ p&ak ;L!ph0؛؛ئئV'()*+ p&c  v,'()*+ p&cuؠ h%'()*+ p&cؘ؏ l '()*+ p&&&c:ؔ؛,Xh*G*'()*+ p&cؕ ^ /Iu '()*+ p&fP2ؙؐ2 o E 9'()'  *+ p&alymR , iغغ__UUN'()*+ p&&c}5 pK Z>'()*+ p&f+2ؗ2 t 'l()'> *+ p&f02ؘؒ2 N O P ltxU'V0V0 G()'   *+ p&&a\-=m@4>hPI}'()*+ p&a-Jk@4 f'()*+ p&c^  Q| ,l'()*+ p&ceW %ؔ؜/ ?.ز/.X rd |4j:'()*+ p&agB[{`f@ U#?'()*+ p&cXN Kv!R'()*+ p&c   *Iq'()*+ p&aIUk@ :T[[cc==6 '()*+ p&a3 EQZZ66\\7'()*+ p&c0؟| V#Ko'()*+ p&c-tCtD{A'()*+ p&cCؔ 6jQew'()*+ p&a Q@ V'()*+ p&&aqX-Jd@}lSi'()*+ p&c=2-؛ؕ@<LՔNtF'()*+ p&abxب, iXX؆؆ss?6'()*+ p&a؈dVد @ b'()*+ p&ciLt؟ؕY K/va'()*+ p&d {   '()*+ p&c   >'()*+ p&c     '()*+ p&ad x@65_'()*+ p&fJ2؏2 "a4.   O',J,Ji()'p   ?*+ p&ctؔ (IQ'()*+ p&s4؟ؑ y    '()*+ p&cؔ ؛*+r'()*+ p&a~z?؆@0 Y'()*+ p&&aoN ^@ Yt'()*+ p&aR8dd@J\M_i'()*+ p&aoMi@iL{'()*+ p&&au/& `h-022JY?'()*+ p&cY Lw'()*+ p&c '()*+ p&cn% a2'()*+ p&c7ؙؔ *D'()*+ p&&a"@z ''()*+ p&czؕ ؚ\ ; '()*+ p&fd؟ؖg8H m8OY3:'dd()'t*+ p&c0Yؕii;=s'()*+ p&c, J'()*+ p&at-|@4>uo'()*+ p&a(-aؠ/@JcoL\'()*+ p&aa-n؆@Jq7Sf'()*+ p&&&&c` !@&.,.S~9'()*+ p&c;Z@?瓮B#S'()*+ p&s+@& u= '()*+ p&a`T-N 0؜ZZ== Oe'()*+ p&cxx؊ؕ SZ8 '()*+ p&cOT BmR'()*+ p&aq *( @ R'()*+ p&a-NC@ZRR44E"o'()*+ p&c S@B4 = 0'()*+ p&&a$n؀0: i '()*+ p&&aJ j]@lrثثAA8xR'()*+ p&c ؟ ؜  *f'()*+ p&&cY؟,0  L %(<zDQ'()*+ p&cW؟c}  J H!'()*+ p&a {Z@ & S7'()*+ p&&cM؞ @k\'()*+ p&a-N$v p؜??""].'()*+ p&a?{5@ *@J:'()*+ p&c\ؕOKOi'()*+ p&&aIL Z  PPQ'()*+ p&c2 ؕ l?:B'()*+ p&aZT V qؓؓغغQ_'()*+ p&cؕؒ   '()*+ p&c'  G<B'()*+ p&c{hؕ ]O (:'()*+ p&&a?-N؍8C $$))c '()*+ p&a~~"$@0 @'()*+ p&c:xؔ؜@TZJGh'()*+ p&a|Y-{R +؜M'()*+ p&c&R؟2 SC`E@ !'()*+ p&aHI  ZZUU88 VV$M'()*+ p&aR U ZZSS!MZ'()*+ p&c$|xؔ@Tir "1/?'()*+ p&a-U-v NX/؅؅XX#.'()*+ p&&aN n% xPiZZDD% S'()*+ p&a\-'= ؜oo&))'()*+ p&c؛;A[;A\'&'()*+ p&ae [ @g11(Vj'()*+ p&c@[0 J)'()*+ p&c 9@ I@*  %,'()*+ p&s؟+؟ + #&'()*+ p&a*x9&8D`33,66 '()*+ p&a[J-W5@4'-0`'()*+ p&&a."x[U@/)Z'()*+ p&s03ؓ 0 '()*+ p&aYWؠ@Za1.'()*+ p&s{2 ؓ 2ZZ'()*+ p&s83؟ؓ 3nIO'()*+ p&c"@Euvp4'()*+ p&a4T@0 5O'()*+ p&s\6؍ 6LmB'()*+ p&f؟ؔ-Pp/7@ x Oq' H()'h h *+ p&cQ؟Zk 8Do9'()*+ p&ap8. 0h-0؉؉>>944Rk'()*+ p&a:5 0i:|m'()*+ p&c; ; 7 >w'()*+ p&a# } {Pi<SS^'()*+ p&a*xf%@zx=*'()*+ p&crؔ؛x >eSD'()*+ p&a4"x%@i"?)c'()*+ p&cz؟@ @؛i#+'()*+ p&cQ|؟ؕ@؀iAc2gQ'()*+ p&a (iصصBKK'()*+ p&cؔ8@C}\*'()*+ p&c>x@tDh:'()*+ p&a QN'5@ E ky'()*+ p&c2ؕ*U FQ'()*+ p&cs G؛P'()*+ p&aK I dH P'()*+ p&cqؕؒ, I qH '()*+ p&ahAR$aPvvJMm'()*+ p&c3bitat orK&Q ^Vf'()*+ p&cs  L. r'()*+ p&&cX؝ؕeN:e'()*+ p&&aNg@ P}'()*+ p&aF]؞.  ZZثثQ(('()*+ p&aj(" ((R""CU'()*+ p&cPnؕQ QS2]'()*+ p&cYؕ@\T;fS'()*+ p&c  ؕ,D@U]#Z'()*+ p&c%@ VpV"'()*+ p&c*@%B`o[PW C-'()*+ p&a-Yؿ ?@8CXPwwX))v'()*+ p&cb @Yw ZX'()*+ p&aVU-Y؍؂ S@@؞eeZ~~ '()*+ p&c}]ؒ [ije/'()*+ p&cf؜\ ]&'()*+ p&aU-\ +؞SS]'()*+ p&anKV*O  i22rr^((.'()*+ p&c x@<Cm_b'()*+ p&c2؜*V`c'()*+ p&c`ؔ؟ am'()*+ p&a_x@Cb'()*+ p&a`zZZnncs'()*+ p&aR?{5g@ *@dM:l'()*+ p&aJ>{'w.@ *@e)9O|'()*+ p&c{ ؟KؗQ Qfm E'()*+ p&cmؕؕ gLw'()*+ p&cg ؞-;hq'()*+ p&al@0 i g'()*+ p&a2XU & aj::-Z'()*+ p&a@yk'()*+ p&cؕ lz'()*+ p&c>ؔtC@mKop'()*+ p&aDg@0 nbI'()*+ p&a, ma Pببo7'()*+ p&aIm 1a 0oop''Df'()*+ p&ah _@"q'()*+ p&a V'pʁhh,,r11!'()*+ p&f)؟F ؓ s 4& t'))()'r*+ p&d-s s1 t$6'()*+ p&c؟ؓ u'()*+ p&cjؔ v3w'()*+ p&ag; ZZw'()*+ p&a-@V @0 x(%['()*+ p&aMS*O@0 yH%TX'()*+ p&cW ؕQ R\z9dC'()*+ p&c#ؔrf {1'()*+ p&cV+rf |It:)'()*+ p&ac-؁0{ 8CXPSS++}qq2uQ'()*+ p&&f S2'؎2 g%&P'  ()' ((*+ p&f N5h ?  & kN'GG ()'  *+ p&d %   '()*+ p&a\-Nb *`]'()*+ p&c,1ؔ؝rf p9'()*+ p&al; ۈZZ'()*+ p&au;8 ZZ99'()*+ p&aE9{!a@ *@4Jf'()*+ p&fUG } GHJK'()'mm*+ p&d '()*+ p&avv;  ZZؿؿ<<CCq'()*+ p&a]#: ]]&&771O'()*+ p&a|; .ZZAA<<w"'()*+ p&a1 r/a ^^VV<<4'()*+ p&aph 22l'()*+ p&c9dؕR@9F'()*+ p&aa ~ \'()*+ p&c؟؞Q R Qj'()*+ p&a O+ص@Bh'()*+ p&c:*kؔؕ<Vd G7>~'()*+ p&a:ظ@B6S'()*+ p&c؟"?< L '()*+ p&d'()*+ p&s+ؔ..0{'('()*+ p&aC6{dd@ *@1Hi'()*+ p&aC2{#@ *@H7'()*+ p&c-9ؕE_ Q'()*+ p&ca:ؔE_ Tr2'()*+ p&c,CE_ J'()*+ p&cDؕE_ '()*+ p&d?6H'()*+ p&f5؎?# '7()'*+ p&c ؕ@-'()*+ p&aN-''g6 HH\'()*+ p&&aBD VcpJJ'()*+ p&ax:؁m @7`dd||}d'()*+ p&a@g|  cp؞؞##T7'()*+ p&aXSnJi؋؋wwSzs'()*+ p&cbpؔ؞: op'()*+ p&cVH6I}4'()*+ p&aY-\ x0؟bb'()*+ p&cQ!ؕxj%>a'()*+ p&c%J/NNC'5'()*+ p&cK P h P(7'()*+ p&a؞ع" cp؁؁%%'()*+ p&a~ cp~~5'()*+ p&aFY-ظE / TGG55/Q'()*+ p&&aay؍#@ ҈cp؏؏$$>>Y,'()*+ p&a$6C|{{,,{{R'()*+ p&c;-a{\ (PW'()*+ p&&c؟Rؑwg'()*+ p&a1z؟4R֥immMM^^Q'wr'()*+ p&afl&iغغ-- '()*+ p&ac؂i '()*+ p&aZg+x @BU}'()*+ p&a6t ӈcprr#rF'()*+ p&&cGؕ ~(dc]'()*+ p&aQ6+s@B Vx'()*+ p&c(ؕzCzDs'()*+ p&c -@ (Py'()*+ p&&&a(W?$ 4cp''()*+ p&a- ' Zcp؇؇CC"2'()*+ p&atؠ[ill~v'()*+ p&c ؚؕy   '()*+ p&aO'lvi'()*+ p&a*@-+)@4$&/E'()*+ p&c ؕ m89:'()*+ p&c1Sؗxj%@A^'()*+ p&cYw@e JL`w'()*+ p&c,SXxj%8Jn2\'()*+ p&c p ؛L>'()*+ p&a$[X~Q %@aDDBBVV'()*+ p&&cJY=qh01!'()*+ p&cZxj%_'()*+ p&a+xH @ZX` 0 '()*+ p&ca(/ؔ؞ Tn5'()*+ p&aU-M/ {p0--%%\\ 4R'()*+ p&c1gP$O#'()*+ p&d  "'()*+ p&c8[ؕE)-'()*+ p&a+ @B|'()*+ p&aRmjP>@ z9KWr'()*+ p&aMx()  c==$#R'()*+ p&sT؟1؎ ?^tu&cW3'()*+ p&a@؃شk  cp66OOhh)N'()*+ p&fRؗؓ@ة,, jY1Gg1'f)()'*+ p&azBr6I cpeeRR44'()*+ p&aJ؆ا@BT'()*+ p&aD+ (0cpVV ,'()*+ p&a+S?9 cp``  P'()*+ p&c$Vؙ؍ 1%,'()*+ p&ag_@BGW'()*+ p&aM 6I @BRN'()*+ p&a+F, cpLL**7\'()*+ p&aU:@BP'()*+ p&cFؿ ,'()*+ p&cH)Lmr2;fB'()*+ p&c!OHhHn?'()*+ p&c@P/ p?(*$_'()*+ p&a  icp'()*+ p&a3 aJ gcp__))JJA0'()*+ p&a  n cp '()*+ p&aV #اk cp:'()*+ p&a_^`?آ3 ??آآ33'>C'()*+ p&a]_y yy 6'()*+ p&aBf^ؾ` ؾؾ``0A'()*+ p&a5] e= ee==R'()*+ p&a[_a؇" ؇؇""<T'()*+ p&aab. ..Qko'()*+ p&agE`bI# II##AJlm'()*+ p&anra<؆< <<؆؆<<3U'()*+ p&ahmb|> uuBB9hm'()*+ p&aNQ 0cpS'()*+ p&a،6 ،،66Y'()*+ p&aA d hh<i$'()*+ p&awcG؄" ((qv'()*+ p&aIEL i+EE&&MM '()*+ p&aTR BKKW'()*+ p&av q!'()*+ p&c]3MH P{'()*+ p&aG1 1KTGG11^'()*+ p&a4 44q*'()*+ p&a~I!, 8R'()*+ p&agO pcp++88b '()*+ p&c؟؞[  '()*+ p&cd3؟f؝@ef JW&qh.'()*+ p&c[u0  3'()*+ p&c؟"?< Z js'()*+ p&a@o cp++W6'()*+ p&aB Lcp'()*+ p&aؑ T cpؓؓ!!RR*'()*+ p&aRQJ eTT '()*+ p&sP <ؓ  ;f_b]`@T|'()*+ p&amW"" FF%'()*+ p&s&ؓ@UX^  7#$%'()*+ p&aH1Ktjj C'()*+ p&a _c)FF؁؁xxM&'()*+ p&a@B '()*+ p&a<x[V@zx7['()*+ p&amn%'()*+ p&cx K'()*+ p&f(?ؘؒ/ 3F v}\`wZ'((>>()' $  K*+ p&&at- cp صص\'()*+ p&a   @B  '()*+ p&cC2ؔ@?yP'()*+ p&c؟4؝@?x'()*+ p&aHY@8`C^'()*+ p&c2 ؖI]<'()*+ p&&ax MMHH '()*+ p&aX[؉@ '()*+ p&gddؚرG''($Phytoplankton)VV[;Cr֏CCsr|7,o?-o? ?o,o??!72!7o?P, Geneva .(yalgal light limPh(Algal nut limitP( Alg Gross PP"b"bp"""Pz(my Alg consumed"8": u8u w:wpp|upv|w"z"u"sP(q Algae X out D Dp  """33PxJb(? t Al MortPl(b t Pl RespP( Alg Net prodP[is(u%Algal prod limP(zAlg resp"&"( &h (hp}bhhhb}hh"""P %(4incident lightPBZ( A Alg Temp lim UUUUP (SIS conc 0E("ALGAE UUUUPl+@ Algal concP(DIN sf wt conc33 PF^(m}ic algaePGW_o(nValg PCP $ + Algal conc UUUU0/(> SURFACE WATERP2J(, cell size P-=EU(T1 Algae mort :y: 8y8py3?y:y?9y3y8"9="98"66`)3`.8؉`.8<`)3 UUUUP]u+0 Sf wt E inPo(X Sf wt E outff PWo+ Algae conc EP0H((W Algae conc W UUUUP( Sf wt W outPj %+[ Sf wt W in ``؉`#<`# UUUUPg(`PO4 sf wt concP\'t(6L mac tot biom P 8(G Algae X in B D"B"Dp=IB=CID" "" UUUUP6(6~(Esunlight ``؉`<`` `؉`<` `{R\`{Wa؉`Wa<`R\`fp`fp؉`ku<`ku33PH+`C(o*alg NC UUUUP4L+MH2O TempP%=(LSOS concP!9( Sf wt depthP(L of Day PLd+c turbidityP>V(#tot max mac biom UUUUPM@eX(FCons ingest algae P74OL(^/alg code UUUUP( alg codeP(Current vel ave P1+k pKakjyx2P".3aZxh2P+1`.e'a "2P*0`2؎(a2P`pa:2P `%pqa2P "}a2P""xSa>CLQ2PuP{V`Y :a;2P`؞a2P`nNnax2P`maiw,2Pms"a2P`aL25a&C4Q\2P9?`0Za/2P-3`S8؍%a2P`رar2P`X ؉a2P"^}ba2P[za"oa`zn2Plr`6g{a_cmqض2P CI`7m؀a>L2P8>`,=a+2P"("a] k2P",a\j2P)/`a2P La2PIO""a2P%"(a"2P%+`\+na'g2Po%u"UV,a 2PRX`#_،ma7E2P17`ء9|a{2P_ReX`,bva&42PX^`&J4Xدa"052P3 9`(:3ؠa)7Q2Pio`._ػa(6J2PY_`+_La#12P"슚a2P`+Ra!/b2P`Kia12P"/a,2P".a\j2P`5ySLa\j2Pu{`>a:2P"a2P} o?2!378"""$2D~ABABABABA~Aإ""!:1Z?""!14p00>>p",",,"(,Algae0!g7./*+ p&a M[G@ ֈ/'()*+ p&aSz N'()*+ p&a[t@ ֈ!UF'()*+ p&a@ "'()*+ p&a-  S@ X^ #(X'()*+ p&a5 g@ X^ $0'()*+ p&a !>@ X^ %a'()*+ p&a@x؄%   $wDD& E'()*+ p&aص  ''()*+ p&at Acphh( '()*+ p&a^|?+ fHmЮ?<??)**.P'()*+ p&a HDD*'()*+ p&a @Y!!+x'()*+ p&asQm\ A cpxx,NN*w'()*+ p&a>Qز J cp؀؀->\'()*+ p&fEBؖؒ- .P #sP'EE()' K *+ p&amMزAMM؄؄/FF'()*+ p&aVQN fHmЮ?<bb0I['()*+ p&a{yM Gcptt1H ~'()*+ p&aX H0 Ecp؁؁&&2>>+]M'()*+ p&cBؕAf3G'()*+ p&cؕAf4&H'()*+ p&c>ؕAf5 K0K'()*+ p&&c^؟؜Af7Q kLJ'()*+ p&ae^!= !!8==Ny'()*+ p&a] 9'()*+ p&a }ؘi" ؘؘii:""V2z'()*+ p&&a hjX0 ZZ$$<v'()*+ p&&a;^< 55>A'()*+ p&cv@ l P/(?Z'()*+ p&c3j؛ؕ&.@g@A'()*+ p&a5@4 cpررA/:'()*+ p&c\ؕ RBt3'()*+ p&a-|78(Z||++C++i'()*+ p&a{-c(TpZbbDv'()*+ p&a-/JZ//ئئE33v'()*+ p&aQZ @F,sN'()*+ p&ajV3؇ f؄؄GffS'()*+ p&a4M +PؾؾH@@R'()*+ p&a~؆7??I11w'()*+ p&a7~ +P؆؆JddJ_'()*+ p&a5F +PظظK224U'()*+ p&&dm M'()*+ p&ssؓ<NRPE'()*+ p&dGa O>XPj'()*+ p&dPs PGY'()*+ p&auA$@VQLo'()*+ p&cY؟؜tCtDRwL'()*+ p&&clؔؿd6 T+y'()*+ p&&&&&&a xbB@z8Z''()*+ p&&aEx. Pkk\)J'()*+ p&aP xy\ iLYY22]6G'()*+ p&s xؓ ^#'()*+ p&sC x _3dde'()*+ p&a xdCR sXXLL`rr )Q'()*+ p&cO؟ at'cd'()*+ p&a+ x;) ##b==A'()*+ p&a@xa$ cIIc33E'()*+ p&aR_xa{ oH؜؜^^dM'()*+ p&a=_x$ |WW__eB '()*+ p&c0`ؔ3M@ f=hs='()*+ p&a[Y"@ ֈgv'()*+ p&axf؃ d~Oثثh77'7'()*+ p&&g  j'($Salt)*+ p&apZku'()*+ p&aUx5 NN22l0Z'()*+ p&aL30ZRR؀؀m v'()*+ p&&cؕ;UX o'()*+ p&&&&a5fx.B M gyrrs\\0G'()*+ p&a x* eet\\0W'()*+ p&ahx ppuPPc'()*+ p&axQ pp))v'()*+ p&cQ| -@L?w;g'()*+ p&a &3x" !!xzz-K <'()*+ p&g#  y'( $Nutrients)*+ p&&cO3 -%@Y,{B&mx\ de'()*+ p&cL,(U ؕ@@x|.YQ"Av'()*+ p&cd@}$ L'()*+ p&c!@~~'()*+ p&ajM؀2ZRRv'()*+ p&a!^x; <<@'()*+ p&ax.#@ PUUdd(M'()*+ p&cuEؕ mh  b(*'()*+ p&f4c2ؗؑ2 k?  8Mt'44()' O*+ p&c$ ~L@0'()*+ p&c{ ؟Va'()*+ p&c V'()*+ p&c|  !'()*+ p&c-  !G'()*+ p&axw5S ?tt '()*+ p&cdo A  6# h?'()*+ p&cPxt? Cn'()*+ p&a CZZ77؃؃]],'()*+ p&a +0Z'()*+ p&cJ؟ؙtCtD=vW'()*+ p&cgt؟tCtDZgt'()*+ p&f؟ ؓ c  56'  k ()'{z{*+ p&fؘؒ* ?A$ 'G()'1q1*+ p&c3؈ؕ@@7Eow@'()*+ p&c'؜@'()*+ p&a _؂!oll'()*+ p&g)/ )N ? V G'( $Epiphytes)*+ p&a ! Z""&'()*+ p&c5@D'()*+ p&cX@2O'()*+ p&a~U*@ /'()*+ p&aU+@ 0'()*+ p&a?/c@O:h'()*+ p&c6 ؜ؿ %4G='()*+ p&a R؊4v ;Zؕؕ&&qS'()*+ p&aEx6 ";;331J'()*+ p&&cO  Bm.'()*+ p&&a hi gZhhdm'()*+ p&c4Nhh eg'()*+ p&cIOi& 3_'()*+ p&cv5xO i'()*+ p&a xN ZxxI}'()*+ p&cY;ؗؕ@ ;fH'()*+ p&aA s$~LZhh^^((--'()*+ p&cVvؕ@ 8^c'()*+ p&a"R si4iw`Zii%%oo-U'()*+ p&ah sGPZ..cc,lm'()*+ p&a`T'(/-p&q'()*+ p&amo\ث@0g.~`\\ةة//y'()*+ p&cSؕ@ 5`'()*+ p&ab sAV`Z))vv+U{'()*+ p&cZ >ؐؕ<Vd <g-s'()*+ p&cGؿؕ=WW , T'()*+ p&am(jjأvZeeةةX&v>'()*+ p&c)R ؕؕHmn_Hn!6_F'()*+ p&a7 jP$Z @{{zz9za'()*+ p&aj؁@raå22I/'()*+ p&al]jwز>fHmЮ?<77=8'()*+ p&cj؟ؿ ]wY'()*+ p&cK؟؜ >XG'()*+ p&aOj]Z aågg*0G'()*+ p&a2jO؃k`aå??,(`'()*+ p&&&&&&c)@؛ؕ@B4ȁ 6K'()*+ p&&&&a6w LSZcc4S'()*+ p&cS؜ؕ@ 5`'()*+ p&aY Rmأ(e@Zkkجج_'()*+ p&a: r{#Zrr{{M'()*+ p&ad P2 Zee؂؂KKhy'()*+ p&d[ X pakRd L'()*+ p&c؟؝N"  '()*+ p&aka daRRg'()*+ p&a cauu;;'()*+ p&c ؕؔHmn_HnQCk'()*+ p&cpRؔHNpc}_'()*+ p&ai\X{[ accVWn'()*+ p&cؓؕHmn_Hn|oDm'()*+ p&cd4ؕ@ F qJ'()*+ p&c؟ @'()*+ p&c؟ؐHmn_HnE'()*+ p&aaD@ u\?'()*+ p&a, )@ uX|'()*+ p&a_@ uZ'()*+ p&a{BYF@ uZ'()*+ p&aXgHh@ j6'()*+ p&a ']@ &ZB#'()*+ p&c؟؞ ؚ<R'()*+ p&aZd)p8@ {Yn'()*+ p&a~4!MK@ }[m'()*+ p&c=!@00[ '()*+ p&a<5BV@ #Zj '()*+ p&a?sL@:Q'()*+ p&a:c@ Pi^'()*+ p&ai'ZD`@ &*'()*+ p&adu@ {p'()*+ p&aw!tCV@ }n'()*+ p&af5@Xt@ #i'()*+ p&a2o+Z@ P'()*+ p&c@ؕ@;;a '()*+ p&ch(ؕ@Ju '()*+ p&a~y@t'()*+ p&a@|'()*+ p&c@86!'()*+ p&a C@1 '()*+ p&c DxؗؕM }  QL'()*+ p&ad x ;!Uww_'()*+ p&c`  x@PG`S ~7k'()*+ p&s z @ة_Je-qN'()*+ p&cxؔ$m> q'()*+ p&ax (xB55 '()*+ p&c!lؕ w.'()*+ p&cq uk'()*+ p&a0 ss88<<'()*+ p&c2wؔ؝_ y x ?'()*+ p&c ؟z؟C^]   oKLR'()*+ p&&&cؕ@86'()*+ p&a&mعw@ o'()*+ p&c@ '()*+ p&aA؞8@  a'()*+ p&a7@'()*+ p&c")V "'()*+ p&c[ؔ؞@sc[NJh'()*+ p&ahMs @UVuw'()*+ p&a#%J@  k#'()*+ p&f؟oN5h r 1'()'4  *+ p&c{;U 5&Gub'()*+ p&cZ @86Mx$'()*+ p&d% z'()*+ p&aWL@G\'()*+ p&c؟؟@86G '()*+ p&c+-a{\ wIN'()*+ p&a@U '()*+ p&ci@P Q '()*+ p&c@PÙ <'()*+ p&c}ؕ3j  ~nci'()*+ p&au R@c p'()*+ p&a  RY`@ '()*+ p&a1-L4 0++,9Q'()*+ p&cI c @ d'()*+ p&a,a-5 8C=P{{\1'()*+ p&&cy8<l X}'()*+ p&&a-\ P8CXPأأII|'()*+ p&c2[ؕ@A?o'()*+ p&ac[@ 1aC'()*+ p&c@`|'()*+ p&a@z'()*+ p&c  @86&'()*+ p&a-](خ@  )'()*+ p&cؚؔ  v B'()*+ p&s6ؓ@ug}"'()*+ p&apUII ww؂؂uJ5bu'()*+ p&a$s49@'()*+ p&cZ>*  Mx#Q#\UX'()*+ p&cؕ@86!('()*+ p&&aU+A>\ FF#.'()*+ p&ci؟ $؛\y:T'()*+ p&cCؔ؞@P%P*'()*+ p&&&a!@({'()*+ p&&aL%n@!*Gs'()*+ p&&f`d2ؗؓ2 ,~1 yi')  ()' l} l} (*+ p&fre2ؘؔ2 -1  wCDE')()' w >*+ p&&&c@hijk01'()*+ p&a0@ 1'()*+ p&c]tؕ 2?gjF`'()*+ p&&&cjK3rM 5] dVh'()*+ p&c ؚؔ@~&69'()*+ p&dx 7'()*+ p&d ~x 8  '()*+ p&a6H@ N9M'()*+ p&cJ؊ؕ S:r W K'()*+ p&&c!ؚؔ@ w<B'()*+ p&&dx > '()*+ p&&cz#@mc@iA'()*+ p&ao@S@ ANt'()*+ p&a<B@ BG'()*+ p&&aUP4@ D2'()*+ p&&c"5&ؔ؜@gkwF/BG'()*+ p&a*dF 1ز>.@G5f'()*+ p&ct'ؔ@gkwHUI'()*+ p&a H ؼG@mIz'()*+ p&&aJ:,@BK)D'()*+ p&a}UM@ L 6'()*+ p&&aR0M""N::$'()*+ p&a&U7؇4@ pO6U'()*+ p&c؟'؝@ PP{i'()*+ p&c VQ*'()*+ p&cp=VRc'S'()*+ p&&cKv$VT>\i'()*+ p&d0 Ox U' F9 X'()*+ p&cj[ V[X'()*+ p&a=F M{{VVW8K'()*+ p&adVN e;;ررX##Ii'()*+ p&&a{2jA} ددZZZv<'()*+ p&c8ؚؕ@./[od\^_a!g'()*+ p&a[،ػo@ ,\@%'()*+ p&c9lؕ=WW ]'F'()*+ p&at[ زH@ ^4'()*+ p&a[@ _'()*+ p&&aR[uظY@ a'()*+ p&d x b  '()*+ p&aSx  4NvvcN'()*+ p&a>axc:@F/d#l'()*+ p&awg[xJ  * *ززDDeXXEb|'()*+ p&a: j,Aaåززfoo w'()*+ p&cU ؕ@Pg?l'()*+ p&a@xZ@Xwh;_'()*+ p&aP,V/!@^ib'()*+ p&a@Bj '()*+ p&ax؄(@ؐJPk_}'()*+ p&&&&a&V `ZZVVo )'()*+ p&&&c P0lJ r<`c'()*+ p&amg4؈@Bs'()*+ p&fDؖ@t6u{}'()'&&*+ p&dt@u'()*+ p&cY@@v}}'()*+ p&c;3V5xO w.!YE'()*+ p&c: W) x-X'()*+ p&a.U2@ @y)7'()*+ p&c xؕ@*Rzri'()*+ p&a.tc @O {'()*+ p&cg؟؜ |ؚZtw'()*+ p&avtص&@}'()*+ p&cdS@~A'()*+ p&a0{ FFv'()*+ p&&frؘؔ3  Uhe('M(()'oXX1*+ p&cP ؔ@d]jkt'()*+ p&a%%0L2 1KtRR '()*+ p&a hx~@ '()*+ p&ai0U[ ==ooU'()*+ p&c)^ؔ@C 6k'()*+ p&a$i@@d)'()*+ p&aZ_ OOv'()*+ p&aVu ؃؃Qz'()*+ p&akm pr'()*+ p&c؟ؿ@3R/u '()*+ p&aK?*@1`p'()*+ p&ajB> ,n'()*+ p&c؟5@{{ '()*+ p&c |@vu *'()*+ p&aQ?@8`'()*+ p&aC@0 '()*+ p&a؟u@/'()*+ p&a0@BI5#'()*+ p&az:C-@BX$J'()*+ p&chؔ؞ w 89'()*+ p&af<#@/'()*+ p&&a@/'()*+ p&ao~@rj'()*+ p&ciYiؔ \zf M_'()*+ p&aQQعZZQQعع'()*+ p&fx2؟؎2 F1 '9x ')N()' Z Z*+ p&c8&y" +Vf_'()*+ p&&s8|؆ؓ I '&N'()*+ p&&aG@0 B'()*+ p&&c ؕ iin'()*+ p&&aS|ج6W@0 N('()*+ p&&cq~ dcFkm9'()*+ p&aGg&&Bl'()*+ p&c؟ؚ ؚQ'()*+ p&&&cOZk{ B8m|'()*+ p&adR ] ZZؼؼ((.8iQ'()*+ p&sCؓ@y<q'()*+ p&a>0q@0 r'()*+ p&f=8G@7yc "'k{()'q*+ p&&g, ) ,2 _ : g'( $Detritus)*+ p&as6ة@1g;'()*+ p&&a4@19'()*+ p&&cؕ]a '"V'()*+ p&cؔ؜}!% '()*+ p&&&cX@pHp!'()*+ p&&&&aH-J[@4HC`'()*+ p&&&a7T'?@ZZ..==oC'()*+ p&&c<$j؉ؕIW@I'()*+ p&&cjؾؕIW@'()*+ p&&&&a) df@f"4a_'()*+ p&fN5h ?@ |&'G()' l l*+ p&a vx*2(@k^p  :'()*+ p&&fya@ime/(m^.pP~[Oo'XX()'X[[*+ p&d% '()*+ p&cajؙ؍IW@Gn'()*+ p&c|j؏ؕV'()*+ p&a/jؗ+>m`XhDD A'()*+ p&a:V@B}Q'()*+ p&&aj]-0X66   2'()*+ p&axj؏.GpX&&//,=$'()*+ p&a`j[XMMVe'()*+ p&aW+lY\77UU%YL'()*+ p&dS,jn} wJ\'()*+ p&aDZZI'()*+ p&aBZZ G'()*+ p&aBZZkkG'()*+ p&a2j@ ,'()*+ p&as=j(b@ ,`'()*+ p&&&&a#@B ( '()*+ p&caؕؕd6  n'()*+ p&aL@B '()*+ p&aL@B '()*+ p&a"@B''()*+ p&c؟;UX u'()*+ p&c&ؔ؛~p?c3,4'()*+ p&c8ixؔ؞@TY+Ev'()*+ p&af4/QD,@B m'()*+ p&&&c'ؙؔ  s 84'()*+ p&cؔ؝ 9 t 2'()*+ p&&a>x[S q  JJ}}77~9X'()*+ p&a9xb  < ..YY4g'()*+ p&a.x7؆ Vp-ee{{E'()*+ p&a1xv Eؖؖ||,{'()*+ p&ax: H b tt ?'()*+ p&&ajاا11IIe'()*+ p&&a/x< t }   WW*A'()*+ p&ax-Jm@4U h}'()*+ p&dj '()*+ p&a1cx)5V  ~ 22;'()*+ p&sb -@ߦ4`Rn!'()*+ p&&cFAaؕ@JAJA(S- )lp'()*+ p&c'L ؟   ?8 6!'()*+ p&aAK3@2{".F'()*+ p&&&&dU '()*+ p&dOI FX '()*+ p&d '()*+ p&d9Sv '()*+ p&f^@!S 'K()'n*+ p&dSj@ J\'()*+ p&aO?j,-aåثث \\.v'()*+ p&a"+XOQ^kll zT'()*+ p&c`[  Gr '()*+ p&ae r e;;:: ##`w'()*+ p&ad j,؄ئaåخخ00Pt'()*+ p&aqj,aådv '()*+ p&d+ "4'()*+ p&&d '()*+ p&dD ;M'()*+ p&aqj@l'()*+ p&&aX2j{ظI@'()*+ p&&&&aS 880A'()*+ p&a2H-0X  21 GG-+]M'()*+ p&&fv؟ؒIc ?3  'G()'# A# v v*+ p&&&&aSu-T0@P!+OXz'()*+ p&&&&&aG1Kt&B'()*+ p&aAMk ^^'GG.'()*+ p&&cؕ@A40M)+H'()*+ p&aCm@ *>r'()*+ p&aP)nؑ l@2{"+'()*+ p&fx ~ ,gS~ f 'xKx()'nn,,*+ p&g?@I-Q'($Light)*+ p&sw9.bC{'()*+ p&fkH /vZvPOB'kk~~()'@@*+ p&f1Iؖ0vk1CI'c()'}}*+ p&d01'()*+ p&&aUxc@3j`3Ph'()*+ p&aZ-؂ا; #p=4iW'()*+ p&c؟ؕYs6 5'()*+ p&abJ-P@}6K!g'()*+ p&&a+ rr80'()*+ p&a؀ ==669'()*+ p&a~DxH@3j:e'()*+ p&&a@4 P<|'()*+ p&a9#P@4w=<r'()*+ p&fGn >Kc .?@E'SGSG()' K  *+ p&dI>x@ ?@R'()*+ p&a x> ;U@OF'()*+ p&c_n AR} E'()*+ p&a<6/44@066B55cu'()*+ p&asT.0 J0--C336h'()*+ p&cN// N#D/I'()*+ p&a[Ax>J bHHQQEE`'()*+ p&c  F)G'()*+ p&aFxK (B9ؼؼ00G00'()*+ p&c ؟?Yd H PQ'()*+ p&aD0+Jx`DDI0'()*+ p&&ax  4)33K  '()*+ p&a <x $ 6_ؔؔL  A'()*+ p&aaa6 ؤؤVVMll1\f'()*+ p&a2'@LN'()*+ p&dx@O'()*+ p&c 9.xC! P  TQR'()*+ p&a .PxH  G,$$Q\\  3'()*+ p&a /Px  ~|؅؅R<<  4'()*+ p&a/T-)K@tS]~'()*+ p&s ؓ@6MˡTUV['()*+ p&aT d@MU_'()*+ p&a)T5؈@N Vj'()*+ p&s؟ؓ@6O+WXZ\'()*+ p&aW e@ODX`'()*+ p&&a=W؈&@ORZ'()*+ p&aTg )i@P['()*+ p&aWg@P\'()*+ p&cW@SE0]09`'()*+ p&&aDL& ؾؾii_BB?!oR'()*+ p&a]rkX@`'()*+ p&cp@@\SȈaEbfd'()*+ p&caؕ@O`bHc'()*+ p&a br@c'()*+ p&c-W@&Pd Ne'()*+ p&adZ@&Pe2'()*+ p&ca@&Q`f"g'()*+ p&afت}@&Q g'()*+ p&c 5xؕ@T*R h D'()*+ p&azx @i'()*+ p&aqU-R@4DjlMZ'()*+ p&a~-[@4RkyV'()*+ p&&&cx@xno'()*+ p&ainx!}@9@o|'()*+ p&cj@~bpus~'()*+ p&a7ج3k@28qv'()*+ p&aLXQNDhyyrV'()*+ p&apjc@ s]'()*+ p&d3x t*o<'()*+ p&aR؆@Lu'()*+ p&c, @vo%|}'()*+ p&d- w'()*+ p&dF x=O'()*+ p&d, y'()*+ p&&&a` xv@8|['()*+ p&a xvؕ&@8}'()*+ p&a{pj A@ ~<%'()*+ p&c@ Vp'()*+ p&amy qmmVVh'()*+ p&ac3 UoL((x'()*+ p&c@ Vpt&'()*+ p&cm HW@؃_` ?'()*+ p&ckHa@؃_Ң^_'()*+ p&&ajdj@&4eo'()*+ p&ak di@&4dp '()*+ p&&a^gL @BY'()*+ p&caS@10T0v'()*+ p&a[]}O@ JX`'()*+ p&c؟ؕWvq upp'()*+ p&&&&&cF=. #i'()*+ p&cIG. <%gj'()*+ p&cI+ &l'()*+ p&cxN؜   k%j'()*+ p&&&&&c@eؕ@)\'()*+ p&sT-ؗؓ@J~;j'()*+ p&c4ؕ@B4̥K'()*+ p&g)%)C 2 K'($Global)*+ p&aBpj]صA{@ {2'()*+ p&c:  ؕ*D+ -p'()*+ p&a5@0  1'()*+ p&cvؕ ؛V<='()*+ p&cG; ؜:'eO'()*+ p&a/7r@4G;ER'()*+ p&cyX-@؈ lDl'()*+ p&a}\AM|y@itW'()*+ p&c$@ػ 6 '()*+ p&aD>x @zPI'()*+ p&ax@z'()*+ p&c@e.z`! !k'()*+ p&&cq2@_֒d|'()*+ p&cg_B@ZLr'()*+ p&auء1?@6j2Z'()*+ p&ankh@6cfs'()*+ p&cؕ@)_`0a'()*+ p&a(d@F6[g}'()*+ p&cLgw   Zn'()*+ p&cN 6 | Al J'()*+ p&a-W!@s1 &'()*+ p&s0ؓ(xB v'()*+ p&ay-"l@}g~''()*+ p&aWpj`@ |@'()*+ p&cG @Br'()*+ p&av |0J@m 'E'()*+ p&cq؟؞( dd'()*+ p&cDؚؔ@UQ'()*+ p&aWxt@Hi'()*+ p&a:Wxf@i/5k'()*+ p&ct؟x؟@߿g'()*+ p&cw؟ j '()*+ p&cM؋ؕG '%'()*+ p&axx@:xzs'()*+ p&cq؟؞@86d'()*+ p&cr؟ؿ@{c{e'()*+ p&ajxx%@|x'()*+ p&avx@T|q'()*+ p&ct؟))x@(40g'()*+ p&avx@:ɺ}q'()*+ p&c 33@0.] * a'()*+ p&a (x@ -'()*+ p&a;uxB@z6pG'()*+ p&aKj@zFo'()*+ p&c؟W@ 2?-'()*+ p&a@?v'()*+ p&cؕ@#@Vp'()*+ p&&a[-@41 V'()*+ p&c@]3 &'()*+ p&a @3I{'()*+ p&c@!eVp '()*+ p&a Rv؈o@m H_'()*+ p&c@ Vp'()*+ p&c)؟d؞6 '()*+ p&c3@؜4P'()*+ p&ai2&@ o8['()*+ p&c-a{\ }b'()*+ p&a-@47'()*+ p&c؇ؕ@EnAf l{i'()*+ p&axؔ)@|@c'()*+ p&aN@ m1ee#'()*+ p&a| p@m nK{'()*+ p&cDؕ@\&Q'()*+ p&&svS-@G3 f9i'()*+ p&aT-jE@7P.Kl^'()*+ p&a> W@m PN'()*+ p&&&&&&&a-g *@BZ('()*+ p&aT6L@B Y'()*+ p&aW6/@B* \'()*+ p&c b7ؕ@` @ H-wxi'()*+ p&cL_@` @ ?Gjwyn'()*+ p&a6KjF@0 1eK'()*+ p&a {7**'()*+ p&a`SkO@0 Jfe'()*+ p&a%nW@0  \'()*+ p&&a'X(H@0 "#M]'()*+ p&aobq@V?l]t'()*+ p&&&&&&&&&&&&&&&a?!}@3 xD'()*+ p&&&&aV-@~3 '()*+ p&c`\؞ؕ@55> '()*+ p&aX@ '()*+ p&a%@   '()*+ p&aR ea@f"4 M`f'()*+ p&&&cw@ػ j'()*+ p&c{@ػn8'()*+ p&&&&&&&&&dj@WP'()*+ p&fػ؋ؗ@ػp_D'()'j*+ p&aoj@ '()*+ p&apj@ k'()*+ p&fb\ؖ@p?1m 'b9bG()' *+ p&&&&&&&&&&&&&d  ()'()*+ p&&&c~tC + g4  '()*+ p&&&&c؟u؝ / 5'()*+ p&c`vؔ؝N" 0m6'()*+ p&&&&&a/=jj5YY'()*+ p&aX01Skk6%%f'()*+ p&&&&&&&&&&&&&&&&&&&&&&&&&cb؟nؚrf OUo]'()*+ p&&cZ  SPm QMx-0,'()*+ p&&&&&&&&&&&&&&&&&a| 0M''ررb!!w'()*+ p&&ck0lJ d^'()*+ p&c &5xO e " 9'()*+ p&c< %) f/ Z 8'()*+ p&c` g$hP$[\'()*+ p&&&&&c*%ؕLf ls'()*+ p&&&c}$ op@'()*+ p&&&&&&&&&&&&&&&&&&c:) {'()*+ p&&&&&&&&&&&&&&&&&gZ0 Z A= I'( $seaweed)*+ p&&&&&&&&&&&&&&&&alg JJ#'()*+ p&a.g  bةةDD3'()*+ p&&&cbId U3_'()*+ p&&&&&&&&&&&&&dI @R%'()*+ p&fMؘؔ 3X 'M$M()'*+ pdrmdΏ2 *66ňHbxp** 0ZZ ***@ᢋ@*2Untitled Graph1D2D #+?+ I^@4$I$:?/lW;?* 3+&" !0*nI7c#Έp** 0**pl@ᢋ@*:Untitled Graph1 2 #6@fffff;@Rϖvv 36&2 #9@[;?d`39&" !0*hp xhppP*d 0*pHpD@ᢋ@*DUntitled Graph1 2 #<?+jf@3<&2 #1?E˰@_@31&0*p(p$@2@*DUntitled Graph1T2 #A@>rZ1>f3&2 #@K?^!t_* @$trTz3&2 # @vUUUUV?I/t Pm@.C!3&2 #B@888@R3B&2 # @B @J7Pq3&0*pp@2@*DUntitled Graph1T2 #uA@>Z1>f3u&2 #w?QR@]s3w&2 #? @fffff?<ˣ@T%vI3?&2 #@B ?]-V@QQ#|3&2 #z @UUUUV?҅!sw@堋193z&0*pp @2@*DUntitled Graph12 #?@sTI$I?]?"  3&2 #?@sTI$I?(\)? 3&2 # ?@k?؜y?؜y3&2 #/?@{?gFPM?V{v?3/&" !0*,q00jp`ppj*( 0pToo@I@*sUntitled Graph12 #@"@D5UUUUU>@"$'@-I53&0pToo@ᢋ@*stide12 #%@I?++>@NPc@reg3+3%&2 #@I?++>@LԆɬ@#DB3&2 # @I?++>@K K@3&2 #@I?>@I ]9T@=K~R3&0pTo@I@*s Untitled1"j !0* qo* 0po@x@*Untitled Graph1" !0*"rHpo* 0ooo@ᢍ@*Untitled Graph12 #|?E˰@@3|&0ooo@ᢍ@* Untitled1 2 #c@R 3c&2 #i@k?3i&0oox@јX@* Untitled1" !0*s@`L8pop~* 0ott@ᢋ@*Untitled Graph12 #@m~?`bM񖩟3&"~ !0*;u :0 pt*7 0ottt@ᢋ@*Untitled Graph12 #0@@?ǦE˰@l.@30&0ottt@ᢋ@* Untitled12 #@K3&2 #?@sTI$I?P)b2+?Q/@&3&2 #o ?@sTI$I?xws_?xtj~3o&0ottt@ᢋ@* Untitled12 #@k?3&2 #@k?3&2 # @k?3&2 #@k?3&" !0*,u`.ptX*( 0ttPtL@ᢋ@*Untitled Graph12 #@fffff@-@?Gp3&0tt8t4@ᢋ@* Untitled12 #@j?htj~3&0tt t@јX@* Untitled12 #@$?v.3&0ttt@јX@* Untitled12 #M@k?~ t3M&" !0*v`@H24pt* 0t\ts@ᢋ@*Untitled Graph12 #3?~K 4a:Kx@n33&" !0*x ڥP6<px* 0t xw@ᢋ@*Untitled Graph12 #@k?zG{?zG{3&" !0*VxL8٬pw* 0x w@I@*Untitled Graph10x ww@ᢋ@* Untitled12 #-?E˰@$?){?3-&" !0./ 7~6 5 /ww`w@w`  ww|     drmd N 4 ww\  4   drmd ~  ww<      drmd  www 56T, #%)/027bbw$}tF L(89:-  , Geneva.+ $Setup #1($t11/10/99 10:53 AM(?Input Variables(URun #)iswitch air temp lim(k1)_0.00({2)_1.000 **r ;v2<=>Geneva>Geneva;? @ /k6@'f4A Q k*A:&*2AMHj*sA*$D*A$>*A'1A*A#*A2*A*A ;%'*DAUA*:AL8*