Changeset 7339
 Timestamp:
 11/06/21 13:15:01 (6 months ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/slowproc.f90
r7338 r7339 2833 2833 INTEGER(i_std), DIMENSION(:), ALLOCATABLE :: vecpos 2834 2834 REAL(r_std) :: sgn !! sum of fractions excluding glaciers and ocean 2835 2836 ! For the calculation of field capacity and wilting point 2837 REAL(r_std),DIMENSION (nbpt) :: mvan !! Van Genuchten parameter m 2838 REAL(r_std),DIMENSION (nbpt) :: psi_w !! Matrix potential characterizing the wilting point (mm) 2839 REAL(r_std),DIMENSION (nbpt) :: psi_fc !! Matrix potential characterizing the field capacity (mm) 2840 2835 2841 !_ ================================================================================================================================ 2836 2842 … … 3466 3472 WRITE(numout,*) 'thetas map read _______' 3467 3473 3468 variablename = 'thetapwpvg' 3474 variablename = 'thetapwpvg' ! mcw 3469 3475 CALL interpweight_2Dcont(nbpt, 0, 0, lalo, resolution, neighbours, & 3470 3476 contfrac, filename, variablename, lonname, latname, vmin, vmax, nonegative, maskingtype, & … … 3473 3479 WRITE(numout,*) 'thetapwpvg map read _______' 3474 3480 3475 variablename = 'thetafcvg' 3481 variablename = 'thetafcvg' !mcfc 3476 3482 CALL interpweight_2Dcont(nbpt, 0, 0, lalo, resolution, neighbours, & 3477 3483 contfrac, filename, variablename, lonname, latname, vmin, vmax, nonegative, maskingtype, & … … 3495 3501 mcs(:) = mcs_usda(njsc(:)) 3496 3502 ks(:) = ks_usda(njsc(:)) 3497 mcfc(:) = mcf_usda(njsc(:)) 3498 mcw(:) = mcw_usda(njsc(:)) 3503 !!$ mcfc(:) = mcf_usda(njsc(:)) 3504 !!$ mcw(:) = mcw_usda(njsc(:)) 3505 3506 !! Calculation of FC and WP based on above 5 parameters 3507 mvan(:) = un  (un / nvan(:)) 3508 ! Define matrix potential in mm for wilting point and field capacity (with sand vs claysilt variation) 3509 psi_w(:) = 150000. 3510 DO ib=1, nbpt 3511 IF ( ks(ib) .GE. 560 ) THEN ! Sandy soils (560 is equivalent of 2.75 at log scale of Ks, mm/d) 3512 psi_fc(ib) = 1000. 3513 ELSE ! Finer soils 3514 psi_fc(ib) = 3300. 3515 ENDIF 3516 ENDDO 3517 mcfc(:) = mcr(:) + (( mcs(:)  mcr(:)) / (un + ( avan(:) * psi_fc(:))** nvan(:))** mvan(:)) 3518 mcw(:) = mcr(:) + (( mcs(:)  mcr(:)) / (un + ( avan(:) * psi_w(:))** nvan(:))** mvan(:)) 3519 3499 3520 ELSE ! if spmipexp is not among texture or maps or unif 3500 3521 WRITE(numout,*) "Unsupported spmipexp=",spmipexp
Note: See TracChangeset
for help on using the changeset viewer.