Not all hard-coded constants in the CLM model are double precision...

**Erik Kluzek < erik > - 2016-07-25 10:56:16 -0600**
**Bugzilla Id:** [2336](http://bugz.cgd.ucar.edu/show_bug.cgi?id=2336)
**Bugzilla Depends:** 2141,
**Bugzilla CC:** andre, dlawren, oleson, rfisher, sacks,

There's a bunch of hard-coded constants in the model that should be double precision, but are single precision...

Here's a sampling...

biogeochem/CNC14DecayMod.F90: spinup_term = 1.
biogeochem/CNDVDriverMod.F90: long_name='land/ocean mask (0.=ocean and 1.=land)')
biogeochem/CNDVDriverMod.F90: long_name='land/ocean mask (0.=ocean and 1.=land)')
biogeochem/CNDVEstablishmentMod.F90: real(r8), parameter :: ramp_agddtw = 300.0
biogeochem/CNDVEstablishmentMod.F90: greffic(p) = 0.
biogeochem/CNDVEstablishmentMod.F90: heatstress(p) = 0.
biogeochem/CNFUNMod.F90: if(delta_CN .gt.0.and. frac_ideal_C_use.lt.1.0)then
biogeochem/CNFUNMod.F90: if(delta_CN.lt.0.0)then
biogeochem/CNFUNMod.F90: do while(exitloop==0.and.cost_retrans_temp .lt. total_n_resistance.and. &
biogeochem/CNFireBaseMod.F90: baf_peatf(c))/2.0*dt))/(dtrotr_col(c)*dayspyr*secspday/dt)/dt
biogeochem/CNGapMortalityMod.F90: real(r8):: k_mort = 0.3 ! coeff of growth efficiency in mortality equation
biogeochem/CNProductsMod.F90: kprod1 = 7.2e-8
biogeochem/CNProductsMod.F90: kprod10 = 7.2e-9
biogeochem/CNProductsMod.F90: kprod100 = 7.2e-10
biogeochem/DUSTMod.F90: write(iulog,*) 'p2l_1d error: sumwt is greater than 1.0 at l= ',index
biogeochem/DUSTMod.F90: write(iulog,*) 'dstmbl: ryn_nbr_frc_thr_prx_opt < 0.03'
biogeochem/DryDepVelocity.F90: if (elai(pi) < (minlai+0.05*(maxlai-minlai))) then
biogeochem/DryDepVelocity.F90: if(ispec.eq.index_o3.or.ispec.eq.index_o3a.or.ispec.eq.index_so2) then
biogeochem/DryDepVelocity.F90: rs=(fsun(pi)*rssun(pi)/elai(pi))+((rssha(pi)/elai(pi))*(1.-fsun(pi)))
biogeochem/DryDepVelocity.F90: if(ispec.ne.index_o3.and.ispec.ne.index_o3a.and.ispec.ne.index_so2) then
biogeochem/NutrientCompetitionCLM45defaultMod.F90: f3 = (2.7/(1.0+exp(-0.004*(annsum_npp(p) - 300.0)))) - 0.4
biogeochem/NutrientCompetitionCLM45defaultMod.F90: f3 = (2.7/(1.0+exp(-0.004*(annsum_npp(p) - 300.0)))) - 0.4
biogeochem/NutrientCompetitionFlexibleCNMod.F90: f3 = (2.7/(1.0+exp(-0.004*(annsum_npp(p) - 300.0)))) - 0.4
biogeochem/NutrientCompetitionFlexibleCNMod.F90: f3 = (2.7/(1.0+exp(-0.004*(annsum_npp(p) - 300.0)))) - 0.4
biogeochem/VOCEmissionMod.F90: co2_ppmv = 1.e6*forc_pco2(g)/forc_pbot(c)
biogeochem/ch4Mod.F90: write(iulog,*)'CH4 Conservation Error in CH4Mod driver, nstep, c, errch4 (gC /m^2.timestep)', &
biogeochem/ch4Mod.F90: write(iulog,*)'CH4 Conservation Error in CH4Mod driver for lake column, nstep, c, errch4 (gC/m^2.timestep)', &
biogeochem/ch4Mod.F90: write(iulog,*)'CH4 Conservation Error in CH4Mod during diffusion, nstep, c, errch4 (mol /m^2.timestep)', &
biogeochem/dynConsBiogeochemMod.F90: if(patch%itype(p)==4.or.patch%itype(p)==6)then
biogeophys/ActiveLayerMod.F90: if ( grc%lat(g) > 0. ) then
biogeophys/ActiveLayerMod.F90: altmax(c) = 0.
biogeophys/ActiveLayerMod.F90: if ( grc%lat(g) <= 0. ) then
biogeophys/ActiveLayerMod.F90: altmax(c) = 0.
biogeophys/BalanceCheckMod.F90: forc_rain_col(c) = 0.
biogeophys/BalanceCheckMod.F90: forc_snow_col(c) = 0.
biogeophys/BandDiagonalMod.F90: ab=0.0
biogeophys/CanopyHydrologyMod.F90: qflx_prec_grnd_snow(p) = 0.
biogeophys/CanopyHydrologyMod.F90: qflx_prec_grnd_rain(p) = 0.
biogeophys/CanopyHydrologyMod.F90: accum_factor=0.1
biogeophys/CanopyHydrologyMod.F90: if (h2osno(c) > 0.0) then
biogeophys/CanopyHydrologyMod.F90: int_snow(c) = 0. !reset prior to adding newsnow below
biogeophys/CanopyHydrologyMod.F90: d=0.0
biogeophys/CanopyHydrologyMod.F90: sigma=1.0e3 * micro_sigma(c) ! convert to mm
biogeophys/CanopyHydrologyMod.F90: +sigma/sqrt(2.0*shr_const_pi)*exp(-d**2/(2.0*sigma**2)) &
biogeophys/CanopyTemperatureMod.F90: qred = (1.-frac_sno(c))*hr_road_perv + frac_sno(c)
biogeophys/LakeFluxesMod.F90: real(r8), parameter :: prn = 0.713 ! Prandtl # for air at neutral stability
biogeophys/LakeFluxesMod.F90: real(r8), parameter :: sch = 0.66 ! Schmidt # for water in air at neutral stability
biogeophys/LakeTemperatureMod.F90: hc_soisno(c) = hc_soisno(c) + cv_lake(c,j)*t_lake(c,j)/1.e6
biogeophys/LakeTemperatureMod.F90: hc_soisno(c) = hc_soisno(c) + cv(c,j)*t_soisno(c,j)/1.e6
biogeophys/LakeTemperatureMod.F90: if (j >= 1) hc_soi(c) = hc_soi(c) + cv(c,j)*t_soisno(c,j)/1.e6
biogeophys/LunaMod.F90: if(tlai(p) > 0.0)then
biogeophys/LunaMod.F90: PARTop = RadTop*4.6 !conversion from w/m2 to umol/m2/s. PAR is still in umol photons, not electrons. Also the par240d_z is only for radiation at visible range. Hence 4.6 not 2.3 multiplier.
biogeophys/LunaMod.F90: if(gpp_day(p)>0.0 )then !only optimize if there is growth and it is C3 plants
biogeophys/LunaMod.F90: write(iulog, *) 'Warning: leaf nitrogen content become unrealistically high (>15.0 g N/m2 leaf) ', &
biogeophys/LunaMod.F90: if(enzs_z(p,z)<1.0) then
biogeophys/LunaMod.F90: if(ndaysteps(p)>0.0) then
biogeophys/LunaMod.F90: if(ndaysteps(p)>0.0)then
biogeophys/LunaMod.F90: chg_per_step = 0.02* FNCa
biogeophys/PhotosynthesisMod.F90: long_name='enzyme decay status during stress: 1.0-fully active; 0.0-all decayed', units='unitless', &
biogeophys/PhotosynthesisMod.F90: if(abs(xm) <= tol1 .or. fb == 0.)then
biogeophys/PhotosynthesisMod.F90: real(r8), parameter :: tolf=1.e-6,toldx=1.e-9
biogeophys/RootBiophysMod.F90: rootfr(p,:) = 0.
biogeophys/RootBiophysMod.F90: real(r8), parameter :: rootprof_exp = 3. ! how steep profile is for root C inputs (1/ e-folding depth) (1/m)
biogeophys/RootBiophysMod.F90: rootfr(p,1) = 0.
biogeophys/SnowHydrologyMod.F90: real(r8), parameter :: min_snow_to_keep = 1.e-9 ! fraction of bottom snow layer to keep with capping
biogeophys/SnowSnicarMod.F90: mu_one = 0.5
biogeophys/SnowSnicarMod.F90: gamma1(i) = (3**0.5)*(2-(omega_star(i)*(1+g_star(i))))/2
biogeophys/SnowSnicarMod.F90: gamma2(i) = omega_star(i)*(3**0.5)*(1-g_star(i))/2
biogeophys/SnowSnicarMod.F90: gamma3(i) = (1-((3**0.5)*g_star(i)*mu_not))/2
biogeophys/SnowSnicarMod.F90: mu_one = 1/(3**0.5)
biogeophys/SnowSnicarMod.F90: gamma3(i) = (1-((3**0.5)*g_star(i)*mu_not))/2
biogeophys/SnowSnicarMod.F90: mu_one = 0.5
biogeophys/SnowSnicarMod.F90: n=floor(i/2.0)
biogeophys/SnowSnicarMod.F90: if (flx_abs_lcl(i,bnd_idx) < -0.00001) then
biogeophys/SnowSnicarMod.F90: write (iulog,"(a,e12.6,a,i6,a,i6)") "SNICAR ERROR: Energy conservation error of : ", energy_sum, &
biogeophys/SnowSnicarMod.F90: if (albout_lcl(bnd_idx) > 1.0) then
biogeophys/SnowSnicarMod.F90: write (iulog,*) "SNICAR ERROR: Albedo > 1.0 at c: ", c_idx, " NSTEP= ",nstep
biogeophys/SoilFluxesMod.F90: if ((h2osoi_liq(c,j)+h2osoi_ice(c,j)) > 0.) then
biogeophys/SoilFluxesMod.F90: qflx_evap_grnd(p) = 0.
biogeophys/SoilFluxesMod.F90: + (1-frac_veg_nosno(p))*(1.-emg(c))*forc_lwrad(c) &
biogeophys/SoilHydrologyInitTimeConstMod.F90: if ( soil_layerstruct /= '10SL_3.5m' )then
biogeophys/SoilHydrologyInitTimeConstMod.F90: d = 0.0
biogeophys/SoilHydrologyInitTimeConstMod.F90: dfdd = exp(-d**2/(2.0*micro_sigma(c)**2))/(micro_sigma(c)*sqrt(2.0*shr_const_pi))
biogeophys/SoilHydrologyInitTimeConstMod.F90: micro_sigma(c)/sqrt(2.0*shr_const_pi)*exp(-d**2/(2.0*micro_sigma(c)**2))
biogeophys/SoilHydrologyInitTimeConstMod.F90: xksat = 0.0070556 *( 10.**(-0.884+0.0153*sandvic(i)) )
biogeophys/SoilHydrologyMod.F90: if (xs(c) > 0.) then
biogeophys/SoilHydrologyMod.F90: if (h2osfc(c) < 0.0) then
biogeophys/SoilHydrologyMod.F90: h2osfc(c) = 0.0
biogeophys/SoilHydrologyMod.F90: * ( 1. - (1.+1.e3*zwt(c)/sucsat(c,nlevsoi))**(-1./bsw(c,nlevsoi)))
biogeophys/SoilHydrologyMod.F90: if(qcharge_tot > 0.) then !rising water table
biogeophys/SoilHydrologyMod.F90: * ( 1. - (1.+1.e3*zwt(c)/sucsat(c,j))**(-1./bsw(c,j)))
biogeophys/SoilHydrologyMod.F90: qcharge_layer=min(qcharge_tot,(s_y*(zwt(c) - zi(c,j-1))*1.e3))
biogeophys/SoilHydrologyMod.F90: if (qcharge_tot <= 0.) exit
biogeophys/SoilHydrologyMod.F90: * ( 1. - (1.+1.e3*zwt(c)/sucsat(c,j))**(-1./bsw(c,j)))
biogeophys/SoilHydrologyMod.F90: qcharge_layer=max(qcharge_tot,-(s_y*(zi(c,j) - zwt(c))*1.e3))
biogeophys/SoilHydrologyMod.F90: if (qcharge_tot >= 0.) then
biogeophys/SoilHydrologyMod.F90: sat_lev=0.9
biogeophys/SoilHydrologyMod.F90: if (rsub_top_tot >= 0.) then
biogeophys/SoilHydrologyMod.F90: sat_lev=0.9
biogeophys/SoilHydrologyMod.F90: if (rsub_top_tot >= 0.) then
biogeophys/SoilHydrologyMod.F90: * ( 1. - (1.+1.e3*zwt(c)/sucsat(c,nlevsoi))**(-1./bsw(c,nlevsoi)))
biogeophys/SoilHydrologyMod.F90: if(rsub_top_tot > 0.) then !rising water table
biogeophys/SoilHydrologyMod.F90: * ( 1. - (1.+1.e3*zwt(c)/sucsat(c,j))**(-1./bsw(c,j)))
biogeophys/SoilHydrologyMod.F90: rsub_top_layer=max(rsub_top_tot,-(s_y*(zi(c,j) - zwt(c))*1.e3))
biogeophys/SoilHydrologyMod.F90: if (rsub_top_tot >= 0.) then
biogeophys/SoilHydrologyMod.F90: sat_lev=0.9
biogeophys/SoilHydrologyMod.F90: * ( 1. - (1.+1.e3*zwt(c)/sucsat(c,k))**(-1./bsw(c,k)))
biogeophys/SoilHydrologyMod.F90: if (drainage_tot >= 0.) then
biogeophys/SoilHydrologyMod.F90: sat_lev=0.9
biogeophys/SoilHydrologyMod.F90: if(rsub_top_tot > 0.) then !rising water table
biogeophys/SoilHydrologyMod.F90: * ( 1. - (1.+1.e3*zwt(c)/sucsat(c,j))**(-1./bsw(c,j)))
biogeophys/SoilHydrologyMod.F90: rsub_top_layer=max(rsub_top_tot,-(s_y*(zi(c,j) - zwt(c))*1.e3))
biogeophys/SoilHydrologyMod.F90: if (rsub_top_tot >= 0.) then
biogeophys/SoilStateInitTimeConstMod.F90: if ( soil_layerstruct /= '10SL_3.5m' )then ! apply soil texture from 10 layer input dataset
biogeophys/SoilStateInitTimeConstMod.F90: soilstate_inst%bsw_col(c,lev) = 2.91 + 0.159*clay
biogeophys/SoilStateInitTimeConstMod.F90: xksat = 0.0070556 *( 10.**(-0.884+0.0153*sand) ) ! mm/s
biogeophys/SoilStateInitTimeConstMod.F90: xksat = 0.0070556 *( 10.**(-0.884+0.0153*sand) ) ! mm/s
biogeophys/SoilTemperatureMod.F90: xmf_h2osfc(c) = 0.
biogeophys/SoilTemperatureMod.F90: zh2osfc=1.0e-3*(0.5*h2osfc(c)) !convert to [m] from [mm]
biogeophys/SoilTemperatureMod.F90: if (h2osoi_ice(c,j) > 0. .AND. t_soisno(c,j) > tfrz) then
biogeophys/SoilTemperatureMod.F90: dzm=(0.5*dz_h2osfc(c)+col%z(c,1))
biogeophys/SoilTemperatureMod.F90: dzm=(0.5*dz_h2osfc(c)+col%z(c,1))
biogeophys/SoilTemperatureMod.F90: dzm=(0.5*dz_h2osfc(c)+col%z(c,1))
biogeophys/SoilTemperatureMod.F90: dzm=(0.5*dz_h2osfc(c)+col%z(c,1))
biogeophys/SoilTemperatureMod.F90: dzm=(0.5*dz_h2osfc(c)+col%z(c,1))
biogeophys/SoilWaterMovementMod.F90: baseflow_sink = 0.
biogeophys/SoilWaterMovementMod.F90: vwc_zwt(c) = min(vwc_zwt(c), 0.5*(watsat(c,nlevsoi) + h2osoi_vol(c,nlevsoi)) )
biogeophys/SoilWaterMovementMod.F90: dwat_temp = 0.
biogeophys/SoilWaterMovementMod.F90: qout(j) = 0.
biogeophys/SoilWaterMovementMod.F90: dqodw1(j) = 0.
biogeophys/SoilWaterMovementMod.F90: dqodw2(j) = 0.
biogeophys/SoilWaterMovementMod.F90: qout(j) = 0.
biogeophys/SoilWaterMovementMod.F90: dqodw1(j) = 0.
biogeophys/SoilWaterMovementMod.F90: qcharge(c) = -ka * (wh_zwt-wh) /((zwt(c)+1.e-3)*m_to_mm)
biogeophys/SoilWaterMovementMod.F90: qcharge(c) = -ka * (wh_zwt-wh)/((zwt(c)-z(c,jwt))*m_to_mm*2.0)
biogeophys/SurfaceAlbedoMod.F90: albgrd_bc(c,ib) = albsod(c,ib)*(1.-frac_sno(c)) + albsnd_bc(c,ib)*frac_sno(c)
biogeophys/SurfaceAlbedoMod.F90: albgri_bc(c,ib) = albsoi(c,ib)*(1.-frac_sno(c)) + albsni_bc(c,ib)*frac_sno(c)
biogeophys/SurfaceAlbedoMod.F90: albgrd_oc(c,ib) = albsod(c,ib)*(1.-frac_sno(c)) + albsnd_oc(c,ib)*frac_sno(c)
biogeophys/SurfaceAlbedoMod.F90: albgri_oc(c,ib) = albsoi(c,ib)*(1.-frac_sno(c)) + albsni_oc(c,ib)*frac_sno(c)
biogeophys/SurfaceAlbedoMod.F90: albgrd_dst(c,ib) = albsod(c,ib)*(1.-frac_sno(c)) + albsnd_dst(c,ib)*frac_sno(c)
biogeophys/SurfaceAlbedoMod.F90: albgri_dst(c,ib) = albsoi(c,ib)*(1.-frac_sno(c)) + albsni_dst(c,ib)*frac_sno(c)
biogeophys/SurfaceAlbedoMod.F90: albgrd_pur(c,ib) = albsod(c,ib)*(1.-frac_sno(c)) + albsnd_pur(c,ib)*frac_sno(c)
biogeophys/SurfaceAlbedoMod.F90: albgri_pur(c,ib) = albsoi(c,ib)*(1.-frac_sno(c)) + albsni_pur(c,ib)*frac_sno(c)
biogeophys/SurfaceAlbedoMod.F90: ((1.-frac_sno(c))*(1-albsod(c,ib))*(flx_absd_snw(c,i,ib)/(1.-albsnd(c,ib))))
biogeophys/SurfaceAlbedoMod.F90: ((1.-frac_sno(c))*(1-albsoi(c,ib))*(flx_absi_snw(c,i,ib)/(1.-albsni(c,ib))))
biogeophys/SurfaceAlbedoMod.F90: ((1.-frac_sno(c))*(1-albsod(c,ib))*(flx_absd_snw(c,i,ib)/(1.-albsnd(c,ib))))
biogeophys/SurfaceAlbedoMod.F90: ((1.-frac_sno(c))*(1-albsoi(c,ib))*(flx_absi_snw(c,i,ib)/(1.-albsni(c,ib))))
biogeophys/SurfaceAlbedoMod.F90: flx_absdv(c,i) = flx_absd_snw(c,i,ib)*(1.-albsnd(c,ib))
biogeophys/SurfaceAlbedoMod.F90: flx_absiv(c,i) = flx_absi_snw(c,i,ib)*(1.-albsni(c,ib))
biogeophys/SurfaceAlbedoMod.F90: flx_absdn(c,i) = flx_absd_snw(c,i,ib)*(1.-albsnd(c,ib))
biogeophys/SurfaceAlbedoMod.F90: flx_absin(c,i) = flx_absi_snw(c,i,ib)*(1.-albsni(c,ib))
biogeophys/SurfaceRadiationMod.F90: write(iulog,*)"flx_absdv1 = ",trd(p,1)*(1.-albgrd(c,1))
biogeophys/SurfaceRadiationMod.F90: write(iulog,*)"flx_absiv1 = ",tri(p,1)*(1.-albgri(c,1))
biogeophys/SurfaceRadiationMod.F90: write(iulog,*)"flx_absdn1 = ",trd(p,2)*(1.-albgrd(c,2))
biogeophys/SurfaceRadiationMod.F90: write(iulog,*)"flx_absin1 = ",tri(p,2)*(1.-albgri(c,2))
biogeophys/SurfaceResistanceMod.F90: d0 = 2.12e-5*(t_soisno(c,1)/273.15)**1.75 ![Bitelli et al., JH, 08]
biogeophys/TemperatureType.F90: this%t_soisno_col(c,j) = 297.56 - (j-1) * ((297.56-293.16)/(nlevgrnd-1))
biogeophys/TemperatureType.F90: this%t_soisno_col(c,1:nlevurb) = 297.56
biogeophys/TemperatureType.F90: this%t_soisno_col(c,j) = 289.46 - (j-1) * ((289.46-295.16)/(nlevgrnd-1))
biogeophys/TemperatureType.F90: this%t_soisno_col(c,1:nlevurb) = 289.46
biogeophys/TemperatureType.F90: this%t_veg_patch(p) = 297.56
biogeophys/TemperatureType.F90: this%t_veg_patch(p) = 289.46
biogeophys/TemperatureType.F90: this%t_ref2m_patch(p) = 297.56
biogeophys/TemperatureType.F90: this%t_ref2m_patch(p) = 289.46
biogeophys/TemperatureType.F90: this%t_ref2m_u_patch(p) = 297.56
biogeophys/TemperatureType.F90: this%t_ref2m_u_patch(p) = 289.46
biogeophys/TemperatureType.F90: this%t_ref2m_r_patch(p) = 297.56
biogeophys/TemperatureType.F90: this%t_ref2m_r_patch(p) = 289.46
biogeophys/UrbanFluxesMod.F90: write(iulog,*)'clm model is stopping - error is greater than 4.e-9 kg/m**2/s'
biogeophys/WaterStateType.F90: avgflag='A', long_name='soil water as frac. of whc for top 0.05 m', &
biogeophys/WaterStateType.F90: if(this%snow_depth_col(c) > 0.0) then
biogeophys/WaterStateType.F90: fmelt = (snowbd/100.)**1.
biogeophys/WaterStateType.F90: this%frac_sno_col(c) = tanh( this%snow_depth_col(c) /(2.5 * zlnd * fmelt) )
cpl/lnd_comp_esmf.F90: "The Community Land Model version 4.0 is " // &
cpl/lnd_comp_esmf.F90: "the land model used in the CESM1.0. " // &
dyn_subgrid/dynLandunitAreaMod.F90: real(r8), parameter :: tol = 1.e-14 ! tolerance for making sure sum of landunit weights equals 1
main/FuncPedotransferMod.F90: bsw = 3.10+0.157*clay-0.003*sand
main/FuncPedotransferMod.F90: xksat = 0.0070556 *(10.**(-0.60+0.0126*sand-0.0064*clay) ) !mm/s now use table 4.
main/FuncPedotransferMod.F90: bsw = 2.91 + 0.159*clay
main/FuncPedotransferMod.F90: xksat = 0.0070556 *( 10.**(-0.884+0.0153*sand) ) ! mm/s, from table 5
main/FuncPedotransferMod.F90: watsat = -0.00108*sand+0.494305
main/FuncPedotransferMod.F90: bsw = 0.137*clay + 3.501
main/clm_driver.F90:1000 format (1x,'nstep = ',i10,' TS = ',f21.15)
main/clm_initializeMod.F90: max_decl = 0.409571
main/clm_varcon.F90: real(r8), parameter :: n_melt=0.7 ! fsca shape parameter
main/clm_varcon.F90: real(r8), parameter :: e_ice=6.0 ! soil ice impedance factor
main/clm_varcon.F90: real(r8), parameter :: pc = 0.4 ! threshold probability
main/clm_varcon.F90: real(r8), parameter :: mu = 0.13889 ! connectivity exponent
main/clm_varcon.F90: real(r8), parameter :: preind_atm_del13c = -6.0 ! preindustrial value for atmospheric del13C
main/clm_varcon.F90: real(r8), parameter :: preind_atm_ratio = SHR_CONST_PDB + (preind_atm_del13c * SHR_CONST_PDB)/1000.0 ! 13C/12C
main/clm_varcon.F90: real(r8) :: c13ratio = preind_atm_ratio/(1.0+preind_atm_ratio) ! 13C/(12+13)C preind atmosphere
main/clm_varcon.F90: real(r8) :: vent_ach = 0.3 ! ventilation rate (air exchanges per hour)
main/clm_varctl.F90: character(len=256), public :: source = "Community Land Model CLM4.0"
main/clm_varctl.F90: character(len=256), public :: conventions = "CF-1.0"
main/clm_varctl.F90: character(len=16), public :: soil_layerstruct = '10SL_3.5m'
main/clm_varpar.F90: if ( soil_layerstruct == '10SL_3.5m' ) then
main/clm_varpar.F90: else if ( soil_layerstruct == '23SL_3.5m' ) then
main/clm_varpar.F90: else if ( soil_layerstruct == '20SL_8.5m' ) then
main/controlMod.F90: write(iulog,*)'ERROR: nlevsno = ',nlevsno,' is not supported, must be in range 3-12.'
main/controlMod.F90: write(iulog,*)'ERROR: h2osno_max = ',h2osno_max,' is not supported, must be greater than 0.0.'
main/controlMod.F90: 'as compared with 0.60, 0.40 for cold frozen lakes with no snow.'
main/histFileMod.F90: long_name='land/ocean mask (0.=ocean and 1.=land)', ncid=nfid(t), &
main/histFileMod.F90: long_name='land/ocean mask (0.=ocean and 1.=land)', ncid=nfid(t), &
main/histFileMod.F90: long_name='pft real/fake mask (0.=fake and 1.=real)', ncid=nfid(t), &
main/histFileMod.F90: long_name='pft real/fake mask (0.=fake and 1.=real)', ncid=nfid(t), &
main/initGridCellsMod.F90: wtcol2lunit = (1. - wtlunit_roof)/3
main/initGridCellsMod.F90: wtcol2lunit = (1. - wtlunit_roof)/3
main/initGridCellsMod.F90: wtcol2lunit = ((1. - wtlunit_roof)/3) * (1.-wtroad_perv)
main/initGridCellsMod.F90: wtcol2lunit = ((1. - wtlunit_roof)/3) * (wtroad_perv)
main/initVerticalMod.F90: real(r8) :: thick_equal = 0.2
main/initVerticalMod.F90: if ( soil_layerstruct == '10SL_3.5m' ) then
main/initVerticalMod.F90: else if ( soil_layerstruct == '23SL_3.5m' )then
main/initVerticalMod.F90: zsoi(j) = 0.5*(zisoi(j-1) + zisoi(j))
main/initVerticalMod.F90: else if ( soil_layerstruct == '20SL_8.5m' ) then
main/initVerticalMod.F90: zsoi(j) = 0.5*(zisoi(j-1) + zisoi(j))
main/initVerticalMod.F90: dzsoi_decomp(1) = 1.
main/initVerticalMod.F90: ziurb_wall(l,0) = 0.
main/initVerticalMod.F90: ziurb_roof(l,0) = 0.
main/initVerticalMod.F90: ziurb_wall(l,0) = 0.
main/initVerticalMod.F90: ziurb_roof(l,0) = 0.
main/initVerticalMod.F90: zurb_wall(l,j) = (j-0.5)*(thick_wall(l)/float(nlevurb)) !node depths
main/initVerticalMod.F90: zurb_roof(l,j) = (j-0.5)*(thick_roof(l)/float(nlevurb)) !node depths
main/initVerticalMod.F90: dzurb_roof(l,1) = 0.5*(zurb_roof(l,1)+zurb_roof(l,2)) !thickness b/n two interfaces
main/initVerticalMod.F90: dzurb_roof(l,j)= 0.5*(zurb_roof(l,j+1)-zurb_roof(l,j-1))
main/initVerticalMod.F90: dzurb_wall(l,1) = 0.5*(zurb_wall(l,1)+zurb_wall(l,2)) !thickness b/n two interfaces
main/initVerticalMod.F90: dzurb_wall(l,j)= 0.5*(zurb_wall(l,j+1)-zurb_wall(l,j-1))
main/initVerticalMod.F90: ziurb_wall(l,0) = 0.
main/initVerticalMod.F90: ziurb_wall(l,j) = 0.5*(zurb_wall(l,j)+zurb_wall(l,j+1)) !interface depths
main/initVerticalMod.F90: ziurb_wall(l,nlevurb) = zurb_wall(l,nlevurb) + 0.5*dzurb_wall(l,nlevurb)
main/initVerticalMod.F90: ziurb_roof(l,0) = 0.
main/initVerticalMod.F90: ziurb_roof(l,j) = 0.5*(zurb_roof(l,j)+zurb_roof(l,j+1)) !interface depths
main/initVerticalMod.F90: ziurb_roof(l,nlevurb) = zurb_roof(l,nlevurb) + 0.5*dzurb_roof(l,nlevurb)
main/subgridAveMod.F90: write(iulog,*)'p2c_1d error: sumwt is greater than 1.0'
main/subgridAveMod.F90: write(iulog,*)'p2c_2d error: sumwt is greater than 1.0 at c= ',index,' lev= ',j
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: write(iulog,*)'p2l_1d error: sumwt is greater than 1.0 at l= ',index
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: write(iulog,*)'p2l_2d error: sumwt is greater than 1.0 at l= ',index,' j= ',j
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: write(iulog,*)'p2g_1d error: sumwt is greater than 1.0 at g= ',index
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: write(iulog,*)'p2g_2d error: sumwt gt 1.0 at g/sumwt = ',index,sumwt(index)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: write(iulog,*)'c2l_1d error: sumwt is greater than 1.0 at l= ',index
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: write(iulog,*)'c2l_2d error: sumwt is greater than 1.0 at l= ',index,' lev= ',j
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: write(iulog,*)'c2g_1d error: sumwt is greater than 1.0 at g= ',index
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 * lun%canyon_hwr(l)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = (3.0 * lun%canyon_hwr(l)) / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: scale_c2l(c) = 3.0 / (2.*lun%canyon_hwr(l) + 1.)
main/subgridAveMod.F90: write(iulog,*)'c2g_2d error: sumwt is greater than 1.0 at g= ',index
main/subgridAveMod.F90: write(iulog,*)'l2g_1d error: sumwt is greater than 1.0 at g= ',index
main/subgridAveMod.F90: write(iulog,*)'l2g_2d error: sumwt is greater than 1.0 at g= ',index,' lev= ',j
main/subgridRestMod.F90: real(r8), parameter :: tol = 5.e-3 ! tolerance for checking weights
main/surfrdMod.F90: if (ldomain%lonc(n)-surfdata_domain%lonc(n) > 300.) then
main/surfrdMod.F90: elseif (ldomain%lonc(n)-surfdata_domain%lonc(n) < -300.) then
main/surfrdUtilsMod.F90: write(iulog,*) trim(caller), ' ERROR: sum of ', trim(name), ' not 1.0 at nl=', nindx
soilbiogeochem/SoilBiogeochemCarbonStateType.F90: if (zsoi(j) < 0.3 ) then !! only initialize upper soil column
soilbiogeochem/SoilBiogeochemCarbonStateType.F90: m = 1. / decomp_cascade_con%spinup_factor(k)
soilbiogeochem/SoilBiogeochemCompetitionMod.F90: sminn_tot(c) = 0.
soilbiogeochem/SoilBiogeochemCompetitionMod.F90: if (sminn_tot(c) > 0.) then
soilbiogeochem/SoilBiogeochemCompetitionMod.F90: sminn_tot(c) = 0.
soilbiogeochem/SoilBiogeochemCompetitionMod.F90: if (sminn_tot(c) > 0.) then
soilbiogeochem/SoilBiogeochemDecompCascadeBGCMod.F90: tau_l1 = 1./18.5
soilbiogeochem/SoilBiogeochemDecompCascadeBGCMod.F90: tau_l2_l3 = 1./4.9
soilbiogeochem/SoilBiogeochemDecompCascadeBGCMod.F90: tau_s1 = 1./7.3
soilbiogeochem/SoilBiogeochemDecompCascadeBGCMod.F90: tau_s2 = 1./0.2
soilbiogeochem/SoilBiogeochemDecompCascadeBGCMod.F90: tau_s3 = 1./.0045
soilbiogeochem/SoilBiogeochemDecompCascadeBGCMod.F90: tau_cwd = 1./0.3
soilbiogeochem/SoilBiogeochemDecompCascadeCNMod.F90: rf_decomp_cascade(bounds%begc:bounds%endc,1:nlevdecomp,i_s4atm) = 1.
soilbiogeochem/SoilBiogeochemLittVertTranspMod.F90: epsilon = 1.e-30
soilbiogeochem/SoilBiogeochemLittVertTranspMod.F90: trcr_tendency_ptr(c,j,s) = 0.-(conc_trcr(c,j) + source(c,j,s))
soilbiogeochem/SoilBiogeochemNLeachingMod.F90: real(r8), parameter :: depth_runoff_Nloss = 0.05 ! (m) depth over which runoff mixes with soil water for N loss to runoff
soilbiogeochem/SoilBiogeochemNitrifDenitrifMod.F90: pH(bounds%begc:bounds%endc) = 6.5 !!! set all soils with the same pH as placeholder here
soilbiogeochem/SoilBiogeochemNitrifDenitrifMod.F90: k_nitr_ph_vr(c,j) = 0.56 + atan(rpi * 0.45 * (-5.+ pH(c)))/rpi
soilbiogeochem/SoilBiogeochemNitrifDenitrifMod.F90: n2_n2o_ratio_denit_vr(c,j) = max(0.16*ratio_k1(c,j), ratio_k1(c,j)*exp(-0.8 * ratio_no3_co2(c,j))) * fr_WFPS(c,j)
soilbiogeochem/SoilBiogeochemNitrogenStateType.F90: m = 1. / decomp_cascade_con%spinup_factor(k)
soilbiogeochem/SoilBiogeochemNitrogenUptakeMod.F90: sminn_tot(c) = 0.
soilbiogeochem/SoilBiogeochemNitrogenUptakeMod.F90: if (sminn_tot(c) > 0.) then
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: real(r8), public :: rootprof_exp = 3. ! how steep profile is for root C inputs (1/ e-folding depth) (1/m)
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: real(r8), public :: surfprof_exp = 10. ! how steep profile is for surface components (1/ e_folding depth) (1/m)
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: real(r8) :: delta = 1.e-10
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: cinput_rootfr(p,1) = 0.
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: froot_prof(p,1) = 1./dzsoi_decomp(1)
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: croot_prof(p,1) = 1./dzsoi_decomp(1)
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: leaf_prof(p,1) = 1./dzsoi_decomp(1)
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: stem_prof(p,1) = 1./dzsoi_decomp(1)
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: nfixation_prof(c,1) = 1./dzsoi_decomp(1)
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: ndep_prof(c,1) = 1./dzsoi_decomp(1)
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: ndep_prof_sum = 0.
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: nfixation_prof_sum = 0.
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: froot_prof_sum = 0.
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: croot_prof_sum = 0.
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: leaf_prof_sum = 0.
soilbiogeochem/SoilBiogeochemVerticalProfileMod.F90: stem_prof_sum = 0.
utils/clm_time_manager.F90: if ( (get_curr_calday > 366.0) .and. (get_curr_calday <= 367.0) .and. &
utils/clm_time_manager.F90: if ( (get_curr_calday < 1.0) .or. (get_curr_calday > 366.0) )then
utils/clm_time_manager.F90: if ( (get_calday > 366.0) .and. (get_calday <= 367.0) .and. &
utils/clm_time_manager.F90: if ( (get_calday < 1.0) .or. (get_calday > 366.0) )then
utils/domainMod.F90: domain%frac = -1.0e36