#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Mon Dec 11 17:36:45 2023

@author: mzannese
"""

import matplotlib.pyplot as plt
import numpy as np

fig, ax = plt.subplots(2, figsize=(30, 10), sharex=True)
ax[0].set_ylim(-60, 300)
ax[1].set_ylim(-60,300)

ax[0].set_xlim(2.84,3.25)

data = np.loadtxt('/PDRs4All_d203-506_fig3_Zannese24.txt')

wave = data[:, 0]
intensity=data[:,3]
modelOH=data[:,4]

ax[0].step(wave, intensity, c='black',where='mid')
ax[0].fill_between(wave,modelOH,step='mid',lw=1.5,color='r',zorder=0,alpha=0.8)
ax[1].step(wave,modelOH,c='r',where='mid')


fig.supylabel('Intensity (MJy/sr)')
fig.supxlabel('Wavelength ($\mu$m)')