Subversion Repositories lagranto.ecmwf

Rev

Rev 39 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 39 Rev 44
Line 1972... Line 1972...
1972 
               traint(j,k,ncol+i) = mdv
 1972 
               traint(j,k,ncol+i) = mdv
1973 
            else
 1973 
            else
1974 
               call calc_HEAD (traint(j,k,ncol+i), &
 1974 
               call calc_HEAD (traint(j,k,ncol+i), &
1975 
                    traint(j,k  ,ind1),traint(j,k  ,ind2),traint(j,k+1,ind1),traint(j,k+1,ind2) )
 1975 
                    traint(j,k  ,ind1),traint(j,k  ,ind2),traint(j,k+1,ind1),traint(j,k+1,ind2) )
1976 
            endif
 1976 
            endif
- 
 
 1977
- 
 
 1978 
        ! Directional change (DANGLE)
- 
 
 1979 
         elseif  ( varsint(i+ncol).eq.'DANGLE' ) then
- 
 
 1980 
 
- 
 
 1981 
            varname='lon'
- 
 
 1982 
            call list2ind (ind1,varname,varsint,fok,ncol+ntrace1)
- 
 
 1983 
            varname='lat'
- 
 
 1984 
            call list2ind (ind2,varname,varsint,fok,ncol+ntrace1)
- 
 
 1985 
 
- 
 
 1986 
            if (k.eq.ntim) then
- 
 
 1987 
               traint(j,k,ncol+i) = mdv
- 
 
 1988 
            else
- 
 
 1989 
               if ( k.eq.1 ) then
- 
 
 1990 
                   traint(j,k,ncol+i) = mdv
- 
 
 1991 
               elseif ( k.eq.ntim ) then
- 
 
 1992 
                   traint(j,k,ncol+i) = mdv
- 
 
 1993 
               else
- 
 
 1994 
                   call calc_DANGLE (traint(j,k,ncol+i),      &
- 
 
 1995 
                       traint(j,k-1,ind1),traint(j,k-1,ind2), &
- 
 
 1996 
                       traint(j,k  ,ind1),traint(j,k  ,ind2), &
- 
 
 1997 
                       traint(j,k+1,ind1),traint(j,k+1,ind2) )
- 
 
 1998 
               endif
- 
 
 1999 
            endif
1977 
 
 2000 
 
1978 
 
 2001 
 
1979 
         ! Invalid tracing variable
 2002 
         ! Invalid tracing variable
1980 
         else
 2003 
         else
1981 
 
 2004