Skip to content
This repository was archived by the owner on May 5, 2026. It is now read-only.

Add station IDs with txxxxyyyy for processed IMS snow DA#1

Merged
CoryMartin-NOAA merged 5 commits into
NOAA-EMC:developfrom
YoulongXia-NOAA:feature/addIDs
May 7, 2025
Merged

Add station IDs with txxxxyyyy for processed IMS snow DA#1
CoryMartin-NOAA merged 5 commits into
NOAA-EMC:developfrom
YoulongXia-NOAA:feature/addIDs

Conversation

@YoulongXia-NOAA
Copy link
Copy Markdown

PR Instructions:

  1. Provide details under all headings below.
  2. Assign Clara and one other person as reviewers.
  3. If the PR is not ready for merging, add the "DRAFT/DO NOT MERGE" label.
  4. When a PR is ready to merge, remove the "DRAFT/DO NOT MERGE" and email Clara.
  5. Before requesting that the PR be merged, complete the checklist below.

Notes on preparing PR, using git can be found in README_git

Describe your changes

Summarise all code changes included in PR:

List any associated PRs in the submodules.

Issue ticket number and link

List the git Issue that this PR addresses:

Test output

Is this PR expected to pass the DA_IMS_test (ie., does it change the output)?

Does it pass the DA_IMS_test?

If changes to the test results are expected, what are these changes? Provide a link to the output directory when running the test:

Checklist before requesting a review

  • My branch being merged is up to date with the latest develop.
  • I have performed a self-review of my code by examining the differences that will be merged.
  • I have not made any unnecessary code changes / changed any default behavior.
  • My code passes the DA_IMS_test, or differences can be explained.

Comment thread sorc/IMSaggregate_mod.f90 Outdated
Comment on lines +175 to +187
if(i<=9) then
write(istr1,'(i1)') i
xind='000'//istr1
else if (i >= 10 .AND. i <= 99) then
write(istr2,'(i2)') i
xind='00'//istr2
else if (i >= 100 .AND. i <= 999) then
write(istr3,'(i3)') i
xind='0'//istr3
else
write(istr4,'(i4)') i
xind=istr4
end if
Copy link
Copy Markdown

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

can't you specify leading zeros in the write statement format?

Copy link
Copy Markdown

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

'(I4.4)'

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Corrected as suggested. I test both nc and ascii input. Thank you for simplifying the code.

Comment thread sorc/IMSaggregate_mod.f90 Outdated
character(len=3) :: jstr3,istr3
character(len=4) :: xind, yind, jstr4,istr4
character(len=9) :: staids(idim,jdim,6) ! station ID string
integer :: sid(idim,jdim,6) ! staiion ID integer
Copy link
Copy Markdown

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
integer :: sid(idim,jdim,6) ! staiion ID integer
integer :: sid(idim,jdim,6) ! station ID integer

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Corrected.

Comment thread sorc/IMSaggregate_mod.f90 Outdated
character(len=2) :: jstr2,istr2
character(len=3) :: jstr3,istr3
character(len=4) :: xind, yind, jstr4,istr4
character(len=9) :: staids(idim,jdim,6) ! station ID string
Copy link
Copy Markdown

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggest to use standard name stid to replace staids

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Corrected as suggested.

Comment thread sorc/IMSaggregate_mod.f90
Comment on lines +173 to +177
write(istr4,'(i4.4)') i
xind=istr4
! create yindex string for final station identification
write(jstr4,'(i4.4)') j
yind=jstr4
Copy link
Copy Markdown

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
write(istr4,'(i4.4)') i
xind=istr4
! create yindex string for final station identification
write(jstr4,'(i4.4)') j
yind=jstr4
write(istr4,'(i4.4)') i
xind=istr4
! create yindex string for final station identification
write(jstr4,'(i4.4)') j
yind=jstr4

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Corrected.

Comment thread sorc/IMSaggregate_mod.f90 Outdated
do i=1,idim
do j=1,jdim
write(tile_str, '(i1)') t ! assuming <10 tiles.
! create xindex string for final station identification
Copy link
Copy Markdown

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
! create xindex string for final station identification
! create xindex string for final station identification

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Corrected.

Copy link
Copy Markdown

@CoryMartin-NOAA CoryMartin-NOAA left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants