Changeset d7935de in flexpart.git for src/readreceptors.f90

Timestamp:

Sep 11, 2018, 6:06:31 PM (6 years ago)

Author:

pesei <petra seibert at univie ac at>

Branches:

univie

Children:

93786a1

Parents:

34f1452

Message:

modify most input read subroutines

changed some variable names (mostly for I-N reasons)
includes two names appearing also in timemanager, com_mod
corrected a few mistakes
simplified some parts of code
changed options/RELEASES which is in nml fmt correspondingly

File:

• src/readreceptors.f90 (modified) (9 diffs)

src/readreceptors.f90

@@ -27 +27 @@
   !                                                                            *
   !     Author: A. Stohl                                                       *
-  !     1 August 1996                                                          *
-  !     HSO, 14 August 2013
-  !     Added optional namelist input
+  !     1 August 1996        
+  !                                                  *
+  !     HSO, 14 August 2013: Added optional namelist input
+  !     PS, 2/2015: access= -> position=
+  !     PS, 6/2015: variable names, simplify code
   !                                                                            *
   !*****************************************************************************
@@ -52 +54 @@
   character(len=16) :: receptor
 
-  integer :: readerror
-  real :: lon,lat   ! for namelist input, lon/lat are used instead of x,y
+  integer :: ios
+  real :: xlon,ylat   ! for namelist input, lon/lat are used instead of x,y
   integer,parameter :: unitreceptorout=2
 
   ! declare namelist
-  namelist /receptors/ &
-    receptor, lon, lat
+  namelist /nml_receptors/ receptor, xlon, ylat
 
-  lon=-999.9
-
-  ! For backward runs, do not allow receptor output. Thus, set number of receptors to zero
+!CPS I comment this out - why should we not have receptor output in bwd runs? 
+  ! For backward runs, do not allow receptor output. Thus, set number of
+  ! receptors to zero
   !*****************************************************************************
 
-  if (ldirect.lt.0) then
-    numreceptor=0
-    return
-  endif
+!  if (ldirect.lt.0) then
+!    numreceptor=0
+!    return
+!  endif
 
 
@@ -74 +75 @@
   !************************************************************
 
