Changeset 7239 for branches/ORCHIDEE_2_2/ORCHIDEE
- Timestamp:
- 2021-06-28T14:38:08+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/hydrol.f90
r7200 r7239 95 95 REAL(r_std), PARAMETER :: allowed_err = 2.0E-8_r_std 96 96 REAL(r_std), PARAMETER :: EPS1 = EPSILON(un) !! A small number 97 97 98 ! one dimension array allocated, computed, saved and got in hydrol module 98 99 ! Values per soil type 99 100 100 REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:) :: pcent !! Fraction of saturated volumetric soil moisture above 101 !! which transpir is max (0-1, unitless) 102 !$OMP THREADPRIVATE(pcent) 101 !! which transpir is max (0-1, unitless) 102 !$OMP THREADPRIVATE(pcent) 103 103 REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:) :: mc_awet !! Vol. wat. cont. above which albedo is cst 104 104 !! @tex $(m^{3} m^{-3})$ @endtex … … 440 440 kjpindex, index, rest_id, & 441 441 njsc, soiltile, veget, veget_max, & 442 humrel, vegstress, drysoil_frac,&442 humrel, vegstress, drysoil_frac, & 443 443 shumdiag_perma, qsintveg, & 444 444 evap_bare_lim, evap_bare_lim_ns, snow, snow_age, snow_nobio, & … … 450 450 !! 0.1 Input variables 451 451 452 453 !salma: added soil params as input variables 452 INTEGER(i_std), INTENT(in) :: kjit !! Time step number 453 INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size 454 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: index !! Indeces of the points on the map 455 INTEGER(i_std),INTENT (in) :: rest_id !! Restart file identifier 456 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class in the 457 !! grid cell (1-nscm, unitless) 458 ! 2D soil parameters 454 459 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: ks !! Hydraulic conductivity at saturation (mm {-1}) 455 460 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: nvan !! Van Genuchten coeficients n (unitless) … … 459 464 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcfc !! Volumetric water content at field capacity (m^{3} m^{-3}) 460 465 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcw !! Volumetric water content at wilting point (m^{3} m^{-3}) 461 !end salma: added soil params as input variables 462 463 INTEGER(i_std), INTENT(in) :: kjit !! Time step number 464 INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size 465 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: index !! Indeces of the points on the map 466 INTEGER(i_std),INTENT (in) :: rest_id !! Restart file identifier 467 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class in the grid cell (1-nscm, unitless) 466 468 467 REAL(r_std),DIMENSION (kjpindex,nstm), INTENT (in) :: soiltile !! Fraction of each soil tile within vegtot (0-1, unitless) 469 468 REAL(r_std),DIMENSION (kjpindex,nvm), INTENT (in) :: veget !! Fraction of vegetation type 470 469 REAL(r_std),DIMENSION (kjpindex,nvm), INTENT (in) :: veget_max !! Max. fraction of vegetation type (LAI -> infty) 471 470 471 472 472 !! 0.2 Output variables 473 473 REAL(r_std),DIMENSION (kjpindex,nvm), INTENT (out) :: humrel !! Relative humidity … … 492 492 REAL(r_std),DIMENSION (kjpindex) :: soilwetdummy !! Temporary variable never used 493 493 494 495 494 !! 0.4 Local variables 496 495 INTEGER(i_std) :: jsl 497 496 498 !salma: added soil params which depend on soil texture499 REAL(r_std), DIMENSION (nscm) :: nvan_text500 REAL(r_std), DIMENSION (nscm) :: avan_text501 REAL(r_std), DIMENSION (nscm) :: mcr_text502 REAL(r_std), DIMENSION (nscm) :: mcs_text503 REAL(r_std), DIMENSION (nscm) :: ks_text504 REAL(r_std), DIMENSION (nscm) :: pcent_text505 REAL(r_std), DIMENSION (nscm) :: mcfc_text506 REAL(r_std), DIMENSION (nscm) :: mcw_text507 REAL(r_std), DIMENSION (nscm) :: mc_awet_text508 REAL(r_std), DIMENSION (nscm) :: mc_adry_text509 !end salma: added soil params which depend on soil texture510 497 !_ ================================================================================================================================ 511 498 … … 622 609 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: reinf_slope !! Slope coef 623 610 624 !salma: added input soil params:625 611 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: ks !! Hydraulic conductivity at saturation (mm {-1}) 626 612 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: nvan !! Van Genuchten coeficients n (unitless) … … 630 616 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcfc !! Volumetric water content at field capacity (m^{3} m^{-3}) 631 617 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcw !! Volumetric water content at wilting point (m^{3} m^{-3}) 632 !end salma: added input soil params633 618 634 619 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: evapot !! Soil Potential Evaporation … … 695 680 !! at the present time-step @tex ($K$) @endtex 696 681 697 698 682 !! 0.4 Local variables 699 !salma: added variable kfact_root_type700 683 INTEGER(i_std) :: jst !! Index of soil tiles (unitless, 1-3) 701 684 INTEGER(i_std) :: jsl !! Index of soil layers (unitless) 702 685 INTEGER(i_std) :: ji, jv 703 CHARACTER(LEN=80) :: kfact_root_type !! read from run.def: when equal to 'cons', it indicates that ks does not increase 704 !! in the rootzone, ie, kfact_root=1; else, kfact_root defined as usual 686 CHARACTER(LEN=80) :: kfact_root_type !! read from run.def: when equal to 'cons', it indicates that 687 !! ks does not increase in the rootzone, ie, kfact_root=1; 688 !! else, kfact_root defined as usual 705 689 REAL(r_std),DIMENSION (kjpindex) :: soilwet !! A temporary diagnostic of soil wetness 706 690 REAL(r_std),DIMENSION (kjpindex) :: snowdepth_diag !! Depth of snow layer containing default values, only for diagnostics 707 REAL(r_std),DIMENSION (kjpindex, nsnow) :: snowdz_diag !! Depth of snow layer on all layers containing default values, only for diagnostics 691 REAL(r_std),DIMENSION (kjpindex, nsnow) :: snowdz_diag !! Depth of snow layer on all layers containing default values, 692 !! only for diagnostics 708 693 REAL(r_std),DIMENSION (kjpindex) :: njsc_tmp !! Temporary REAL value for njsc to write it 709 694 REAL(r_std), DIMENSION (kjpindex) :: snowmelt !! Snow melt [mm/dt_sechiba] … … 911 896 ! mcs etc are identical in all layers (no normalization by vegtot to be comparable to mc) 912 897 DO jsl=1,nslm 913 !salma: replaced mcs(njsc(:)) by mcs(:) and same for other variables914 898 land_mcs(:,jsl) = mcs(:) 915 899 land_mcfc(:,jsl) = mcfc(:) … … 1332 1316 !_ ================================================================================================================================ 1333 1317 !!_ hydrol_init 1334 !salma: added variables and njsc in arguments and input variables1335 1318 1336 1319 SUBROUTINE hydrol_init(ks, nvan, avan, mcr, mcs, mcfc, mcw, njsc,& … … 1344 1327 1345 1328 !! 0.1 Input variables 1346 !salma introduced njsc variable1347 1329 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class in the grid cell (1-nscm, unitless) 1348 1330 INTEGER(i_std), INTENT (in) :: kjit !! Time step number … … 1355 1337 !! 0.2 Output variables 1356 1338 1357 !salma: added input variables1358 1339 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: ks !! Hydraulic conductivity at saturation (mm {-1}) 1359 1340 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: nvan !! Van Genuchten coeficients n (unitless) … … 1389 1370 INTEGER(i_std) :: jsc !! Index of soil texture (1) 1390 1371 INTEGER(i_std), PARAMETER :: error_level = 3 !! Error level for consistency check 1391 !! Switch to 2 tu turn fatal errors into warnings1372 !! Switch to 2 tu turn fatal errors into warnings 1392 1373 REAL(r_std), ALLOCATABLE, DIMENSION (:) :: free_drain_max !! Temporary var for initialization of free_drain_coef 1393 1374 REAL(r_std), ALLOCATABLE, DIMENSION (:) :: zwt_default !! Temporary variable for initialization of zwt_force 1394 1375 LOGICAL :: zforce !! To test if we force the WT in any of the soiltiles 1376 1395 1377 1396 1378 !_ ================================================================================================================================ … … 1463 1445 !! 2.1 array allocation for soil texture 1464 1446 1465 1466 ALLOCATE (pcent(kjpindex),stat=ier) 1447 ALLOCATE (pcent(nscm),stat=ier) 1467 1448 IF (ier /= 0) CALL ipslerr_p(3,'hydrol_init','Problem in allocate of variable pcent','','') 1468 1449 1469 ALLOCATE (mc_awet( kjpindex),stat=ier)1450 ALLOCATE (mc_awet(nscm),stat=ier) 1470 1451 IF (ier /= 0) CALL ipslerr_p(3,'hydrol_init','Problem in allocate of variable mc_awet','','') 1471 1452 1472 ALLOCATE (mc_adry( kjpindex),stat=ier)1453 ALLOCATE (mc_adry(nscm),stat=ier) 1473 1454 IF (ier /= 0) CALL ipslerr_p(3,'hydrol_init','Problem in allocate of variable mc_adry','','') 1474 1455 1475 !! __2.2 Soil texture choose1456 !! 2.2 Soil texture choose 1476 1457 1477 1458 SELECTCASE (nscm) 1478 CASE (3) 1479 1480 pcent(:) = pcent_fao(njsc(:)) 1481 mc_awet(:) = mc_awet_fao(njsc(:)) 1482 mc_adry(:) = mc_adry_fao(njsc(:)) 1483 CASE (13) 1484 1485 pcent(:) = pcent_usda(njsc(:)) 1486 mc_awet(:) = mc_awet_usda(njsc(:)) 1487 mc_adry(:) = mc_adry_usda(njsc(:)) 1488 1459 CASE (3) 1460 pcent(:) = pcent_fao(:) 1461 mc_awet(:) = mc_awet_fao(:) 1462 mc_adry(:) = mc_adry_fao(:) 1463 CASE (13) 1464 pcent(:) = pcent_usda(:) 1465 mc_awet(:) = mc_awet_usda(:) 1466 mc_adry(:) = mc_adry_usda(:) 1489 1467 CASE DEFAULT 1490 1468 WRITE (numout,*) 'Unsupported soil type classification. Choose between zobler and usda according to the map' … … 1495 1473 1496 1474 !! 2.3 Read in the run.def the parameters values defined by the user 1497 1498 1499 1500 1501 1502 1503 1504 1505 1475 1506 1476 !Config Key = WETNESS_TRANSPIR_MAX … … 1520 1490 & "Please, check parameter value in run.def. ") 1521 1491 END IF 1522 1523 1524 1492 1525 1493 … … 2716 2684 2717 2685 !! 0.1 Input variables 2718 !salma: added the following soil parameters2719 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: ks !! Hydraulic conductivity at saturation (mm {-1})2720 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: nvan !! Van Genuchten coeficients n (unitless)2721 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: avan !! Van Genuchten coeficients a (mm-1})2722 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcr !! Residual volumetric water content (m^{3} m^{-3})2723 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcs !! Saturated volumetric water content (m^{3} m^{-3})2724 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcfc !! Volumetric water content at field capacity (m^{3} m^{-3})2725 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcw !! Volumetric water content at wilting point (m^{3} m^{-3})2726 2727 2686 ! input scalar 2728 2687 INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size (number of grid cells) (1) … … 2733 2692 !! in the grid cell (1-nscm, unitless) 2734 2693 REAL(r_std), DIMENSION (kjpindex,nstm), INTENT (in) :: soiltile !! Fraction of each soil tile within vegtot (0-1, unitless) 2735 2694 2695 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: ks !! Hydraulic conductivity at saturation (mm {-1}) 2696 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: nvan !! Van Genuchten coeficients n (unitless) 2697 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: avan !! Van Genuchten coeficients a (mm-1}) 2698 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcr !! Residual volumetric water content (m^{3} m^{-3}) 2699 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcs !! Saturated volumetric water content (m^{3} m^{-3}) 2700 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcfc !! Volumetric water content at field capacity (m^{3} m^{-3}) 2701 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcw !! Volumetric water content at wilting point (m^{3} m^{-3}) 2702 2736 2703 !! 0.2 Output variables 2737 2704 … … 2788 2755 !! (unitless) 2789 2756 !! need special treatment 2790 !salma: added local variable ii and jiref replaced with iiref2791 2757 INTEGER(i_std) :: ii 2792 2758 INTEGER(i_std) :: iiref !! To identify the mc_lins where k_lin and d_lin … … 2962 2928 !! 3 Compute the profile for a and n 2963 2929 !- 2964 !salma: for every grid cell2965 2930 DO ji = 1, kjpindex 2966 2931 DO jsl=1,nslm … … 2974 2939 ENDDO 2975 2940 ENDDO 2976 !end salma2977 2941 2978 2942 ! For every grid cell … … 3043 3007 ( frac**(-un/m) -un ) ** (-m) 3044 3008 ENDIF 3045 ENDDO !Salma end loop over landpoints3009 ENDDO 3046 3010 3047 3011 ! Special case for ii=imin … … 3138 3102 tmc_litter_field(ji,jst) = dz(2) * mcfc(ji) / deux 3139 3103 tmc_litter_sat(ji,jst) = dz(2) * mcs(ji) / deux 3140 tmc_litter_awet(ji,jst) = dz(2) * mc_awet( ji) / deux3141 tmc_litter_adry(ji,jst) = dz(2) * mc_adry( ji) / deux3104 tmc_litter_awet(ji,jst) = dz(2) * mc_awet(njsc(ji)) / deux 3105 tmc_litter_adry(ji,jst) = dz(2) * mc_adry(njsc(ji)) / deux 3142 3106 ENDDO 3143 3107 END DO … … 3163 3127 tmc_litter_awet(ji,jst) = tmc_litter_awet(ji,jst) + & 3164 3128 &(dz(jsl)+ dz(jsl+1))* & 3165 & mc_awet( ji)/deux3129 & mc_awet(njsc(ji))/deux 3166 3130 tmc_litter_adry(ji,jst) = tmc_litter_adry(ji,jst) + & 3167 3131 & (dz(jsl)+ dz(jsl+1))* & 3168 & mc_adry( ji)/deux3132 & mc_adry(njsc(ji))/deux 3169 3133 END DO 3170 3134 END DO … … 3660 3624 3661 3625 !! 0.1 Input variables 3662 !salma added soil params to arguments and input variables 3626 3627 INTEGER(i_std), INTENT(in) :: kjpindex 3628 REAL(r_std), DIMENSION (kjpindex,nvm), INTENT (in) :: veget_max !! Map of max vegetation types [-] 3629 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class 3630 !! in the grid cell (1-nscm, unitless) 3631 3663 3632 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: ks !! Hydraulic conductivity at saturation (mm {-1}) 3664 3633 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: nvan !! Van Genuchten coeficients n (unitless) … … 3668 3637 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcfc !! Volumetric water content at field capacity (m^{3} m^{-3}) 3669 3638 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcw !! Volumetric water content at wilting point (m^{3} m^{-3}) 3670 3671 3672 INTEGER(i_std), INTENT(in) :: kjpindex 3673 REAL(r_std), DIMENSION (kjpindex,nvm), INTENT (in) :: veget_max !! Map of max vegetation types [-] 3674 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class 3675 !! in the grid cell (1-nscm, unitless) 3639 3676 3640 REAL(r_std), DIMENSION (kjpindex,nstm), INTENT (in) :: soiltile !! Fraction of each soil tile within vegtot (0-1, unitless) 3677 3641 REAL(r_std), DIMENSION (kjpindex,nvm), INTENT(in) :: transpir !! Transpiration … … 4393 4357 ! sm_nostress = soil moisture of each layer at which us reaches 1, here at the middle of [smw,smf] 4394 4358 DO jsl = 1,nslm 4395 sm_nostress(:,jsl) = smw(:,jsl) + pcent( :) * (smf(:,jsl)-smw(:,jsl))4359 sm_nostress(:,jsl) = smw(:,jsl) + pcent(njsc(:)) * (smf(:,jsl)-smw(:,jsl)) 4396 4360 END DO 4397 4361 … … 4979 4943 4980 4944 ! GLOBAL (in or inout) 4981 ! salma added input soil hydraulic parameters4982 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: ks !! Hydraulic conductivity at saturation (mm {-1})4983 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: nvan !! Van Genuchten coeficients n (unitless)4984 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: avan !! Van Genuchten coeficients a (mm-1})4985 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcr !! Residual volumetric water content (m^{3} m^{-3})4986 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcs !! Saturated volumetric water content (m^{3} m^{-3})4987 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcfc !! Volumetric water content at field capacity (m^{3} m^{-3})4988 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcw !! Volumetric water content at wilting point (m^{3} m^{-3})4989 4990 4991 4945 INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size 4992 4946 INTEGER(i_std), INTENT(in) :: ins 4993 4947 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class in the grid cell 4994 4948 !! (1-nscm, unitless) 4949 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: ks !! Hydraulic conductivity at saturation (mm {-1}) 4950 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: nvan !! Van Genuchten coeficients n (unitless) 4951 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: avan !! Van Genuchten coeficients a (mm-1}) 4952 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcr !! Residual volumetric water content (m^{3} m^{-3}) 4953 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcs !! Saturated volumetric water content (m^{3} m^{-3}) 4954 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcfc !! Volumetric water content at field capacity (m^{3} m^{-3}) 4955 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcw !! Volumetric water content at wilting point (m^{3} m^{-3}) 4995 4956 REAL(r_std), DIMENSION (kjpindex), INTENT (in) :: flux_infilt !! Water to infiltrate 4996 4957 !! @tex $(kg m^{-2})$ @endtex 4997 4998 4958 !! 0.2 Output variables 4999 4959 REAL(r_std), DIMENSION(kjpindex,nstm), INTENT(out) :: check !! delta SM - flux (mm/dt_sechiba) … … 5163 5123 !! 0.1 Input variables 5164 5124 5165 !salma: added parameter mcr5166 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcr !! Residual volumetric water content (m^{3} m^{-3})5167 5125 INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size 5168 5126 INTEGER(i_std), INTENT(in) :: ins !! Soiltile index (1-nstm, unitless) 5169 5127 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class in grid cell 5170 !! (1-nscm, unitless) 5128 !! (1-nscm, unitless) 5129 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcr !! Residual volumetric water content (m^{3} m^{-3}) 5171 5130 5172 5131 !! 0.2 Output variables … … 5321 5280 5322 5281 !! 0.1 Input variables 5323 !salma: added mcs5324 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcs !! Saturated volumetric water content (m^{3} m^{-3})5325 5282 INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size 5326 5283 INTEGER(i_std), INTENT(in) :: ins !! Soiltile index (1-nstm, unitless) 5327 5284 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class in grid cell 5328 5285 !! (1-nscm, unitless) 5286 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcs !! Saturated volumetric water content (m^{3} m^{-3}) 5329 5287 5330 5288 !! 0.2 Output variables … … 5454 5412 5455 5413 !! 0.1 Input variables 5456 !salma: added mcs5457 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcs !! Saturated volumetric water content (m^{3} m^{-3})5458 5414 INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size 5459 5415 INTEGER(i_std), INTENT(in) :: ins !! Soiltile index (1-nstm, unitless) 5460 5416 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class in grid cell 5461 5417 !! (1-nscm, unitless) 5418 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcs !! Saturated volumetric water content (m^{3} m^{-3}) 5462 5419 5463 5420 !! 0.2 Output variables … … 5740 5697 5741 5698 !! 0.1 Input variables 5742 !salma: added mcr and mcs 5699 INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size 5700 INTEGER(i_std), INTENT(in) :: ins !! Index of soil type 5701 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class 5702 !! in the grid cell (1-nscm, unitless) 5743 5703 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcr !! Residual volumetric water content (m^{3} m^{-3}) 5744 5704 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcs !! Saturated volumetric water content (m^{3} m^{-3}) 5745 INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size5746 INTEGER(i_std), INTENT(in) :: ins !! Index of soil type5747 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class in the grid cell (1-nscm, unitless)5748 5705 5749 5706 !! 0.2 Output variables … … 5832 5789 5833 5790 !! 0.1 Input variables 5834 !salma: added the following soil variables 5791 5792 INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size 5793 INTEGER(i_std), INTENT(in) :: ins !! Index of soil type 5794 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class 5795 !! in the grid cell (1-nscm, unitless) 5835 5796 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: nvan !! Van Genuchten coeficients n (unitless) 5836 5797 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: avan !! Van Genuchten coeficients a (mm-1}) 5837 5798 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcr !! Residual volumetric water content (m^{3} m^{-3}) 5838 5799 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcs !! Saturated volumetric water content (m^{3} m^{-3}) 5839 5840 INTEGER(i_std), INTENT(in) :: kjpindex !! Domain size5841 INTEGER(i_std), INTENT(in) :: ins !! Index of soil type5842 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class in the grid cell (1-nscm, unitless)5843 5800 5844 5801 !! 0.2 Output variables … … 6349 6306 6350 6307 !! 0.1 Input variables 6351 !salma: added the following soil variables 6308 6309 ! input scalar 6310 INTEGER(i_std), INTENT(in) :: kjpindex 6311 REAL(r_std),DIMENSION (kjpindex,nvm), INTENT (in) :: veget_max !! Max. vegetation type 6312 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class in the grid cell (1-nscm, unitless) 6313 REAL(r_std), DIMENSION (kjpindex,nstm), INTENT (in) :: soiltile !! Fraction of each soil tile within vegtot (0-1, unitless) 6352 6314 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: ks !! Hydraulic conductivity at saturation (mm {-1}) 6353 6315 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: nvan !! Van Genuchten coeficients n (unitless) … … 6357 6319 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcfc !! Volumetric water content at field capacity (m^{3} m^{-3}) 6358 6320 REAL(r_std),DIMENSION (kjpindex), INTENT (in) :: mcw !! Volumetric water content at wilting point (m^{3} m^{-3}) 6359 6360 ! input scalar6361 INTEGER(i_std), INTENT(in) :: kjpindex6362 REAL(r_std),DIMENSION (kjpindex,nvm), INTENT (in) :: veget_max !! Max. vegetation type6363 INTEGER(i_std),DIMENSION (kjpindex), INTENT (in) :: njsc !! Index of the dominant soil textural class in the grid cell (1-nscm, unitless)6364 REAL(r_std), DIMENSION (kjpindex,nstm), INTENT (in) :: soiltile !! Fraction of each soil tile within vegtot (0-1, unitless)6365 6321 REAL(r_std), DIMENSION (kjpindex), INTENT(in) :: evapot !! 6366 6322 REAL(r_std), DIMENSION (kjpindex), INTENT(in) :: returnflow !! Water returning to the deep reservoir
Note: See TracChangeset
for help on using the changeset viewer.
