Skip to content

Conversation

@ADCollard
Copy link

Description

Changes required to use netCDF files by default. Added missing routines for antenna correction coefficients.

Issue(s) addressed

Partially Resolves #192

Dependencies

None

Impact

Expected impact on downstream repositories: None

Checklist

  • I have performed a self-review of my own code
  • I have made corresponding changes to the documentation
  • I have run the unit tests before creating the PR

@ADCollard
Copy link
Author

Putting this Draft PR in for initial review. It is still a draft PR as:

  1. I am still getting seg faults when running the code from the GSI
  2. The changes to src/Coefficients/FitCoeff/FitCoeff_Define.f90 that were reversed are still in this version as we still need to address this.

Right now inclusion of CRTMv3 with GFSv17 is unlikely given the issues with seg faults.

@BenjaminTJohnson
Copy link
Contributor

@ADCollard I was in Italy last week, so just now getting back to this today.

(RTV(1)%n_Stokes > 1 .AND. CloudC%N_PHASE_ELEMENTS < 6 )) THEN
Error_Status = FAILURE
WRITE( Message,'("N_PHASE_ELEMENTS OF CLOUD LUT NOT RIGHT ",i0)' ) CloudC%N_PHASE_ELEMENTS
WRITE( Message,'("N_PHASE_ELEMENTS OF CLOUD LUT NOT RIGHT ",i0)' ) CloudC%N_PHASE_ELEMENTS, AeroC%N_PHASE_ELEMENTS, RTV(1)%n_Stokes
Copy link
Contributor

Choose a reason for hiding this comment

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

It makes sense to add RTV(1)%n_Stokes, but there is no need to print AeroC%N_PHASE_ELEMENTS here.

(RTV(1)%n_Stokes > 1 .AND. AeroC%N_PHASE_ELEMENTS < 6 )) THEN
Error_Status = FAILURE
WRITE( Message,'("N_PHASE_ELEMENTS OF AEROSOL LUT NOT RIGHT ",i0)' ) AeroC%N_PHASE_ELEMENTS
WRITE( Message,'("N_PHASE_ELEMENTS OF AEROSOL LUT NOT RIGHT ",i0)' ) AeroC%N_PHASE_ELEMENTS, AeroC%N_PHASE_ELEMENTS, RTV(1)%n_Stokes
Copy link
Contributor

Choose a reason for hiding this comment

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

Same here. It makes sense to add RTV(1)%n_Stokes, but there is no need to print CloudC%N_PHASE_ELEMENTS here. (AeroC%N_PHASE_ELEMENTS is added twice, likely a typo)

The IF statement at line 618 makes sure CloudC%N_PHASE_ELEMENTS == AeroC%N_PHASE_ELEMENTS if they are both loaded.

! ...Default filenames
Default_Aerosol_Model = 'CRTM'
Default_AerosolCoeff_File = 'AerosolCoeff.bin'
Default_AerosolCoeff_File = 'AerosolCoeff.nc'
Copy link
Contributor

Choose a reason for hiding this comment

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

Haven't tested this branch yet, waiting for fix folder updates on coefficient filenames (nc4 --> nc)

Copy link
Contributor

Choose a reason for hiding this comment

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

@chengdang working on it now -- not ready yet, so ignore pushes to this PR until I've given the thumbs up.

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.

Some bug fixes to ensure CRTMv3 can run with netCDF files

4 participants