var_new = new((/ny,nx/),"float",getFillValue(u)) ; variável que conterá os valores máximos, tem que ser do
                                                 ; tipo var(lat,lon)
copy_VarCoords(u(0,:,:), var_new )

tempo = var_new                                  ; variável que conterá os valores máximos, tem que ser do
                                                 ; tipo tempo(lat,lon)
do lat = 0,ny-1
   do lon = 0,nx-1
      var              = u(:,lat,lon)            ; varredura de todos os tempos (:) para cada lat/lon
      max_var          = max(var)                ; obtenho para todos os tempos o valor máximo da variável em cada lat/lon.
      indMax          := ind(var.eq.max_var)     ; overwrite  :=  (6.1.2)
      index            = mod( indMax(0),12 ) + 1 ; obtenho o mês para cada valor máximo e somo 1 pq os meses começam em zero (índii
ces).
      var_new(lat,lon) = max_var                 ; preencho a variável com os máximos valores.
      tempo(lat,lon)   = index                   ; preencho o vetor do tempo com a variável "index" .

      if (dimsizes(indMax).gt.1) then            ; are there multiple places where var(:).eq.max_var?
          print("The following had "+dimsizes(indMax)+" index values = max_var")
      end if

      print("lat_i:  " + lat + "     lon_j:  " + lon + "     valor:  " + max_var + "    mes:  " + index)
   end do
end do