From 8a522e05c3eabbcbf6172c2a6fa59a32ca2d6651 Mon Sep 17 00:00:00 2001
From: Sabine <sabine.eckhardt@nilu.no>
Date: Thu, 27 Jun 2024 11:56:24 +0200
Subject: [PATCH] added specias bca - black carbon aerosol - concentration
instead of mr calculated
---
source/simulate.f90 | 18 +++++++++++++++---
1 file changed, 15 insertions(+), 3 deletions(-)
diff --git a/source/simulate.f90 b/source/simulate.f90
index be349b7..6f7759c 100644
--- a/source/simulate.f90
+++ b/source/simulate.f90
@@ -151,7 +151,7 @@ subroutine simulate(iter, files, config, fluxes, obs, states, grad_o, cost_o)
else
lsatellite = .false.
endif
- print*, 'simulate: lsatellite = ',lsatellite
+! print*, 'simulate: lsatellite = ',lsatellite
! read correction factor for dry dir (ratio rho_wet/rho_dry)
if ( month.ne.prevmonth ) then
@@ -238,7 +238,13 @@ subroutine simulate(iter, files, config, fluxes, obs, states, grad_o, cost_o)
! convert s.m3/kg to s.m2/kg
grid = grid/outheight(1)
! convert from equivalent ppt to observation units (e.g. ppmv)
- grid = grid*config%coeff*mmair/config%molarmass
+ ! sec not do it for bca
+ if ( config%spec.eq.'bca' ) then
+ grid = grid*config%coeff
+ else
+ grid = grid*config%coeff*mmair/config%molarmass
+ endif
+
! apply numerical scaling
grid = grid/numscale
! if ( lsatellite ) then
@@ -306,7 +312,13 @@ subroutine simulate(iter, files, config, fluxes, obs, states, grad_o, cost_o)
! convert s.m3/kg to s.m2/kg
gridnest = gridnest/outheight(1)
! convert from equivalent ppt to observation units (e.g. ppmv)
- gridnest = gridnest*config%coeff*mmair/config%molarmass
+ ! sec not do it for bca
+ if ( config%spec.eq.'bca' ) then
+ gridnest = gridnest*config%coeff
+ else
+ gridnest = gridnest*config%coeff*mmair/config%molarmass
+ endif
+
! apply numerical scaling
gridnest = gridnest/numscale
! if ( lsatellite ) then
--
GitLab