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.

Sign up for GitHub

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

Feature/add_CrISobserror #1

Open
wants to merge 12 commits into
base: main
Choose a base branch
from
Open

Feature/add_CrISobserror #1

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
12 commits
Select commit Hold shift + click to select a range
6d65ac6
Update README.md
ibanos90 Jul 12, 2022
67ae00d
Modify ahi reading and bugs to write the ioda files
ibanos90 Jul 19, 2022
aa50c6e
Fix reading of data for radiance calculation
ibanos90 Aug 6, 2022
7e576a2
Use '&' when reading the variable for better code readability. Add so…
ibanos90 Aug 6, 2022
5afca03
Fixing bug from previous commit. The reading of idata was deleted uni…
ibanos90 Aug 8, 2022
1bdf499
Add superobbing for AHI data
ibanos90 Aug 8, 2022
048bd5f
Merge pull request #3 from NCAR/feature/ahiFix
ibanos90 Aug 12, 2022
f44308d
Delete unnecessary print
ibanos90 Aug 12, 2022
a3d537b
Merge branch 'main' into feature/ahiSuperob
ibanos90 Aug 12, 2022
fc8b074
Remove debugging build flags
ibanos90 Aug 12, 2022
8797470
Merge pull request #4 from NCAR/feature/ahiSuperob
ibanos90 Aug 15, 2022
dbf2a3e
Adding obserror for CrIS
zhumingying Jul 18, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
5 changes: 5 additions & 0 deletions README.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -29,3 +29,8 @@ File types that obs2ioda can handle:
* https://rda.ucar.edu/datasets/ds735.0/
* https://registry.opendata.aws/noaa-goes/
* https://registry.opendata.aws/noaa-himawari/

## Contributors
* Primary original developer: Jamie Bresch (@jamiebresch)

Original repository can be found at https://github.com/jamiebresch/obs2ioda
2 changes: 1 addition & 1 deletion goes_abi/src/Makefile
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@
#INTEL compiler
#--------------
FC = ifort
FFLAGS = -mcmodel medium -g -traceback -debug all -check all # needed for intel error message "failed to convert GOTPCREL relocation"
FFLAGS = -mcmodel medium # needed for intel error message "failed to convert GOTPCREL relocation" #-g -traceback -debug all -check all

LIBS = -L$(NETCDF)/lib -lnetcdff -lnetcdf
INCS = -I$(NETCDF)/include
Expand Down
4 changes: 2 additions & 2 deletions obs2ioda-v2/src/Makefile
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -14,8 +14,8 @@
FC = ifort
#BUFR_LIB = -L/glade/u/home/hclin/extlib/intel -lbufr
#FFLAGS = -mcmodel medium # needed for intel error message "failed to convert GOTPCREL relocation"
BUFR_LIB = -L/glade/p/mmm/parc/ivette/pandac/converters/WRFDA_3DVAR_dmpar/var/external/bufr -lbufr #-L/glade/u/home/hclin/extlib/intel -lbufr
FFLAGS = -mcmodel medium -g -traceback -debug all -check all # needed for intel error message "failed to convert GOTPCREL relocation"
BUFR_LIB = -L/glade/campaign/mmm/parc/ivette/pandac/converters/WRFDA_3DVAR_dmpar/var/external/bufr -lbufr
FFLAGS = -mcmodel medium # needed for intel error message "failed to convert GOTPCREL relocation" # -g -traceback -debug all -check all

LIBS = -L$(NETCDF)/lib -lnetcdff -lnetcdf ${BUFR_LIB}
INCS = -I$(NETCDF)/include
Expand Down
30 changes: 30 additions & 0 deletions obs2ioda-v2/src/define_mod.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -412,10 +412,40 @@ subroutine set_brit_obserr(name_inst, nchan, obserrors)
case default
return
end select