-  open(unitreceptor,file=path(1)(1:length(1))//'RECEPTORS',form='formatted',status='old',err=999)
+  open(unitreceptor,file=path(1)(1:length(1))//'RECEPTORS',form='formatted',&
+    status='old',err=999)
 
   ! try namelist input
-  read(unitreceptor,receptors,iostat=readerror)
+  read(unitreceptor,nml_receptors,iostat=ios)
 
   ! prepare namelist output if requested
-  if (nmlout.and.lroot) then
-    open(unitreceptorout,file=path(2)(1:length(2))//'RECEPTORS.namelist',&
-         &access='append',status='replace',err=1000)
-  endif
+  if (nmlout) open(unitreceptorout,file=path(2)(1:length(2))// &
+    'RECEPTORS.namelist', position='append',status='new',err=1000)
 
-  if ((lon.lt.-900).or.(readerror.ne.0)) then
+  if (ios .ne. 0) then ! read as regular text file
 
     close(unitreceptor)
-    open(unitreceptor,file=path(1)(1:length(1))//'RECEPTORS',status='old',err=999)
+    open(unitreceptor,file=path(1)(1:length(1))// &
+     'RECEPTORS',status='old',err=999)
+ 
     call skplin(5,unitreceptor)
 
@@ -101 +103 @@
     read(unitreceptor,'(4x,a16)',end=99) receptor
     call skplin(3,unitreceptor)
-    read(unitreceptor,'(4x,f11.4)',end=99) x
+      read(unitreceptor,'(4x,f11.4)',end=99) xlon
     call skplin(3,unitreceptor)
-    read(unitreceptor,'(4x,f11.4)',end=99) y
-    if ((x.eq.0.).and.(y.eq.0.).and. &
+      read(unitreceptor,'(4x,f11.4)',end=99) ylat
+      if ((xlon.eq.0.).and.(ylat.eq.0.).and. &
          (receptor.eq.'                ')) then
       j=j-1
@@ -116 +118 @@
     endif
     receptorname(j)=receptor
-    xreceptor(j)=(x-xlon0)/dx       ! transform to grid coordinates
-    yreceptor(j)=(y-ylat0)/dy
-    xm=r_earth*cos(y*pi/180.)*dx/180.*pi
+      xreceptor(j)=(xlon-xlon0)/dx       ! transform to grid coordinates
+      yreceptor(j)=(ylat-ylat0)/dy
+      xm=r_earth*cos(ylat*pi/180.)*dx/180.*pi
     ym=r_earth*dy/180.*pi
     receptorarea(j)=xm*ym
-
-    ! write receptors file in namelist format to output directory if requested
-    if (nmlout.and.lroot) then
-      lon=x
-      lat=y
-      write(unitreceptorout,nml=receptors)
-    endif
+  ! write receptors file in namelist format to output directory if requested
+    if (nmlout.and.lroot) write(unitreceptorout,nml=nml_receptors)
 
     goto 100
@@ -136 +133 @@
 
     j=0
-    do while (readerror.eq.0)
+    do while (ios .eq. 0)
       j=j+1
-      lon=-999.9
-      read(unitreceptor,receptors,iostat=readerror)
-      if ((lon.lt.-900).or.(readerror.ne.0)) then
-        readerror=1
-      else
+      read(unitreceptor,nml_receptors,iostat=ios)
+      if (ios .eq. 0) then
         if (j.gt.maxreceptor) then
           write(*,*) ' #### FLEXPART MODEL ERROR! TOO MANY RECEPTOR #### '
@@ -150 +144 @@
         endif
         receptorname(j)=receptor
-        xreceptor(j)=(lon-xlon0)/dx       ! transform to grid coordinates
-        yreceptor(j)=(lat-ylat0)/dy
-        xm=r_earth*cos(lat*pi/180.)*dx/180.*pi
+        xreceptor(j)=(xlon-xlon0)/dx       ! transform to grid coordinates
+        yreceptor(j)=(ylat-ylat0)/dy
+        xm=r_earth*cos(ylat*pi/180.)*dx/180.*pi
         ym=r_earth*dy/180.*pi
         receptorarea(j)=xm*ym
       endif
 
-      ! write receptors file in namelist format to output directory if requested
-      if (nmlout.and.lroot) then
-        write(unitreceptorout,nml=receptors)
-      endif
+    ! write receptors file in namelist format to output directory if requested
+      if (nmlout.and.lroot) write(unitreceptorout,nml=nml_receptors)
 
     end do
@@ -166 +158 @@
     close(unitreceptor)
 
-  endif
+  endif ! end no-nml / nml bloc
 
-  if (nmlout.and.lroot) then
-    close(unitreceptorout)
-  endif
+  if (nmlout .and. lroot) close(unitreceptorout)
 
   return
@@ -176 +166 @@
 
 999 write(*,*) ' #### FLEXPART MODEL ERROR! FILE "RECEPTORS"  #### '
-  write(*,*) ' #### CANNOT BE OPENED IN THE DIRECTORY       #### '
-  write(*,'(a)') path(1)(1:length(1))
+  write(*,*)   ' #### CANNOT BE OPENED IN THE DIRECTORY       #### '
+  write(*,'(a)') trim(path(1))
   stop
 
-1000 write(*,*) ' #### FLEXPART MODEL ERROR! FILE "RECEPTORS"    #### '
-  write(*,*) ' #### CANNOT BE OPENED IN THE DIRECTORY       #### '
-  write(*,'(a)') path(2)(1:length(2))
+1000 write(*,*) ' #### FLEXPART MODEL ERROR! FILE "RECEPTORS"  #### '
+  write(*,*)    ' #### CANNOT BE OPENED IN THE DIRECTORY       #### '
+  write(*,'(a)') trim(path(2))
   stop