x.f90

来自「CLM集合卡曼滤波数据同化算法」· F90 代码 · 共 24 行

!
! -----------------------------------------------------
            print*, 'sand = ?, clay = ?'
            read(*,*) sand, clay
             porsl = 0.489 - 0.00126*sand
              phi0 = 10. * ( 10.**(1.88-0.0131*sand) )
               bsw = 2.91 + 0.159*clay

                 dkm    = (8.80*sand+2.92*clay) / (sand+clay)       ! W/(m K)
              csol = (2.128*sand+2.385*clay)/(sand+clay)*1.e6  ! J/(m3 K)

                     bd = (1.- porsl)*2.7e3
              dkmg = dkm**(1.-porsl)
            dksatu = dkmg*0.57**porsl
             dkdry = (.135*bd + 64.7) / (2.7e3 - 0.947*bd)

                  xksat = 0.0070556 * ( 10.**(-0.884+0.0153*sand) ) ! mm/s

             print*, 'porsl, phi0, bsw, xksat', porsl, phi0, bsw, xksat
             print*, 0.01**(-bsw)

! -----------------------------------------------------
              end