Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Process pm #44

Merged
merged 3 commits into from
Sep 28, 2023
Merged

Process pm #44

merged 3 commits into from
Sep 28, 2023

Conversation

Ruifang-Li
Copy link
Contributor

This is program to read in PM observation from ASCII file and convert to PrepBUFR format for PM DA.
Test with real time runs on Hera is successful.

PRIVATE ${PROJECT_NAME}::esggrid_util
PRIVATE ${PROJECT_NAME}::pesglib
PRIVATE NetCDF::NetCDF_Fortran
PRIVATE gsi::gsi
Copy link
Contributor

Choose a reason for hiding this comment

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

May not need ncio,esggird_util,pesglib, gsi

Copy link
Contributor

Choose a reason for hiding this comment

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

You can delete them and try again.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

done

! trim(adjustl(pm_aqi_str(i))), trim(adjustl(pm_measured_str(i))), trim(adjustl(pm_str(i)))
!write(*,'(a10,7i15)') 'len ', len(trim(adjustl(sta_id(i)))), len(trim(adjustl(lat_str(i)))), len(trim(adjustl(lon_str(i)))), len(trim(adjustl(elv_str(i)))), &
! len(trim(adjustl(pm_aqi_str(i)))), len(trim(adjustl(pm_measured_str(i)))), len(trim(adjustl(pm_str(i))))
end do
Copy link
Contributor

Choose a reason for hiding this comment

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

remove line 61-64 if they are not needed.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

done


! hdstr='SID XOB YOB DHR TYP T29 SQN PROCN RPT CAT TYPO TSIG'
! obstr='TPHR QCIND COPOPM'

Copy link
Contributor

Choose a reason for hiding this comment

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

remove line 76, 83,84 if they are not needed.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

done

str=adjustl(sta_id(j))
c_sid=str(2:9)
!write(*,'(2a15,4f15.5)') 'Obs ', c_sid, lat, lon, elv, pm
hdr(1)=rstation_id; hdr(2)=lon; hdr(3)=lat; hdr(4)=0.0_8; hdr(5)=102; hdr(10)=6 ! Single level report
Copy link
Contributor

Choose a reason for hiding this comment

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

remove line 109

Copy link
Contributor Author

Choose a reason for hiding this comment

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

done

@Ruifang-Li Ruifang-Li merged commit c5dbfb2 into NOAA-GSL:develop Sep 28, 2023
1 check passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants