routing_irrigmap: nix and njx are not initialized when they do not belong to mpi root processor

[ajornet]

Under routing_irrigmap variables nix and njx are only initialized when they belong to mpi root processor. So they do not get a initial value. But later used printing its value.

7792	    IF (is_root_prc) THEN
7793	       nix=INT(MAXVAL(resolution_g(:,1))/MAXVAL(resol_lu(:,:,1)))+2
7794	       njx=INT(MAXVAL(resolution_g(:,2))/MAXVAL(resol_lu(:,:,2)))+2
7795	       nbpmax = nix*njx*2
7796	    ENDIF
7797	    CALL bcast(nbpmax)
7798	    !
7799	    callsign = 'Irrigation map'
7800	    ok_interpol = .FALSE.
7801	    !
7802	    WRITE(numout,*) "Projection arrays for ",callsign," : "
7803	    WRITE(numout,*) "nbpmax = ",nbpmax, nix, njx

Proposed solution:

7793	    nix=INT(MAXVAL(resolution_g(:,1))/MAXVAL(resol_lu(:,:,1)))+2
7794	    njx=INT(MAXVAL(resolution_g(:,2))/MAXVAL(resol_lu(:,:,2)))+2
7795	    nbpmax = nix*njx*2
7798	    !
7799	    callsign = 'Irrigation map'
7800	    ok_interpol = .FALSE.
7801	    !
7802	    WRITE(numout,*) "Projection arrays for ",callsign," : "
7803	    WRITE(numout,*) "nbpmax = ",nbpmax, nix, njx

This way it ensures nix and njx will get be initialized. It is not the most optimal proposal but it keeps consistency regarding routing_basins.

4236	    ! 
4237	    nix=INT(MAXVAL(resolution_g(:,1))/MAXVAL(resol_lu(:,:,1)))+2
4238	    njx=INT(MAXVAL(resolution_g(:,2))/MAXVAL(resol_lu(:,:,2)))+2
4239	    nbvmax = nix*njx*2
4240	    !
4241	    ! We are on the root processor here as this routine is not in parallel. So no need to broadcast.
4242	    !

[ajornet]

I applied this fix at r2948.

[jgipsl]

Done rev [2963] : write moved to IF is_root_prc