### 
### This is a template to run HSE, assuming that iter1.tmplt has been run.
###


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


VectorRead$QE.FORMAT$( p_abs ) { 
   "" 
}
VectorRead$QE.FORMAT$( t_abs ) { 
   "" 
}
VectorRead$QE.FORMAT$( z_abs ) { 
   "" 
}
MatrixRead$QE.FORMAT$( vmrs ) { 
   "" 
}
h2o_absSet {
}
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


VectorWrite$QE.FORMAT$( z_abs ) { "" }