diff --git a/scm/src/scm_forcing.F90 b/scm/src/scm_forcing.F90 index 44e6026c..eeb97b68 100644 --- a/scm/src/scm_forcing.F90 +++ b/scm/src/scm_forcing.F90 @@ -61,7 +61,7 @@ subroutine interpolate_forcing(scm_input, scm_state, in_spinup) if(scm_state%model_time >= scm_input%input_time(scm_input%input_ntimes)) then !> - If so, hold the forcing terms constant at the last supplied values. The forcing still needs to be interpolated to the grid. write(*,*) "The model_time has exceeded the specifed period of forcing. Forcing will now be held constant at the last & - specified values." + &specified values." if(scm_state%input_type == 0) then diff --git a/scm/src/scm_input.F90 b/scm/src/scm_input.F90 index 4bb6b3c0..3eb67d15 100644 --- a/scm/src/scm_input.F90 +++ b/scm/src/scm_input.F90 @@ -476,7 +476,7 @@ subroutine get_case_init(scm_state, scm_input) else call check(NF90_INQUIRE_DIMENSION(ncid, varID, tmpName, input_nsoilcat),"nf90_inq_dim(nsoilcat)") end if - + !> - Allocate the dimension variables. allocate(input_pres(input_nlev),input_time(input_ntimes), stat=allocate_status) @@ -732,7 +732,7 @@ subroutine get_case_init(scm_state, scm_input) call check(NF90_CLOSE(NCID=ncid),"nf90_close()") call scm_input%create(input_ntimes, input_nlev, input_nsoil, input_nsnow, input_nice, input_nvegcat, input_nsoilcat) - + ! GJF already done in scm_input%create routine !scm_input%input_nlev = input_nlev !scm_input%input_ntimes = input_ntimes @@ -1072,7 +1072,7 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) real(kind=dp), allocatable :: input_lakedepth(:) !< lake depth (m) real(kind=dp), allocatable :: input_vegtype_frac(:,:) !< fraction of horizontal grid area occupied by given vegetation category real(kind=dp), allocatable :: input_soiltype_frac(:,:) !< fraction of horizontal grid area occupied by given soil category - + real(kind=dp), allocatable :: input_tvxy(:) !< vegetation temperature (K) real(kind=dp), allocatable :: input_tgxy(:) !< ground temperature for Noahmp (K) real(kind=dp), allocatable :: input_tahxy(:) !< canopy air temperature (K) @@ -1188,7 +1188,7 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) logical :: needed_for_lsm_ics, needed_for_model_ics, lev_in_altitude integer :: input_n_init_times, input_n_forcing_times, input_n_lev, input_n_snow, input_n_ice, input_n_soil, input_nvegcat, input_nsoilcat - + missing_value_eps = missing_value + 0.01 !> - Open the case input file found in the processed_case_input dir corresponding to the experiment name. @@ -1215,7 +1215,7 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) write(0,'(a,i0,a)') "The variable 'lev' in the case data file had units different than 'm', 'pa', or 'Pa', but it is expected to be altitude in m or pressure in Pa. Stopping..." STOP end if - + !### TO BE USED IF DEPHY-SCM can be extended to include model ICs ### !possible dimensions (if using model ICs) ierr = NF90_INQ_DIMID(ncid,"nsoil",varID) @@ -1247,8 +1247,8 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) input_nsoilcat = missing_nsoilcat else call check(NF90_INQUIRE_DIMENSION(ncid, varID, tmpName, input_nsoilcat),"nf90_inq_dim(nsoilcat)") - end if - + end if + !> - Allocate the dimension variables. allocate(input_t0 (input_n_init_times), & input_time (input_n_forcing_times), & @@ -1501,14 +1501,14 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) input_sfalb_ice ( input_n_init_times), & input_emis_ice ( input_n_init_times), & stat=allocate_status) - + needed_for_lsm_ics = .False. needed_for_model_ics = .False. if (scm_state%lsm_ics .or. trim(input_surfaceForcingLSM) == "lsm") needed_for_lsm_ics = .True. if (scm_state%model_ics) needed_for_model_ics = .True. - + !> - Read in the initial profiles. - + if (lev_in_altitude) then call NetCDF_read_var(ncid, "pa", .True., input_pres) !zh could be defined in addition to lev, use if so @@ -1532,7 +1532,7 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) end do end if end if - + call NetCDF_read_var(ncid, "ps", .True., input_pres_surf) call NetCDF_read_var(ncid, "ua", .True., input_u) call NetCDF_read_var(ncid, "va", .True., input_v) @@ -1553,7 +1553,7 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) call NetCDF_read_var(ncid, "ri", .False., input_ri) call NetCDF_read_var(ncid, "hur", .False., input_rh) call NetCDF_read_var(ncid, "tke", .False., input_tke) - + call NetCDF_read_var(ncid, "o3", .False., input_ozone) call NetCDF_read_var(ncid, "area", .False., input_area) !orographic parameters @@ -1578,7 +1578,7 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) call NetCDF_read_var(ncid, "lakedepth", needed_for_model_ics, input_lakedepth) call NetCDF_read_var(ncid, "vegtype_frac", needed_for_model_ics, input_vegtype_frac) call NetCDF_read_var(ncid, "soiltype_frac", needed_for_model_ics, input_soiltype_frac) - + !NSST variables call NetCDF_read_var(ncid, "tref", needed_for_model_ics, input_tref) call NetCDF_read_var(ncid, "z_c", needed_for_model_ics, input_z_c) @@ -1598,8 +1598,8 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) call NetCDF_read_var(ncid, "ifd", needed_for_model_ics, input_ifd) call NetCDF_read_var(ncid, "dt_cool", needed_for_model_ics, input_dt_cool) call NetCDF_read_var(ncid, "qrain", needed_for_model_ics, input_qrain) - - + + !> - Allocate the forcing variables. !allocate all, but conditionally read forcing variables given global atts; set unused forcing variables to missing @@ -1666,7 +1666,7 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) end do end do end if - + !conditionally read forcing vars (or set to missing); if the global attribute is set to expect a variable and it doesn't exist, stop the model call NetCDF_conditionally_read_var(adv_u, "adv_ua", "tnua_adv", trim(adjustl(scm_state%case_name))//'.nc', ncid, input_force_u_adv) call NetCDF_conditionally_read_var(adv_v, "adv_va", "tnva_adv", trim(adjustl(scm_state%case_name))//'.nc', ncid, input_force_v_adv) @@ -1735,13 +1735,13 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) ! ! Surface forcing Model LSM ICs ! - + call NetCDF_read_var(ncid, "stc", .False., input_stc) call NetCDF_read_var(ncid, "smc", .False., input_smc) call NetCDF_read_var(ncid, "slc", .False., input_slc) - + call NetCDF_read_var(ncid, "tiice", .False., input_tiice) - + call NetCDF_read_var(ncid, "vegsrc", .False., input_vegsrc ) call NetCDF_read_var(ncid, "vegtyp", .False., input_vegtyp ) call NetCDF_read_var(ncid, "soiltyp", .False., input_soiltyp ) @@ -1779,14 +1779,14 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) call NetCDF_read_var(ncid, "zorll", .False., input_zorll) call NetCDF_read_var(ncid, "zorli", .False., input_zorli) call NetCDF_read_var(ncid, "zorlw", .False., input_zorlw) - + !NoahMP parameters call NetCDF_read_var(ncid, "snicexy", .False., input_snicexy) call NetCDF_read_var(ncid, "snliqxy", .False., input_snliqxy) call NetCDF_read_var(ncid, "tsnoxy", .False., input_tsnoxy ) call NetCDF_read_var(ncid, "smoiseq", .False., input_smoiseq) call NetCDF_read_var(ncid, "zsnsoxy", .False., input_zsnsoxy) - + call NetCDF_read_var(ncid, "tvxy", .False., input_tvxy) call NetCDF_read_var(ncid, "tgxy", .False., input_tgxy) call NetCDF_read_var(ncid, "tahxy", .False., input_tahxy) @@ -1836,12 +1836,12 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) call NetCDF_read_var(ncid, "sfalb_lnd_bck", .False., input_sfalb_lnd_bck) call NetCDF_read_var(ncid, "emis_ice", .False., input_emis_ice) call NetCDF_read_var(ncid, "lai", .False., input_lai) - - + + call check(NF90_CLOSE(NCID=ncid),"nf90_close()") - + call scm_input%create(input_n_forcing_times, input_n_lev, input_n_soil, input_n_snow, input_n_ice, input_nvegcat, input_nsoilcat) - + !fill the scm_input DDT !There may need to be logic to control which of the lon, lat, and init_times to use in the future, but for now, just take the first @@ -2090,7 +2090,7 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) scm_input%input_lakedepth= input_lakedepth(active_init_time) scm_input%input_vegtype_frac = input_vegtype_frac(:,active_init_time) scm_input%input_soiltype_frac = input_soiltype_frac(:,active_init_time) - + scm_input%input_tref = input_tref(active_init_time) scm_input%input_z_c = input_z_c(active_init_time) scm_input%input_c_0 = input_c_0(active_init_time) @@ -2240,7 +2240,7 @@ subroutine get_case_init_DEPHY(scm_state, scm_input) scm_input%input_lh_flux_sfc = input_force_sfc_lat_flx(:) end if end if - + if (trim(input_surfaceForcingLSM) == 'lsm') then !these were considered required variables above, so they should not need to be checked for missing scm_input%input_stc = input_stc(:,active_init_time) @@ -2748,7 +2748,7 @@ subroutine get_reference_profile(scm_state, scm_reference) open(unit=1, file='McCProfiles.dat', status='old', action='read', iostat=ioerror) if(ioerror /= 0) then write(*,*) 'There was an error opening the file McCprofiles.dat in the processed_case_input directory. & - Error code = ',ioerror + &Error code = ',ioerror error stop "There was an error opening the file McCprofiles.dat in the processed_case_input directory." endif diff --git a/scm/src/scm_vgrid.F90 b/scm/src/scm_vgrid.F90 index ca7a4ca6..08905fb7 100644 --- a/scm/src/scm_vgrid.F90 +++ b/scm/src/scm_vgrid.F90 @@ -104,7 +104,7 @@ subroutine get_FV3_vgrid(scm_input, scm_state) open(unit=1, file=scm_state%vert_coord_file, status='old', action='read', iostat=ierr) if(ierr /= 0) then write(*,*) 'There was an error opening the file ', scm_state%vert_coord_file, ' in the run directory. & - Error code = ',ierr + &Error code = ',ierr error stop endif