Changeset 6481010 in flexpart.git

Timestamp:

Jan 26, 2017, 3:10:55 PM (7 years ago)

Author:

Sabine <sabine.eckhardt@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie

Children:

7725d96

Parents:

56a2e54

Message:

vertransform includes option to fix the height of the inputfields

File:

• src/verttransform.f90 (modified) (3 diffs)

src/verttransform.f90

@@ -91 +91 @@
 
   logical :: init = .true.
+  logical :: init_w = .false.
+  logical :: init_r = .false.
+
 
   !ZHG SEP 2014 tests  
@@ -123 +126 @@
   if (init) then
 
+
+    if (init_r) then
+
+        open(333,file='heights.txt', &
+          form='formatted')
+        do kz=1,nuvz
+            read(333,*) height(kz)
+        end do
+        close(333)
+        write(*,*) 'height read'
+    else
+
+
 ! Search for a point with high surface pressure (i.e. not above significant topography)
 ! Then, use this point to construct a reference z profile, to be used at all times
@@ -161 +177 @@
     end do
 
+    if (init_w) then
+        open(333,file='heights.txt', &
+          form='formatted')
+        do kz=1,nuvz
+              write(333,*) height(kz)
+        end do
+        close(333)
+    endif
+
+    endif ! init
 
 ! Determine highest levels that can be within PBL