### 
### This is a template to calculate the absorption to a Q. 
###    
### The following fields of QE must exist:
###    
###    QE.ABS_PARTS
###    QE.FORMAT


### Spectroscopy: ############################################################
IF isempty(Q.OTHER_TAGS)
tgsDefine { 
   [$Q.RETRIEVAL_TAGS$] 
}
END
IF isempty(Q.RETRIEVAL_TAGS)
tgsDefine { 
   [$Q.OTHER_TAGS$] 
}
END
IF ~isempty(Q.OTHER_TAGS) & ~isempty(Q.RETRIEVAL_TAGS)
tgsDefine { 
   [$Q.RETRIEVAL_TAGS$,$Q.OTHER_TAGS$] 
}
END
lineshapeDefine { 
   shape               = "$Q.LINESHAPE$"
   normalizationfactor = "$Q.LINESHAPE_FACTOR$" 
   cutoff              =  $Q.LINESHAPE_CUTOFF$
} 
linesReadFrom$Q.LINEFORMAT$ { 
   filename = "$Q.SPECTRO_DIR$$Q.FILESEP$$Q.LINEFILE$"
   fmin     = 0
   fmax     = 10e18 
		}
lines_per_tgCreateFromLines {
}  
#
@qpi_continua



### p_abs ####################################################################
VectorReadAscii( p_abs ) { 
   "$Q.CALCGRIDS_DIR$$Q.FILESEP$$Q.P_ABS$" 
}


### Monochromatic frequency grid #############################################
VectorReadAscii( f_mono ) { 
   "$Q.CALCGRIDS_DIR$$Q.FILESEP$$Q.F_MONO$" 
}


### The atmosphere ###########################################################
MatrixReadAscii ( raw_ptz ) { 
   filename = "$Q.APRIORI_PTZ$" 
}
raw_vmrsReadFromScenario {
   basename = "$Q.APRIORI_VMR$" 
}
AtmFromRaw {}
h2o_absSet {
}
n2_absSet {
}


### The geoid ###############################################################
r_geoidStd {
}


### Hydrostatic eq. ###########################################################
IF Q.HSE_IN_ON
hseSet{
   pref  = $Q.HSE_PREF$
   zref  = $Q.HSE_ZREF$
   g0    = 9.81
   niter = 2 }
hseCalc{}
ELSE
hseOff{}
END


### Absorption ################################################################
absCalc {
}
VectorWrite$QE.FORMAT$( f_mono ) { "" }
VectorWrite$QE.FORMAT$( p_abs ) { "" }
VectorWrite$QE.FORMAT$( t_abs ) { "" }
VectorWrite$QE.FORMAT$( z_abs ) { "" }
MatrixWrite$QE.FORMAT$( abs ) { "" }
ArrayOfMatrixWrite$QE.FORMAT$( abs_per_tg ) { "" }
IF QE.ABS_PARTS
ArrayOfMatrixWrite$QE.FORMAT$( abs_per_tg ) { "" }
END