Skip to content
Snippets Groups Projects
Commit 04c9776e authored by Lucie Bakels's avatar Lucie Bakels
Browse files

updated SPECIES files in test cases

parent c555491c
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
Showing
with 3 additions and 353 deletions
tests/default_options/SPECIES/SPECIES_021 deleted 100644 → 0
+ 0
− 22
View file @ c555491c
&SPECIES_PARAMS
PSPECIES="Xe-133", ! Tracer name
PDECAY=453168.0, ! Species half life
PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas)
PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas)
PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho
PDQUER=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PNDIA=1, ! Dry deposition (particles) - ndia
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=-9.9, ! molweight
POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=-9.9, ! OH Reaction rate - D [K]
POHNCONST=2.0, ! OH Reaction rate - N (no unit)
/
tests/default_options/SPECIES/SPECIES_022 deleted 100644 → 0
+ 0
− 22
View file @ c555491c
&SPECIES_PARAMS
PSPECIES="CO", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas)
PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas)
PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho
PDQUER=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PNDIA=1, ! Dry deposition (particles) - ndia
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=28.0, ! molweight
POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=-9.9, ! OH Reaction rate - D [K]
POHNCONST=2.0, ! OH Reaction rate - N (no unit)
/
tests/default_options/SPECIES/SPECIES_023 deleted 100644 → 0
+ 0
− 22
View file @ c555491c
&SPECIES_PARAMS
PSPECIES="SO2", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=2.0E-05, ! Below-cloud scavenging (gases) - A (weta_gas)
PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas)
PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=-0.9E-9, ! Dry deposition (particles) - rho
PDQUER=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PNDIA=1, ! Dry deposition (particles) - ndia
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=2.0, ! Dry deposition (gases) - D
PHENRY=100000.0, ! Dry deposition (gases) - Henrys const.
PF0=0, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=64.0, ! molweight
POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=-9.9, ! OH Reaction rate - D [K]
POHNCONST=2.0, ! OH Reaction rate - N (no unit)
/
tests/default_options/SPECIES/SPECIES_024
+ 0
− 1
View file @ 04c9776e
......@@ -10,7 +10,6 @@
PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho
PDQUER=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PNDIA=1, ! Dry deposition (particles) - ndia
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
......
tests/default_options/SPECIES/SPECIES_025 deleted 100644 → 0
+ 0
− 22
View file @ c555491c
&SPECIES_PARAMS
PSPECIES="AERO-TRACE", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas)
PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas)
PCRAIN_AERO=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=1400.0, ! Dry deposition (particles) - rho
PDQUER=2.5E-07, ! Dry deposition (particles) - dquer
PDSIGMA=1.25, ! Dry deposition (particles) - dsig
PNDIA=1, ! Dry deposition (particles) - ndia
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=-9.9, ! molweight
POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=-9.9, ! OH Reaction rate - D [K]
POHNCONST=2.0, ! OH Reaction rate - N (no unit)
/
tests/default_options/SPECIES/SPECIES_026 deleted 100644 → 0
+ 0
− 22
View file @ c555491c
&SPECIES_PARAMS
PSPECIES="CH4", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas)
PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas)
PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho
PDQUER=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PNDIA=1, ! Dry deposition (particles) - ndia
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=16.0, ! molweight
POHCCONST=9.65E-20, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=1082.0, ! OH Reaction rate - D [K]
POHNCONST=2.58, ! OH Reaction rate - N (no unit)
/
tests/default_options/SPECIES/SPECIES_027 deleted 100644 → 0
+ 0
− 22
View file @ c555491c
&SPECIES_PARAMS
PSPECIES="C2H6", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas)
PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas)
PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho
PDQUER=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PNDIA=1, ! Dry deposition (particles) - ndia
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=30.0, ! molweight
POHCCONST=15.2E-18, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=498.0, ! OH Reaction rate - D [K]
POHNCONST=2.0, ! OH Reaction rate - N (no unit)
/
tests/default_options/SPECIES/SPECIES_028 deleted 100644 → 0
+ 0
− 7
View file @ c555491c
&SPECIES_PARAMS
PSPECIES="C3H8 ",
PWEIGHTMOLAR= 44.0000000 ,
POHCCONST= 1.55-E17,
POHDCONST= 61.0000000 ,
POHNCONST= 2.00000000 ,
/
tests/default_options/SPECIES/SPECIES_031 deleted 100644 → 0
+ 0
− 22
View file @ c555491c
&SPECIES_PARAMS
PSPECIES="PCB28", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=3.5E-05, ! Below-cloud scavenging (gases) - A (weta_gas)
PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas)
PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho
PDQUER=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PNDIA=1, ! Dry deposition (particles) - ndia
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=4.3, ! Dry deposition (gases) - D
PHENRY=3.1, ! Dry deposition (gases) - Henrys const.
PF0=0.1, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=257.54, ! molweight
POHCCONST=1.07E-11, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=1203.0, ! OH Reaction rate - D [K]
POHNCONST=0.0, ! OH Reaction rate - N (no unit)
/
tests/default_options/SPECIES/SPECIES_034 deleted 100644 → 0
+ 0
− 22
View file @ c555491c
&SPECIES_PARAMS
PSPECIES="G-HCH", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=3.5E-05, ! Below-cloud scavenging (gases) - A (weta_gas)
PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas)
PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho
PDQUER=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PNDIA=1, ! Dry deposition (particles) - ndia
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=4.3, ! Dry deposition (gases) - D
PHENRY=422.0, ! Dry deposition (gases) - Henrys const.
PF0=0.1, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=290.83, ! molweight
POHCCONST=6.21E-11, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=1203.0, ! OH Reaction rate - D [K]
POHNCONST=2.00, ! OH Reaction rate - N (no unit)
/
tests/default_options/SPECIES/SPECIES_040
+ 3
− 4
View file @ 04c9776e
......@@ -10,15 +10,14 @@
PDENSITY=1000.0, ! Dry deposition (particles) - rho
PDQUER=50.0E-06, ! Dry deposition (particles) - dquer (equivalent diameter for shape)
PDSIGMA=3.3, ! Dry deposition (particles) - dsig
PNDIA=1, ! Dry deposition (particles) - ndia
PDRYVEL=-9.9, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=-9.9, ! molweight
POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=-9.9, ! OH Reaction rate - D [K]
POHNCONST=2.0, ! OH Reaction rate - C [cm^3/molecule/sec]
PCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec]
PDCONST=-9.9, ! OH Reaction rate - D [K]
PNCONST=2.0, ! OH Reaction rate - C [cm^3/molecule/sec]
PSHAPE=1, ! 0 for sphere, 1 any shape (defined by axes PLA,PIA,PSA), 2-cylinder, 3-cube, 4-tetrahedron, 5-octahedron, 6-ellipsoid
PASPECTRATIO=20.0 ! Aspect ratio of cylinders: works for PSHAPE=2 only
PLA=940, ! Longest axis in micrometer (Bagheri & Bonadonna 2016): only for PSHAPE=1
......
tests/default_options/SPECIES/SPECIES_050 deleted 100755 → 0
+ 0
− 25
View file @ c555491c
&SPECIES_PARAMS
PSPECIES="DUST-1um ",
PDECAY= -999.900024 ,
PWETA_GAS= -9.89999993E-09,
PWETB_GAS= 0.00000000 ,
PCRAIN_AERO= 1.00000000 ,
PCSNOW_AERO= 1.00000000 ,
PCCN_AERO= 0.150000000 ,
PIN_AERO= 0.020000000 ,
PRELDIFF= -9.89999962 ,
PHENRY= 9.99999972E-10,
PF0= 0.00000000 ,
PDENSITY= 2650.00000 ,
PDQUER= 9.99999997E-07,
PDSIGMA= 1.10000002 ,
PDRYVEL= -9.98999977 ,
PWEIGHTMOLAR= -999.900024 ,
POHCCONST= -9.98999977 ,
POHDCONST= -9.89999993E-09,
POHNCONST= 2.00000000 ,
PAREA_DOW= 7*1.00000000 ,
PAREA_HOUR= 24*1.00000000 ,
PPOINT_DOW= 7*1.00000000 ,
PPOINT_HOUR= 24*1.00000000 ,
/
tests/default_options/SPECIES/SPECIES_051 deleted 100755 → 0
+ 0
− 25
View file @ c555491c
&SPECIES_PARAMS
PSPECIES="DUST-5um ",
PDECAY= -999.900024 ,
PWETA_GAS= -9.89999993E-09,
PWETB_GAS= 0.00000000 ,
PCRAIN_AERO= 1.00000000 ,
PCSNOW_AERO= 1.00000000 ,
PCCN_AERO= 0.150000000 ,
PIN_AERO= 0.020000000 ,
PRELDIFF= -9.89999962 ,
PHENRY= 9.99999972E-10,
PF0= 0.00000000 ,
PDENSITY= 2650.00000 ,
PDQUER= 4.99999987E-06,
PDSIGMA= 1.10000002 ,
PDRYVEL= -9.98999977 ,
PWEIGHTMOLAR= -999.900024 ,
POHCCONST= -9.98999977 ,
POHDCONST= -9.89999993E-09,
POHNCONST= 2.00000000 ,
PAREA_DOW= 7*1.00000000 ,
PAREA_HOUR= 24*1.00000000 ,
PPOINT_DOW= 7*1.00000000 ,
PPOINT_HOUR= 24*1.00000000 ,
/
tests/default_options/SPECIES/specoverview.f90 deleted 100644 → 0
+ 0
− 115
View file @ c555491c
!**********************************************************************
! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010 *
! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa, *
! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann *
! *
! This file is part of FLEXPART. *
! *
! FLEXPART is free software: you can redistribute it and/or modify *
! it under the terms of the GNU General Public License as published by*
! the Free Software Foundation, either version 3 of the License, or *
! (at your option) any later version. *
! *
! FLEXPART is distributed in the hope that it will be useful, *
! but WITHOUT ANY WARRANTY; without even the implied warranty of *
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
! GNU General Public License for more details. *
! *
! You should have received a copy of the GNU General Public License *
! along with FLEXPART. If not, see <http://www.gnu.org/licenses/>. *
!**********************************************************************
implicit none
character(len=11) :: speciesfn
character(len=3) :: aspec
character(len=16) :: pspecies
real :: pdecay, pweta_gas, pwetb_gas, preldiff, phenry, pf0, pdensity, pdquer
real :: pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst
real :: pcrain_aero, pcsnow_aero, pccn_aero, pin_aero
integer :: readerror, unitspecies, specnumber
! declare namelist
namelist /species_params/ &
pspecies, pdecay, pweta_gas, pwetb_gas, &
pcrain_aero, pcsnow_aero, pccn_aero, pin_aero, &
preldiff, phenry, pf0, pdensity, pdquer, &
pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst
unitspecies=4
write(*,*) ' Species | | WetDep(gas) | DryDep(gas) |WetDep(below-C)| WetDep(in-C)|'//&
' DryDepo(particles) Altern| Radioact. | OH Reaction |'
write(*,*) ' Name |molwght| A B | D H f0 | Crain Csnow | ccn in |' //&
' rho dquer dsig vd | Halflife[s]| C** D[T] N*** |'
write(*,*) '--------------|-------|-----------------|------------------|---------------|-------------|'//&
'-------------------------------|------------|----------------------|'
! write(*,*) ' Specie | Radioact. | WetDep(gas) |WetDep(below-C)| WetDep(in-C)| DryDepo(gas) |'//&
! ' DryDepo(particles) Altern| | OH Reaction |'
! write(*,*) ' Name | Halflife[s]| A B | Crain Csnow | ccn in | D H f0 |' //&
! ' rho dquer dsig vd |molwght| C** D[T] N*** |'
! write(*,*) '--------------|------------|-----------------|---------------|-------------|-------------------|'//&
! '-------------------------------|-------|----------------------|'
do specnumber=1,100
write (aspec,'(i0.3)') specnumber
speciesfn='SPECIES_'//aspec
! write(*,*) 'Processing: ',speciesfn
pspecies="" ! read failure indicator value
pdecay=-9.9
pweta_gas=-0.9E-09
pwetb_gas=0.0
pcrain_aero=-9.9
pcsnow_aero=-9.9
pccn_aero=-9.9
pin_aero=-9.9
preldiff=-9.9
phenry=0.0
pf0=0.0
pdensity=-0.9E09
pdquer=0.0
pdsigma=0.0
pdryvel=-9.99
pohcconst=-9.9
pohdconst=-9.9
pohnconst=2.0
pweightmolar=-9.9
! Open the SPECIES file and read species names and properties
!************************************************************
open(unitspecies,file=speciesfn,status='old',form='formatted',err=998)
read(unitspecies,species_params,err=998)
close(unitspecies)
write(*,45) specnumber,' ',pspecies,'|',pweightmolar,'|',pweta_gas,' ',pwetb_gas,'|', &
preldiff,' ',phenry,' ',pf0,'|', &
pcrain_aero,' ',pcsnow_aero,'|',pccn_aero,' ',pin_aero,'|', &
pdensity,pdquer,pdsigma,pdryvel,'|',pdecay,'|',pohcconst,pohdconst,pohnconst,'|'
45 format(i3,a1,a11,a1,f7.1,a1,e8.1,a1,f8.2,a1, &
f4.1,a1,e8.1,a1,f4.1,a1, &
f7.1,a1,f7.1,a1,f6.1,a1,f6.1,a1, &
e8.1,e9.1,f7.1,f7.2,a1,f12.1,a1,e8.1,f7.1,f7.1,a1)
! write(*,45) specnumber,' ',pspecies,'|',pdecay,'|',pweta_gas,' ',pwetb_gas,'|',pcrain_aero,' ', &
! pcsnow_aero,'|',pccn_aero,' ',pin_aero,'|',preldiff,' ',phenry,' ',pf0,'|', &
! pdensity,pdquer,pdsigma,pdryvel,'|',pweightmolar,'|',pohcconst,pohdconst,pohnconst,'|'
!5 format(i3,a1,a11,a1,f12.1,a1,e8.1,a1,f8.2,a1,f7.1,a1,f7.1,a1,f6.1,a1,f6.1,a1,f5.1,a1,e8.1,a1,f4.1,a1, &
! e8.1,e9.1,f7.1,f7.2,a1,f7.1,a1,e8.1,f7.1,f7.1,a1)
998 continue
enddo
write(*,*) '** unit [cm^3/molec/s] (in FLEXPART version 9.2 and below this had unit [cm3/s], note the unit is now changed!)'
write(*,*) '*** no unit'
print*,'rho: density'
end
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment