getfields Subroutine

public subroutine getfields(itime, nstop)

                                                                        *

This subroutine manages the 3 data fields to be kept in memory. * During the first time step of petterssen it has to be fulfilled that the * first data field must have |wftime|<itime, i.e. the absolute value of * wftime must be smaller than the absolute value of the current time in [s]. The other 2 fields are the next in time after the first one. * Pointers (memind) are used, because otherwise one would have to resort the wind fields, which costs a lot of computing time. Here only the pointers * are resorted. * * Author: A. Stohl * * 29 April 1994 * *

Changes, Bernd C. Krueger, Feb. 2001: * Variables tth,qvh,tthn,qvhn (on eta coordinates) in common block. * Function of nstop extended. * * Unified ECMWF and GFS builds * Marian Harustak, 12.5.2017 * - Added passing of metdata_format as it was needed by called routines *

                                                                        *

Variables: * lwindinterval [s] time difference between the two wind fields read in * indj indicates the number of the wind field to be read in * indmin remembers the number of wind fields already treated * memind(2) pointer, on which place the wind fields are stored * memtime(2) [s] times of the wind fields, which are kept in memory * itime [s] current time since start date of trajectory calcu- * lation * nstop > 0, if trajectory has to be terminated * nx,ny,nuvz,nwz field dimensions in x,y and z direction * uu(0:nxmax,0:nymax,nuvzmax,2) wind components in x-direction [m/s] * vv(0:nxmax,0:nymax,nuvzmax,2) wind components in y-direction [m/s] * ww(0:nxmax,0:nymax,nwzmax,2) wind components in z-direction [deltaeta/s]* tt(0:nxmax,0:nymax,nuvzmax,2) temperature [K] * ps(0:nxmax,0:nymax,2) surface pressure [Pa] * metdata_format format of metdata (ecmwf/gfs) * * Constants: * idiffmax maximum allowable time difference between 2 wind * fields * *

Arguments