Using the following information which are model equations for calculating the total N2O flux from denitrification. Write an R code to calculate the total N2O fluxes from denitrification. Define all input variables and parameters. The input variables include soil temperature, soil moisture, soil pH, nitrate concentration, WFPS, and carbon content. The rest of the variables will be parameters. Also, please wrap the code in a function such that it can loop through the input variables when an excel sheet of input data is read into the code. 
 "N2O from denitrification was obtained by a two-part addition: (1) developing a set of equations to model the total denitrification rate (N2+ N2O) and (2) partitioning N2 from N2O. 
D\_Ntotal= min[Fd(NO3), Fd(C)]\*Fd(WFPS)\*Fd(T)\*Fd(pH)

where D\_Ntotal is the denitrification rate per unit area (production of N2 + N2O, g N ha-1 d-1); Fd(NO3) is the maximum total N gas flux per unit area for a given soil NO3 level in g N ha-1 d-1 (assuming total soil C is not limiting); Fd(C) is the maximum total N gas flux
per unit area for a given total soil C level in g N ha-1 d-1 (assuming soil N is not limiting); and Fd(WFPS), Fd(T), and Fd(pH) are functions that represent the effects of soil moisture, soil temperature, and pH on N gas flux, respectively. 

Fd(NO3)=11,000+ (40,000+atan⁡(π\*0.002\*(NO3-180)))/π

Fd(C)= 24,000/(1+ 200/e^(0.35\*C) )-100

Fd(WFPS)= a/b^(b^((c/(d\*θ))) ) 

where Fd(NO3) and Fd(C) have units of g N ha-1 d-1, WFPS is water-filled pore space in units of m3/m3, C is the C content in g N ha-1 d-1, and a, b, c and d are soil-specific fitting parameters to incorporate the effect of soil texture.

Fd(T)=Max[0.9\*((Soil\_Temp)/(Soil\_Temp+exp⁡(9.93-0.312\*Soil\_Temp ) )+0.1),0.1]

