Changes between Version 3 and Version 4 of Branches/MergeOCN/Peng

Timestamp:

2015-05-06T13:28:38+02:00 (9 years ago)

Author:

speng

Comment:

--

Branches/MergeOCN/Peng

@@ -6 +6 @@
 
    line 4972, soiltext ---> soilbulk
+
+
    line 4975, soiltext ---> soilph
 
 2. corrected bug as ticket#109
 
+3. corrected bug for "carbon"initialization, actually, this variable should contain a element dimension in the calculation, at least icarbon and initrogen. So it is better to rename this variable as ele_pools...
+
+in stomate_io.f90
+
+{{{
+
+ 232 !!!!pssdebugok#1+++
+ 233     REAL(r_std), DIMENSION(npts,ncarb,nvm,nelements), INTENT(out)       :: carbon                  !! som pool: active, slow, or passive, surface
+ 234                                                                                                    !! for carbon and nitrogen (gC/m**2)
+ 235 !!!!pssdebugok#1---
+ 236 !    REAL(r_std), DIMENSION(npts,ncarb,nvm), INTENT(out)                 :: carbon                  !! carbon pool: active, slow, or passive,
+ 237 !                                                                                                   !! (gC/m**2)
+
+}}}
+ 
+4. leaching only have iammonium and initrate, other gas (inox, initrous) for leaching should be zero;
+
+in stomate_soilcarbon.f90
+
+debug: initiate leaching(:,:,:) = 0
+
+
+{{{
+ 654 !!!!pssdebug++++
+ 655     leaching(:,:,:) = 0
+ 656 !!!!pssdebug----
+
+}}}
+
+For denitrification calculation, there is no leaching(:,m,inox) and leaching(:,m,initrous) calculation before. I think we do not need to minus leaching for this two species here. 
+
+{{{
+1079 !!!!pssdebugok+          !
+1080           ! NO consumption
+1081           denitrification(:,m,i_nh4_to_no) = MIN(soil_n_min(:,m,inox), &
+1082                anvf(:,m) * ( mu_no(:) / Y_no + M_no * soil_n_min(:,m,inox) / sum_n(:)) * bact(:,m) * &
+1083                24. * dt )
+1084           !
+1085           ! N2O consumption
+1086           denitrification(:,m,i_nh4_to_n2o) = MIN(soil_n_min(:,m,initrous), &
+1087                anvf(:,m) * ( mu_n2o(:) / Y_n2o + M_n2o * soil_n_min(:,m,initrous) / sum_n(:)) * bact(:,m) * &
+1088                24. * dt )
+1089 !!!!pssdebugok-  
+1090 !!!!pssdebug          ! NO consumption
+1091 !!!!pssdebug          denitrification(:,m,i_nh4_to_no) = MIN(soil_n_min(:,m,inox)-leaching(:,m,inox), &
+1092 !!!!pssdebug               anvf(:,m) * ( mu_no(:) / Y_no + M_no * soil_n_min(:,m,inox) / sum_n(:)) * bact(:,m) * &
+1093 !!!!pssdebug               24. * dt )        
+1094 !!!!pssdebug          !
+1095 !!!!pssdebug          ! N2O consumption
+1096 !!!!pssdebug          denitrification(:,m,i_nh4_to_n2o) = MIN(soil_n_min(:,m,initrous)-leaching(:,m,initrous), &
+1097 !!!!pssdebug               anvf(:,m) * ( mu_n2o(:) / Y_n2o + M_n2o * soil_n_min(:,m,initrous) / sum_n(:)) * bact(:,m) * &
+1098 !!!!pssdebug               24. * dt )
+
+}}}
+
+5. nitrogen demand problem