Improved deposition parametrisation

.gitignore

@@ -81,35 +81,17 @@ Snappy.egg-info
 /src/naccident/snap_batch_copy.f90
 /src/naccident/snap_batch_copy.F90
 
-/src/nbomb/*.inc
-/src/nbomb/*.mod
-/src/nbomb/bsnap_nbomb
-/src/nbomb/create_bomb_input
-/src/nbomb/snap_batch_copy.f90
-/src/nbomb/snap_batch_copy.F90
-
 /src/test/testDateCalc
 /src/test/find_parameters_fi_test
 /src/test/snap_batch_copy.F90
 
-/src/traj/*.inc
-/src/traj/*.mod
-/src/traj/bsnap_traj
-/src/traj/create_traj_input
-/src/traj/snap_batch_copy.f90
-/src/traj/snap_batch_copy.F90
-
-/src/volcano/*.mod
-/src/volcano/*.inc
-/src/volcano/snapargos.inc
-/src/volcano/bsnap_volcano
-/src/volcano/snap_batch_copy.f90
-/src/volcano/snap_batch_copy.F90
 
 src/test/snap_testdata
 
-
 docs/html
 docs/latex
 
 bsnap_*
+
+utils/SnapPy/Snappy.egg-info
+src/test/snap_testdata

Dockerfile

@@ -0,0 +1,37 @@
+FROM ubuntu:22.04 AS builder
+
+ENV SNAP_FIMEX_VERSION=1.9
+
+RUN apt-get update -y && \
+    apt-get install -y software-properties-common && \
+    add-apt-repository -y ppa:met-norway/fimex && \
+    apt-get update -y && \
+    apt-get install -y libnetcdff-dev gfortran libfimex-$SNAP_FIMEX_VERSION-dev make && \
+    apt-get remove -y software-properties-common && \
+    apt-get -y autoremove && \
+    rm -rf /var/lib/apt/lists/*
+
+WORKDIR /snap
+ARG VERSION="latest"
+ADD src .
+RUN ln --symbolic --force gcc_pkgconfig.mk current.mk
+ENV SNAP_USE_OMP=1
+RUN make clean && make
+
+FROM ubuntu:22.04
+
+ENV SNAP_FIMEX_VERSION=1.9
+
+RUN apt-get update -y && \
+    apt-get install -y software-properties-common && \
+    add-apt-repository -y ppa:met-norway/fimex && \
+    apt-get update -y && \
+    apt-get install -y libnetcdff7 libfimex-$SNAP_FIMEX_VERSION-0 tini && \
+    apt-get remove -y software-properties-common && \
+    apt-get -y autoremove && \
+    rm -rf /var/lib/apt/lists/*
+
+WORKDIR /snap
+COPY --from=builder /snap/naccident/bsnap_naccident /snap/bsnap
+
+ENTRYPOINT ["/usr/bin/env", "TINI_VERBOSITY=0", "/usr/bin/tini", "-g", "--", "/snap/bsnap"]

src/common/allocateFields.f90

@@ -24,10 +24,14 @@ module allocateFieldsML
       depdry, depwet, accprec, avgprec, avghbl, precip, &
       pmsl1, pmsl2, field1, field2, field3, field4, field3d1, xm, ym, &
       garea, field_hr1, field_hr2, field_hr3, hbl_hr, &
+      precip3d, cw3d, &
       max_column_scratch, max_column_concentration, &
       aircraft_doserate, aircraft_doserate_scratch, t1_abs, t2_abs, &
-      aircraft_doserate_threshold_height, &
-      total_activity_released, total_activity_lost_domain, total_activity_lost_other
+      aircraft_doserate_threshold_height, vd_dep, &
+      xflux, yflux, hflux, t2m, z0, leaf_area_index, &
+      roa, ustar, monin_l, raero, vs, rs, &
+      total_activity_released, total_activity_lost_domain, total_activity_lost_other, &
+      wscav, cloud_cover
   USE snapfilML, only: idata, fdata
   USE snapgrdML, only: ahalf, bhalf, vhalf, alevel, blevel, vlevel, imodlevel, &
       compute_column_max_conc, compute_aircraft_doserate, aircraft_doserate_threshold
@@ -44,6 +48,9 @@ subroutine allocateFields
   USE snapdimML, only: nx, ny, nk, output_resolution_factor, ldata, maxsiz
   USE snapparML, only: ncomp, iparnum
   USE releaseML, only: mplume, iplume, plume_release, mpart
+  USE drydep, only: drydep_scheme, DRYDEP_SCHEME_EMEP, DRYDEP_SCHEME_EMERSON, &
+    DRYDEP_SCHEME_ZHANG
+  USE snapmetML, only: met_params
 
   logical, save :: FirstCall = .TRUE.
   integer :: AllocateStatus
@@ -147,6 +154,14 @@ subroutine allocateFields
 
   ALLOCATE ( precip(nx,ny), STAT = AllocateStatus)
   IF (AllocateStatus /= 0) ERROR STOP errmsg
+  if (met_params%use_3d_precip) then
+    ALLOCATE(precip3d(nx,ny,nk), cw3d(nx,ny,nk), STAT=AllocateStatus)
+    if (AllocateStatus /= 0) ERROR STOP errmsg
+    ALLOCATE(wscav(nx,ny,nk,ncomp),STAT=AllocateStatus)
+    if (AllocateStatus /= 0) ERROR STOP errmsg
+    wscav(:,:,:,:) = 0.0
+    allocate(cloud_cover(nx,ny,nk))
+  endif
 
 ! the calculation-fields
   ALLOCATE ( avghbl(nx,ny), STAT = AllocateStatus)
@@ -222,6 +237,15 @@ subroutine allocateFields
   total_activity_released(:) = 0.0
   total_activity_lost_domain(:) = 0.0
   total_activity_lost_other(:) = 0.0
+  if (drydep_scheme == DRYDEP_SCHEME_EMEP .or. &
+      drydep_scheme == DRYDEP_SCHEME_ZHANG .or. &
+      drydep_scheme == DRYDEP_SCHEME_EMERSON) then
+    allocate(vd_dep(nx,ny,ncomp), STAT=AllocateStatus)
+    if (AllocateStatus /= 0) ERROR STOP errmsg
+    allocate(xflux, yflux, hflux, t2m, z0, leaf_area_index, mold=ps2)
+    allocate(roa(nx, ny))
+    allocate(ustar, monin_l, raero, vs, rs, mold=roa)
+  endif
 
 end subroutine allocateFields
 
@@ -276,6 +300,10 @@ subroutine deAllocateFields
   DEALLOCATE ( pmsl2)
 
   DEALLOCATE ( precip)
+  if (allocated(precip3d)) deallocate(precip3d)
+  if (allocated(cw3d)) deallocate(cw3d)
+  if (allocated(wscav)) deallocate(wscav)
+  if (allocated(cloud_cover)) deallocate(cloud_cover)
 
   DEALLOCATE ( avghbl )
   DEALLOCATE ( avgprec )
@@ -316,6 +344,11 @@ subroutine deAllocateFields
   DEALLOCATE ( plume_release )
 
   DEALLOCATE( total_activity_released, total_activity_lost_domain, total_activity_lost_other )
+  if (allocated(vd_dep)) then
+    deallocate(vd_dep)
+    deallocate(xflux, yflux, hflux, t2m, z0, leaf_area_index)
+    deallocate(roa, ustar, monin_l, raero, vs, rs)
+  endif
 
 end subroutine deAllocateFields
 end module allocateFieldsML