Clear-sky masking in COSP-RTTOV

driver/run/cosp2_swath_rttov_input_nl.um_global.txt

@@ -0,0 +1,109 @@
+! (c) British Crown Copyright 2008, the Met Office.
+! All rights reserved.
+!
+! Redistribution and use in source and binary forms, with or without modification, are permitted
+! provided that the following conditions are met:
+!
+!     * Redistributions of source code must retain the above copyright notice, this list
+!       of conditions and the following disclaimer.
+!     * Redistributions in binary form must reproduce the above copyright notice, this list
+!       of conditions and the following disclaimer in the documentation and/or other materials
+!       provided with the distribution.
+!     * Neither the name of the Met Office nor the names of its contributors may be used
+!       to endorse or promote products derived from this software without specific prior written
+!       permission.
+!
+! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR
+! IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND
+! FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR
+! CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+! DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+! DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER
+! IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
+! OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+
+! Namelist that sets up the main COSP options
+&COSP_INPUT
+  NPOINTS=1728, ! 1728,6912
+  NPOINTS_IT=1000,! Max number of gridpoints to be processed in one iteration
+  NCOLUMNS=20,  ! Number of subcolumns
+  NLEVELS=54,   ! Number of model levels
+  USE_VGRID=.true., ! Use fixed vertical grid for outputs? (if .true. then you need to define number of levels with Nlr)
+  NLVGRID=40,       ! Number of levels in statistical outputs (only used if USE_VGRID=.true.)
+  CSAT_VGRID=.true., ! CloudSat vertical grid? (if .true. then the CloudSat standard grid is used for the outputs.
+                     !  USE_VGRID needs also be .true.)
+  DINPUT='./',  ! Directory where the input files are located. Useful when processing multiple files.
+                ! Leave blank ('') if you are using the full path in FINPUT.
+  FINPUT='../data/inputs/UKMO/cosp_input.um_global.nc', ! List input NetCDF files
+  FOUTPUT='../data/outputs/UKMO/cosp2_swath_output.um_global.nc',
+  !----------------------------------------------------------------------------------
+  !--------------- Inputs related to radar simulations
+  !----------------------------------------------------------------------------------
+  cloudsat_RADAR_FREQ=94.0, ! CloudSat radar frequency (GHz)
+  SURFACE_RADAR=0, ! surface=1, spaceborne=0
+  cloudsat_use_gas_abs=1,   ! include gaseous absorption? yes=1,no=0
+  cloudsat_do_ray=0,        ! calculate/output Rayleigh refl=1, not=0
+  cloudsat_k2=-1,           ! |K|^2, -1=use frequency dependent default
+  use_precipitation_fluxes=.true.,  ! True if precipitation fluxes are input to the algorithm
+  cloudsat_micro_scheme='MMF_v3_single_moment', !'MMF_v3.5_two_moment'
+  !----------------------------------------------------------------------------------
+  !---------------- Inputs related to lidar simulations
+  !----------------------------------------------------------------------------------
+  lidar_ice_type=0,    ! Ice particle shape in lidar calculations (0=ice-spheres ; 1=ice-non-spherical)
+  OVERLAP=3,           !  overlap assumption used by scops: 1=max, 2=rand, 3=max/rand
+  !----------------------------------------------------------------------------------
+  !---------------- Inputs related to ISCCP simulator
+  !----------------------------------------------------------------------------------
+  ISCCP_TOPHEIGHT=1,  !  1 = adjust top height using both a computed
+                       !  infrared brightness temperature and the visible
+                       !  optical depth to adjust cloud top pressure. Note
+                       !  that this calculation is most appropriate to compare
+                       !  to ISCCP data during sunlit hours.
+                      !  2 = do not adjust top height, that is cloud top
+                       !  pressure is the actual cloud top pressure
+                       !  in the model
+                      !  3 = adjust top height using only the computed
+                       !  infrared brightness temperature. Note that this
+                       !  calculation is most appropriate to compare to ISCCP
+                       !  IR only algortihm (i.e. you can compare to nighttime
+                       !  ISCCP data with this option)
+  ISCCP_TOPHEIGHT_DIRECTION=2,   ! direction for finding atmosphere pressure level
+                                 ! with interpolated temperature equal to the radiance
+                                 ! determined cloud-top temperature
+                                 ! 1 = find the *lowest* altitude (highest pressure) level
+                                 ! with interpolated temperature equal to the radiance
+                                 ! determined cloud-top temperature
+                                 ! 2 = find the *highest* altitude (lowest pressure) level
+                                 ! with interpolated temperature equal to the radiance
+                                 ! determined cloud-top temperature. This is the
+                                 ! default value since V4.0 of the ISCCP simulator.
+                                 ! ONLY APPLICABLE IF top_height EQUALS 1 or 3
+  !----------------------------------------------------------------------------------
+  !-------------- Swathing inputs
+  !----------------------------------------------------------------------------------
+  N_SWATHS_ISCCP=2,
+  SWATH_LOCALTIMES_ISCCP=5.25,18.0,
+  SWATH_WIDTHS_ISCCP=2000,2000,
+  N_SWATHS_MISR=2,
+  SWATH_LOCALTIMES_MISR=6.5,18.25,
+  SWATH_WIDTHS_MISR=2000,2000,
+  N_SWATHS_MODIS=2,
+  SWATH_LOCALTIMES_MODIS=5.75,13.5,
+  SWATH_WIDTHS_MODIS=2000,2000,
+  N_SWATHS_CSCAL=2,
+  SWATH_LOCALTIMES_CSCAL=6.75,19.25,
+  SWATH_WIDTHS_CSCAL=2000,2000,
+  N_SWATHS_PARASOL=2,
+  SWATH_LOCALTIMES_PARASOL=6.75,19,
+  SWATH_WIDTHS_PARASOL=2000,2000,
+  N_SWATHS_ATLID=2,
+  SWATH_LOCALTIMES_ATLID=7.25,19.5,
+  SWATH_WIDTHS_ATLID=2000,2000,
+  !----------------------------------------------------------------------------------
+  !-------------- RTTOV inputs
+  !----------------------------------------------------------------------------------
+  rttov_Ninstruments=1,
+  rttov_instrument_namelists='instrument_nls/cosp2_rttov_inst1.txt',
+  rttov_verbose=.false.
+/

