Notice: wiki.ucar.edu has been updated. Please see the announcement for details.

Page tree

Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

  • Count data of different surface types

  • Set relative weight value

  • Interpolate model fields to observation location, call crtm and create jacobians. Output both tsim and tsim_clr for allsky

radiance operator

 

1
2
3
4
5
6
7
8
9
10
11
12
13
if (lcw4crtm) then
   call call_crtm(obstype,dtime,data_s(:,n),nchanl,nreal,ich, &
        tvp,qvp,clw_guess,prsltmp,prsitmp, &
        trop5,tzbgr,dtsavg,sfc_speed, &
        tsim,emissivity,ptau5,ts,emissivity_k, &
        temp,wmix,jacobian,error_status,tsim_clr=tsim_clr)
else
   call call_crtm(obstype,dtime,data_s(:,n),nchanl,nreal,ich, &
        tvp,qvp,clw_guess,prsltmp,prsitmp, &
        trop5,tzbgr,dtsavg,sfc_speed, &
        tsim,emissivity,ptau5,ts,emissivity_k, &
        temp,wmix,jacobian,error_status)
endif

 

  • If the CRTM returns an error flag, do not assimilate any channels for this ob and set the QC flag to ifail_crtm_qc.

...