Fd(pH)={0.001 for pH≤3.5,(pH-3.5)/3 for 3.5 <pH <6.5,1 for pH≥6.5

We developed a ratio method to differentiate the N2O from N2 produced during denitrification:

R\_(N2/N2O)=min⁡[Fr(NO3),Fr(C)]\*Fr(WFPS)\*Fr(pH)

Where R\_(N2/N2O) is the ratio N2 to N2O, and Fr(NO3), Fr(C), Fr(WFPS),and Fr(pH) represent the effects of soil NO3-N, total soil C, soil moisture, and pH on the ratio of N2 to N2O, respectively. Note here that temperature is not included in the ratio prediction as increasing
temperature increases both N2 and N2O production similarly. 
Fr(NO3)=1-[0.5+(1\*atan⁡(π\*0.01\*(NO3-190))/π]\*25

Fr(C)=13+(30.78\*atan⁡(π\*0.07\*(C-13)))/π

Fr(WFPS)= 1.4/13^(13^(17/((2.2\*WFPS))) ) 

Fr(pH)= 1/(1470\*e^(-1.1\*pH) )

Total N2O production during denitrification is given by:

D\_N2O= D\_Ntotal/(1+ R\_(N2/N2O) )
Where D\_N2O is total N2O flux per unit area in g N/ha d.
D\_(N2 )= D\_Ntotal/(1+ 〖1/R〗\_(N2/N2O) )
Where D\_(N2 ) is total N2 flux per unit area in g N/ha d"
탈질로 인한 총 N2O 플럭스를 계산하기 위한 모델 방정식인 다음 정보를 노래합니다. 탈질로 인한 총 N2O 플럭스를 계산하는 R 코드를 작성합니다. 모든 입력 변수와 매개 변수를 정의합니다. 입력 변수에는 토양 온도, 토양 수분, 토양 pH, 질산염 농도, WFPS 및 탄소 함량이 포함됩니다. 나머지 변수는 파라미터가 됩니다. 또한 엑셀 시트의 입력 데이터를 코드로 읽어들일 때 입력 변수를 반복할 수 있도록 코드를 함수로 감싸주세요.  "탈질로 인한 N2O는 (1) 총 탈질률(N2+ N2O)을 모델링하기 위한 방정식 세트를 개발하고 (2) N2O에서 N2를 분할하는 두 부분으로 나누어 구했습니다. D\_Ntotal= min[Fd(NO3), Fd(C)]\*Fd(WFPS)\*Fd(T)\*Fd(pH)여기서 D\_Ntotal은 단위 면적당 탈질 속도(N2 + N2O 생산량, g N ha-1 d-1)이고, Fd(NO3)는 주어진 토양 NO3 수준에 대한 단위 면적당 최대 총 N 가스 플럭스(g N ha-1 d-1)이며(총 토양 C가 제한되지 않는다고 가정), Fd(C)는 주어진 총 토양에 대한 단위 면적당 최대 총 N 기체 플럭스(단위 면적당 주어진 총 토양 C 수준에 대한 단위 면적당 최대 총 N 가스 플럭스(g N ha-1 d-1)(토양 N이 제한되지 않는다고 가정); Fd(WFPS), Fd(T) 및 Fd(pH)는 각각 토양 수분, 토양 온도 및 pH가 N 가스 플럭스에 미치는 영향을 나타내는 함수입니다. Fd(NO3)=11,000+ (40,000+atan⁡(π\*0.002\*(NO3-180)))/πFd(C)= 24,000/(1+ 200/e^(0.35\*C) )-100Fd(WFPS)= a/b^(b^((c/(d\*θ))) ) 여기서 Fd(NO3) 및 Fd(C) 단위는 g N ha-1 d-1, WFPS는 물로 채워진 공극 공간(m3/m3 단위), C는 g N ha-1 d-1 단위의 C 함량, a, b, c 및 d는 토양 질감의 효과를 통합하기 위한 토양별 피팅 파라미터입니다.Fd(T)=Max[0.9\*((Soil\_Temp)/(Soil\_Temp+exp⁡(9.93-0.312\*Soil\_Temp ) )+0.1),0.1]Fd(pH)={0.001 for pH≤3.5,(pH-3.5)/3 for 3.5 <pH <6.5,1 for pH≥6.5탈질 과정에서 생성되는 N2O와 N2를 구별하는 비율 방법을 개발했습니다:R\_(N2/N2O)=min[Fr(NO3),Fr(C)]\*Fr(WFPS)\*Fr(pH)여기서 R\_(N2/N2O)는 N2 대 N2O의 비율이고, Fr(NO3), Fr(C), Fr(WFPS) 및 Fr(pH)는 각각 토양 NO3-N, 총 토양 C, 토양 수분 및 pH가 N2 대 N2O의 비율에 미치는 영향을 나타냅니다. 여기서 온도는 비율 예측에 포함되지 않습니다.온도가 증가하면 N2와 N2O 생산량이 비슷하게 증가하기 때문입니다. Fr(NO3)=1-[0.5+(1\*atan⁡(π\*0.01\*(NO3-190))/π]\*25Fr(C)=13+(30.78\*atan⁡(π\*0.07\*(C-13)))/πFr(WFPS)= 1.4/13^(13^(17/((2.2\*WFPS))) ) Fr(pH)= 1/(1470\*e^(-1.1\*pH) )탈질 중 총 N2O 생성량은 다음과 같습니다:D\_N2O= D\_Ntotal/(1+ R\_(N2/N2O) )여기서 D\_N2O는 단위 면적당 총 N2O 플럭스(g N/ha d)입니다.D\_(N2 )= D\_Ntotal/(1+ 〖1/R〗\_(N2/N2O) )여기서 D\_(N2 )는 단위 면적당 총 N2 플럭스(g N/ha d")입니다.