driver/run/instrument_nls/cosp2_rttov_inst1.txt

@@ -0,0 +1,85 @@
+! (c) British Crown Copyright 2008, the Met Office.
+! All rights reserved.
+!
+! Redistribution and use in source and binary forms, with or without modification, are permitted
+! provided that the following conditions are met:
+!
+!     * Redistributions of source code must retain the above copyright notice, this list
+!       of conditions and the following disclaimer.
+!     * Redistributions in binary form must reproduce the above copyright notice, this list
+!       of conditions and the following disclaimer in the documentation and/or other materials
+!       provided with the distribution.
+!     * Neither the name of the Met Office nor the names of its contributors may be used
+!       to endorse or promote products derived from this software without specific prior written
+!       permission.
+!
+! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR
+! IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND
+! FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR
+! CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+! DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+! DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER
+! IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
+! OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+
+! Namelist that sets up the main COSP options
+&RTTOV_INPUT
+  !-------------- Logical for constructing outputs and determining settings
+  Lrttov_bt=.true.,                 ! Calculate RTTOV brightness temps
+  Lrttov_rad=.true.,                ! Calculate RTTOV radiances
+  Lrttov_refl=.false.,              ! Calculate RTTOV reflectivities
+  Lrttov_cld=.true.,                ! Calculate RTTOV all-sky + cloudy-sky values
+  Lrttov_maskclear = .true.,        ! Mask RTTOV clear-sky values using model cloud fields.
+  Lrttov_aer=.false.,                ! Calculate RTTOV all-sky + cloudy-sky values
+  Lrttov_cldparam=.false.,          ! Use user input cloud optical parameters
+  Lrttov_aerparam=.false.,          ! Use user input aerosol optical parameters
+  Lrttov_gridbox_cldmmr=.true.,     ! Do cloud mixing ratios represent gridcell means?
+  Ldo_nlte_correction=.false.,       ! Correct for non-local thermal equilibrium effects? Requires appropriate coefficients.
+  Lrttov_pc=.false.,                ! Run with PC-RTTOV
+  !-------------- Number of channels to reconstruct. Should match the length of entries in "channel_filepath"
+  nchannels_rec=5, ! 2645 for all AIRS L1C channels, 45 channel subset mostly from Sergio
+  !-------------- RTTOV Filepaths 
+  Lchannel_filepath=.true.,
+  channel_filepath='rttov_channel_input_AIRSL1C_subset.csv'
+  rttov_srcDir='/glade/u/home/jonahshaw/w/RTTOV/',
+  rttov_coefDir='rtcoef_rttov13/',
+  OD_coef_filepath='AIRS_L1C/rtcoef_eos_2_airs_l1c_7gas.H5', ! Full version with N20 and CH4
+  cld_coef_filepath='AIRS_L1C/sccldcoef_eos_2_airs_l1c.H5',
+  !-------------- RTTOV Trace Gas Logicals (must match coefficient files)
+  SO2_data=.false.,
+  N2O_data=.true., ! .true. for 7gas
+  CO_data=.false.,
+  CO2_data=.true.,
+  CH4_data=.true., ! .true. for 7gas
+  ozone_data=.true.,
+  clw_data=.false.,        ! Only used for MW calculations
+  ! -------------- RTTOV Trace Gas Concentrations
+  ! -------------- Specify units in RTTOV options set-up
+  ! ---------------Standard units are kg/kg over moist air
+  rttov_gas_units=3, ! 0: ppmv over dry air, 1: kg/kg over moist air, 2: ppmv over moist air, 3: kg/kg over dry air
+  user_tracegas_input=.false., ! Have RTTOV use the trace gas concentrations below
+  CO2_mr=6.09e-04, 
+  CH4_mr=9.139e-07,
+  CO_mr=2.098e-07,
+  N2O_mr=4.665e-07,
+  SO2_mr=2.0e-11,  ! RTTOV hard lower limit is 1.0e-11 ppmv
+  !-------------- PC-RTTOV ------------------
+  ipcbnd=1,       ! Always 1 for RTTOV v13
+  ipcreg=2,       ! Set the number of predictors (channel outputs)
+  npcscores=100,  ! Determine how many principal components are used (greater accuracy, greater cost)
+  !-------------- RTTOV Other --------------- 
+  ! Control degree of parallelization in RTTOV 1=no parallel, >1 more parallel
+  rttov_nthreads=1,
+  rttov_ZenAng=0.0,               ! Satellite Zenith Angle
+  ! Set local time sampling with swath widths
+  rttov_Nlocaltime=0,            	! Number of local times to sample at
+  rttov_localtime=1,0,     	! Local times to sample at in hrs [0,24]
+  rttov_localtime_width=1550,0,  ! Swath width (km)
+  ! ---------------------- How to parameterize the optical properties of cloud ice and cloud liquid (see RTTOV user manual pg 75)
+  rttov_clw_scheme=2,   ! Should be 2 if model supplies effective radius and not cloud type.
+  rttov_ice_scheme=1,   ! 1 uses effective radii values if rttov_icede_param=0, 2 and 3 ignore icede entirely.
+  rttov_icede_param=4,  ! 0 if icede should be used from model, if ice_scheme = 1 then 1-4 are different options for computing icede separately.
+  ! ---------------------- How to handle atmosphere between the top of supplied values and the RTTOV instrument.  
+  rttov_extendatmos=1,    ! 0: Atmosphere ends at the top of supplied values. 1: CAM-like additional interfaces at 1e-4 Pa and 2.25 Pa.
+/

