Custom Query (210 matches)

I have tried to run flexpart-wrf with different version(from 3.0 to 3.2). However, I found the valuse of potential voricity equal zeros, while I set particle dump for every output interval. Is any bugs in the programes? thanks for your help in advance.

BTW: There seems a bug in convmix_kfeta.f90 programe, if set variable "maxpar" not equal "numpart", which could make the value of "maxval(igrid)" too big and programe stop at the line of "too much x and y grid modify convmix_kfeta.f". This bug could be solved by add another loop behind the line 154.

do ipart=1,maxpart igrid(ipart)=-1 enddo

I've tried running FLEXPART-WRF with Mercator projection and got the error:

need orientation longitude (stdlon) as: 
   -180e <= lon1 <= 180w
STOP map_init

This error should not occur if Mercator projection is used. The problem is in the subroutine map_set. The variable 'proj_merc' is referenced before it is defined. I've marked the changes needed to fix the problem below:

      proj_code = proj_code_in
      
      !!! This is where map_init should be called to fix the problem
      call map_init

      ! first, check for validity of mandatory variables in proj
      if ( abs(lat1) .gt. 90.001 ) then
        print '(a)', 'latitude of origin corner required as follows:'
        print '(a)', '    -90n <= lat1 < = 90.n'
        stop 'map_init'
      endif
      if ( abs(lon1) .gt. 180.) then
        print '(a)', 'longitude of origin required as follows:'
        print '(a)', '   -180e <= lon1 <= 180w'
        stop 'map_init'
      endif
      if ((dx .le. 0.).and.(proj_code .ne. proj_latlon)) then
        print '(a)', 'require grid spacing (dx) in meters be positive!'
        stop 'map_init'
      endif
      if ((abs(stdlon) .gt. 180.).and.(proj_code .ne. proj_merc)) then
        print '(a,i3)', 'proj_code: ', proj_code
        print '(a,i3)', 'proj_merc: ', proj_merc
        print '(a)', 'need orientation longitude (stdlon) as: '
        print '(a)', '   -180e <= lon1 <= 180w'
        stop 'map_init'
      endif
      if (abs(truelat1).gt.90.) then
        print '(a)', 'set true latitude 1 for all projections!'
        stop 'map_init'
      endif

      !!! This is where "map_init" is called in 3.2 BETA

It's also rather confusing when the stop code is 'map_init' but the problem occurs in the subroutine 'map_set'.

Hi, When I compile the flexwrf3.2 beta with netcdf 4.0.1, the errors was accured as follows : pgf90 -c -I/home/weipeng/programfile/lib/netcdf4.0.1-pgi10/include -fastsse -mcmodel=medium write_ncheader.f90 PGF90-S-0038-Symbol, nf_netcdf4, has not been explicitly declared (write_ncheader.f90) PGF90-S-0038-Symbol, nf_def_var_deflate, has not been explicitly declared (write_ncheader.f90)

0 inform, 0 warnings, 2 severes, 0 fatal for write_ncheader

There was no errors while compiling the *.f file but the script stoped when making the excute file. I don't konw how to solve it and any helps will be appreciate. Thanks!