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

Create NetCDF module for writing standalone CRTM outputs #62

Closed
chengdang opened this issue Aug 7, 2023 · 4 comments
Closed

Create NetCDF module for writing standalone CRTM outputs #62

chengdang opened this issue Aug 7, 2023 · 4 comments
Assignees
Labels
CRTM CRTM enhancement New feature or request

Comments

@chengdang
Copy link
Contributor

chengdang commented Aug 7, 2023

Create NetCDF module for writing standalone CRTM outputs

@chengdang chengdang added enhancement New feature or request CRTM CRTM labels Aug 7, 2023
@chengdang chengdang self-assigned this Aug 7, 2023
@chengdang
Copy link
Contributor Author

chengdang commented Aug 8, 2023

Attempting to use RT solution as the only input (without atmosphere data structure) for CRTM_RTSolution_WriteFile_netCDF:

  FUNCTION CRTM_RTSolution_WriteFile_netCDF( &
    Sensor_Filename   , &  ! Input
    RTSolution        , &  ! Input
    Quiet             , &  ! Optional input
    Debug             ) &  ! Optional input (Debug output control)
  RESULT( err_stat )

Noticed CRTM_RTSolution_type%n_Layers has never been assigned in CRTM calculations. This variable is needed for layer outputs. Modified Common_RTSolution.F90 to assign CRTM_RTSolution_type%n_Layers in modules Assign_Common_Output*.

Updates Aug 13:
CRTM_RTSolution_type%n_Layers should be assigned when calling CRTM_RTSolution_Create( rts, N_LAYERS )

Once this is added in check_crtm.F90, there is no need to reassign this variable. Restored all changes in module Assign_Common_Output

@chengdang
Copy link
Contributor Author

chengdang commented Aug 10, 2023

Initial code updates have been pushed to branch: https://github.com/JCSDA/CRTMv3/tree/feature/cd_RTSout_netCDF

This branch allows users to output CRTM RT Solution in netCDF files, for example, see

! 8d. **** OUTPUT THE RESULTS TO NetCDF files ****

New CRTM_RTSolution_WriteFile_netCDF module writes RT solution of the same sensor (multi channels) and same profile to one single file. For multiple profiles, it generates multiple files.

@chengdang
Copy link
Contributor Author

chengdang commented Aug 10, 2023

Remaining work:

@chengdang
Copy link
Contributor Author

chengdang commented Aug 14, 2023

Layered outputs have been added in module CRTM_RTSolution_WriteFile_netCDF. For this module to work properly, the users must call CRTM_RTSolution_Create to allocate n_Layers and layered outputs. Check_CRTM.f90 is updated to reflect such requirement: commit 1ee69fd

The following variables are not added to RTS solution for now, since:

  • Wavelength is stored in the spc coeff file.
  • Input angles should be part of the input profiles.

Closing this issue for now.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
CRTM CRTM enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

1 participant