driver/run/rttov_channel_input_AIRSL1C_subset.csv

@@ -0,0 +1,45 @@
+  41 ,1 ,0
+  54 ,1 ,0
+  72 ,1 ,0
+ 202 ,1 ,0
+ 234 ,1 ,0
+ 294 ,1 ,0
+ 336 ,1 ,0
+ 378 ,1 ,0
+ 485 ,1 ,0
+ 489 ,1 ,0
+ 572 ,1 ,0
+ 794 ,1 ,0
+ 960 ,1 ,0
+ 961 ,1 ,0
+ 967 ,1 ,0
+1055 ,1 ,0
+1075 ,1 ,0
+1089 ,1 ,0
+1110 ,1 ,0
+1130 ,1 ,0
+1158 ,1 ,0
+1325 ,1 ,0
+1511 ,1 ,0
+1520 ,1 ,0
+1697 ,1 ,0
+1723 ,1 ,0
+1805 ,1 ,0
+1852 ,1 ,0
+1862 ,1 ,0
+1866 ,1 ,0
+1937 ,1 ,0
+2088 ,1 ,0
+2138 ,1 ,0
+2164 ,1 ,0
+2165 ,1 ,0
+2166 ,1 ,0
+2176 ,1 ,0
+2186 ,1 ,0
+2383 ,1 ,0
+2411 ,1 ,0
+2588 ,1 ,0
+2592 ,1 ,0
+2600 ,1 ,0
+2606 ,1 ,0
+2620 ,1 ,0