Opened 6 years ago

Last modified 6 years ago

#197 closed Defect

Inconsistent use of ew in vertransform* — at Version 1

Reported by: adingwell Owned by:
Priority: major Milestone: FLEXPART 10
Component: FP coding/compilation Version: FLEXPART 10.0
Keywords: Cc:

Description (last modified by adingwell)

I found that ew is not properly initialized in verttransform.f90

These are all the occurrences of ew in each of the verttransform*-files:

verttransform.f90:

real :: ew,dz1,dz2,dz
...
tvold(ixm,jym)=tt2(ixm,jym,1,n)*(1.+0.378*ew*(td2(ixm,jym,1,n))/ &
     ps(ixm,jym,1,n))
...
tvold(ix,jy)=tt2(ix,jy,1,n)*(1.+0.378*ew(td2(ix,jy,1,n))/ &
     ps(ix,jy,1,n))

verttransform_gfs.f90:

real :: ew,pint,tv,tvold,pold,dz1,dz2,dz,ui,vi
...
tvold=tt2(ixm,jym,1,n)*(1.+0.378*ew(td2(ixm,jym,1,n))/ &
     ps(ixm,jym,1,n))

verttransform_nests.f90:

real :: ew,dz1,dz2,dz
...
tvold(ix,jy)=tt2n(ix,jy,1,n,l)*(1.+0.378*ew*(td2n(ix,jy,1,n,l))/ &
     psn(ix,jy,1,n,l))

Notice that in verttransform.f90 ew is used as a variable, not a function. This variable is defined but not initialized.

Change History (1)

comment:1 Changed 6 years ago by adingwell

  • Description modified (diff)
Note: See TracTickets for help on using tickets.
hosted by ZAMG