Skip to content

TL implementation for aircraft/user-defined level upward/downward #189

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

Open
chengdang wants to merge 2 commits into
base: develop
Choose a base branch
from
Open

TL implementation for aircraft/user-defined level upward/downward #189

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
1a8cd44
TL implementation for aircraft/user-defined level upward/downward
chengdang Sep 25, 2024
bcb6371
Merge branch 'develop' into feature/cd_AirCraft_TL
chengdang Dec 9, 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
42 changes: 38 additions & 4 deletions src/CRTM_Tangent_Linear_Module.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -657,9 +657,9 @@ FUNCTION profile_solution (m, Opt, AncillaryInput) RESULT( Error_Status )
CALL Display_Message( ROUTINE_NAME, Message, Error_Status )
RETURN
END IF
Atm_TL%Height = Atm%Height

Atm_TL%Height = Atm%Height

! ...Check the total number of Atm layers
IF ( Atm%n_Layers > MAX_N_LAYERS .OR. Atm_TL%n_Layers > MAX_N_LAYERS) THEN
Error_Status = FAILURE
Expand Down Expand Up @@ -859,6 +859,40 @@ FUNCTION profile_solution (m, Opt, AncillaryInput) RESULT( Error_Status )
Predictor_TL(nt) , & ! Output
PVar(nt) ) ! Internal variable input

! Process aircraft pressure altitude
IF ( Opt%Aircraft_Pressure > ZERO ) THEN
RTV(nt)%aircraft%rt = .TRUE.
RTV(nt)%aircraft%idx = CRTM_Get_PressureLevelIdx(Atm, Opt%Aircraft_Pressure)
! ...Issue warning if profile level is TOO different from flight level
IF ( ABS(Atm%Level_Pressure(RTV(nt)%aircraft%idx)-Opt%Aircraft_Pressure) > AIRCRAFT_PRESSURE_THRESHOLD ) THEN
WRITE( Message,'("Difference between aircraft pressure level (",es22.15,&
&"hPa) and closest input profile level (",es22.15,&
&"hPa) is larger than recommended (",f4.1,"hPa) for profile #",i0)') &
Opt%Aircraft_Pressure, Atm%Level_Pressure(RTV(nt)%aircraft%idx), &
AIRCRAFT_PRESSURE_THRESHOLD, m
CALL Display_Message( ROUTINE_NAME, Message, WARNING )
END IF
ELSE
RTV(nt)%aircraft%rt = .FALSE.
END IF

! Process observing downward radiance, Obs_4_downward_P = ZERO means at surface
! Obs_4_downward_P > ZERO, sensor at the pressure
IF ( Opt%Obs_4_downward_P > ZERO ) THEN
RTV(nt)%Obs_4_downward%rt = .TRUE.
RTV(nt)%Obs_4_downward%idx = CRTM_Get_PressureLevelIdx(Atm, Opt%Obs_4_downward_P)
! ...Issue warning if profile level is TOO different from flight level
IF ( ABS(Atm%Level_Pressure(RTV(nt)%Obs_4_downward%idx)-Opt%Obs_4_downward_P) > AIRCRAFT_PRESSURE_THRESHOLD ) THEN
WRITE( Message,'("Difference between Obs pressure level (",es22.15,&
&"hPa) and closest input profile level (",es22.15,&
&"hPa) is larger than recommended (",f4.1,"hPa) for profile #",i0)') &
Opt%Obs_4_downward_P, Atm%Level_Pressure(RTV%Obs_4_downward%idx), &
AIRCRAFT_PRESSURE_THRESHOLD, m
CALL Display_Message( ROUTINE_NAME, Message, WARNING )
END IF
ELSE
RTV(nt)%Obs_4_downward%rt = .FALSE.
END IF

! Compute predictors for AtmAbsorption calcs
! ...Allocate the predictor structure
Expand Down Expand Up @@ -1339,7 +1373,7 @@ FUNCTION profile_solution (m, Opt, AncillaryInput) RESULT( Error_Status )
SensorIndex , & ! Input
ChannelIndex , & ! Input
RTSolution(ln,m)) ! Input/Output

CALL CRTM_Compute_Reflectivity_TL(Atm , & ! Input
AtmOptics(nt) , & ! Input
AtmOptics_TL(nt) , & ! Input
Expand Down
Loading