else if ( name_inst(1:4) == 'cris' ) then
select case ( trim(name_inst) )
case ( 'cris_npp' )
obserrors = (/ 1.0, 0.7, 0.7, 0.7, 0.7, 1.359, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 1.0, 0.6, 1.0, 0.6, 1.0, 0.5, 1.0, 0.5, 1.0, 0.6, 1.0, 0.5, 1.0, 0.5, 0.756, 0.5, 1.0, 0.5, 1.0, 0.5, 1.0, 0.5, 1.0, 0.5, 1.0, 0.5, 1.0, 0.5, 1.0, 0.5, 0.6, 0.5, 1.0, 0.5, 1.0, 0.5, 1.0, 0.5, 0.6, 0.5, 1.0, 0.45, 1.0, 0.45, 1.0, 0.45, 0.635, 0.45, 1.0, 0.45, 1.0, 0.45, 0.735, 0.45, 0.878, 0.45, 0.696, 0.4, 2.0, 0.4, 1.0, 0.4, 1.0, 0.4, 1.0, 0.4, 0.6, 0.35, 1.0, 0.35, 0.701, 0.35, 1.0, 0.35, 0.6, 0.35, 0.663, 0.35, 1.0, 0.35, 1.083, 0.35, 0.6, 0.35, 1.0, 0.35, 0.6, 0.35, 0.6, 0.35, 1.0, 0.3, 0.6, 0.3, 0.6, 0.3, 1.0, 0.3, 0.6, 0.3, 0.6, 0.3, 1.0, 0.3, 0.6, 0.3, 1.0, 0.3, 0.773, 0.3, 0.6, 0.6, 0.6, 0.3, 0.813, 0.907, 0.802, 0.3, 1.493, 1.0, 0.856, 0.3, 1.0, 0.6, 1.0, 0.3, 1.0, 1.0, 1.0, 0.3, 1.0, 1.0, 1.0, 0.3, 1.0, 1.0, 1.0, 0.3, 1.0, 1.0, 0.3, 1.0, 1.0, 1.0, 0.3, 1.0, 1.0, 0.3, 2.0, 1.0, 0.3, 1.0, 0.3, 1.0, 0.3, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 0.3, 1.0, 1.0, 0.3, 1.0, 1.0, 1.0, 0.3, 0.3, 0.3, 0.5, 1.0, 1.0, 2.0, 1.0, 0.5, 0.5, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 0.5, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 0.5, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 0.5, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 0.5, 0.5, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 0.5, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /)
case ( 'cris_n20' )
obserrors = (/ 1.227, 1.201, 1.279, 1.333, 1.31, 1.645, 1.265, 1.129, 1.016, 1.023, 0.976, 1.009, 0.963, 0.96, 0.905, 0.93, 0.886, 0.886, 0.883, 0.873, 0.858, 0.4980026, 0.5111945, 0.4921504, 0.495887, 0.4953809, 0.4836317, 0.5139858, 0.5005143, 0.4916807, 0.4880601, 0.4656278, 0.4793087, 0.4637685, 0.4556648, 0.4665926, 0.4467616, 0.4533513, 0.4471499, 0.4448422, 0.4468595, 0.89, 0.4425544, 0.438781, 0.897, 0.4368367, 0.438022, 0.766, 0.4378737, 0.801, 0.4404045, 0.4405309, 0.4409174, 0.4472441, 0.4555138, 0.4433328, 0.443671, 0.4453777, 0.4447534, 0.4465382, 0.4498734, 0.4488066, 0.69, 0.4533533, 0.4471555, 0.4550021, 0.4562328, 0.4519525, 0.4639232, 0.699, 0.4573326, 0.4603723, 0.4533107, 0.4692493, 0.839, 0.4456794, 0.4457273, 0.5153976, 0.5083739, 1.125, 1.082, 0.958, 0.823, 0.788, 0.4950289, 0.480949, 0.4732342, 0.4860787, 0.463246, 0.762, 0.4747573, 0.5007465, 0.5711111, 0.873, 0.5469363, 0.995, 0.776, 0.774, 0.74, 0.733, 0.4694708, 0.682, 0.695, 0.658, 0.712, 0.4680831, 0.701, 0.713, 0.704, 0.746, 0.4805201, 1.157, 1.138, 0.741, 0.843, 0.802, 0.6293901, 0.5885451, 0.789, 0.841, 0.785, 0.951, 0.562651, 0.5675052, 0.94, 0.5166174, 0.95, 0.5291152, 0.589188, 0.5976298, 0.5834491, 0.6512306, 0.6747992, 0.6614965, 0.6003342, 0.566927, 0.558674, 0.5507283, 0.5870586, 1.188, 1.313, 1.315, 1.358, 1.441, 1.39, 1.641, 1.593, 1.484, 1.447, 1.407, 1.503, 1.42, 1.481, 1.457, 1.495, 1.657, 1.626, 1.613, 1.525, 1.497, 1.56, 1.569, 1.596, 1.574, 1.592, 1.719, 1.624, 2.007, 1.672, 1.895, 1.671, 1.938, 2.025, 1.564, 1.769, 1.506, 1.594, 1.708, 1.527, 1.53, 1.523, 1.531, 1.522, 1.654, 1.525, 1.488, 1.47, 1.611, 1.374, 1.441, 0.619922, 0.6223155, 0.6035748, 0.6003346, 0.5991098, 0.5979717, 0.5910122, 0.5764011, 1.027, 0.5593015, 1.19, 1.001, 1.002, 1.02, 1.021, 1.028, 1.013, 1.036, 1.013, 1.038, 0.996, 1.006, 0.5400555, 0.5500141, 0.5575057, 1.023, 0.5635096, 0.5786099, 0.5807203, 1.038, 1.021, 1.005, 1.006, 0.996, 1.022, 1.087, 1.01, 1.059, 1.008, 1.016, 0.968, 1.007, 1.081, 1.564, 1.69, 1.140175, 1.76, 1.712, 1.705, 1.77, 1.689, 1.675, 1.568, 1.574, 1.722, 1.732, 1.140175, 1.077, 1.845, 1.140175, 1.995, 1.922, 1.140175, 0.999, 1.011, 0.95, 1.006, 1.022, 1.023, 1.022, 1.119, 1.071, 1.083, 1.099, 0.8646028, 1.208, 1.354, 1.019, 1.79, 1.1, 1.153, 1.052, 1.037, 1.347, 1.393, 1.206, 1.252, 1.559, 1.484, 1.854, 1.882, 0.9402268, 1.005375, 2.884, 2.561, 2.485, 2.309, 0.9703129, 2.782, 0.9152867, 3.097, 3.023, 3.149, 3.141, 3.18, 3.023, 0.9914418, 3.022, 3.17, 3.077, 1.281697, 3.066, 3.404, 3.103, 3.128, 3.303, 3.362, 1.118811, 3.277, 3.007, 3.297, 3.344, 3.333, 2.912, 2.861, 2.993, 3.345, 3.282, 3.133, 3.236, 2.882, 3.124, 3.65, 3.623, 1.063818, 3.41, 3.319, 3.44, 3.243, 3.4, 3.013, 2.735, 3.609, 3.416, 3.354, 3.125, 3.174, 2.881, 3.553, 3.471, 3.685, 3.225, 3.372, 3.196, 3.464, 3.331, 3.397, 3.406, 3.232, 3.357, 3.388, 3.377, 3.584, 3.485, 3.493, 3.47, 3.272, 3.279, 3.459, 3.446, 4.128, 3.786, 3.752, 3.612, 3.871, 3.727, 3.045, 2.929, 2.856, 2.845, 2.858, 3.996, 3.922, 4.074, 4.014, 4.019, 3.912, 3.924, 3.777, 3.588, 3.277, 2.928, 2.541, 2.245, 1.964, 1.743, 1.573, 1.511, 1.447, 1.403, 1.342, 1.285, 1.226, 1.173, 1.133, 1.086, 1.07, 1.018, 0.99, 0.981, 0.985, 0.97, 0.972, 0.958, 0.962, 0.981, 0.985, 0.977, 0.967, 0.942, 0.947, 0.966, 0.981, 0.98, 0.965, 0.969, 0.976, 0.973, 0.976, 0.97, 1.335, 1.394, 1.419, 1.423, 1.434, 1.451, 1.424, 1.471, 1.454, 1.481, 1.482 /)
case default
return
end select

else
return
end if

end subroutine set_brit_obserr

subroutine set_ahi_obserr(name_inst, nchan, obserrors)
implicit none

character(len=*), intent(in) :: name_inst ! instrument name
integer(i_kind), intent(in) :: nchan ! channel number
real(r_kind), intent(out) :: obserrors(nchan)
obserrors(:) = missing_r
if ( name_inst(1:3) == 'ahi' ) then
select case ( trim(name_inst) )
case ( 'ahi_himawari8' )
obserrors = (/ 2.2, 3.0, 2.5, 2.2, 2.2, 2.2, 2.2, 2.2, 2.2, 2.2 /)
case default
return
end select
else
return
end if
end subroutine set_ahi_obserr

end module define_mod
Loading