Bonjour, une maintenance Gitlab aura lieu mardi prochain le 6 octobre afin de mettre à jour la version de Gitlab ainsi que d'ajouter la possibilité d'utiliser votre compte UQAR pour la connexion au serveur. Un courriel avec plus de détails vous sera envoyé une fois la maintenance complétée. Si vous avez des questions ou des commentaires, veuillez communiquer à l'adresse admin_calcul@uqar.ca. Merci de votre collaboration.

write_output.f90 1.97 KB
Newer Older
Jérémy Baudry's avatar
Jérémy Baudry committed
1 2 3 4 5 6 7
subroutine write_output

use netcdf
use parameters
implicit none
include 'netcdf.inc'
	integer :: ncid,omID,xID,tID,SpectreID
Jérémy Baudry's avatar
Jérémy Baudry committed
8 9
	integer :: om_varID,x_varID,t_varID,stat
	integer :: Dmax_varID,Dave_varID
Jérémy Baudry's avatar
Jérémy Baudry committed
10

Jérémy Baudry's avatar
Jérémy Baudry committed
11 12
stat=nf90_create(namefile,nf90_clobber,ncid)
call handle_err(stat)
Jérémy Baudry's avatar
Jérémy Baudry committed
13

Jérémy Baudry's avatar
Jérémy Baudry committed
14 15
        stat=nf90_def_dim(ncid,"omega", nfreq, omID)
        call handle_err(stat)
Jérémy Baudry's avatar
Jérémy Baudry committed
16

Jérémy Baudry's avatar
Jérémy Baudry committed
17 18 19 20
        stat=nf90_def_dim(ncid,"x_axis", nbin, xID)
        call handle_err(stat)
        stat=nf90_def_dim(ncid,"time", nsteps, tID)
        call handle_err(stat)
Jérémy Baudry's avatar
Jérémy Baudry committed
21

Jérémy Baudry's avatar
Jérémy Baudry committed
22 23 24 25 26 27 28 29 30 31 32 33
        stat=nf90_def_var(ncid,"omega",nf90_double,omID,om_varID)
        call handle_err(stat)
        stat=nf90_def_var(ncid,"x_axis",nf90_double,xID,x_varID)
        call handle_err(stat)
        stat=nf90_def_var(ncid,"time",nf90_double,tID,t_varID)
        call handle_err(stat)
        stat=nf90_def_var(ncid,"Spectrum",nf90_double,(/ tID, xID, omID /),SpectreID)
        call handle_err(stat)
        stat=nf90_def_var(ncid,"Dave",nf90_double,xID,Dave_varID)
        call handle_err(stat)
        stat=nf90_def_var(ncid,"Dmax",nf90_double,xID,Dmax_varID)
        call handle_err(stat)
Jérémy Baudry's avatar
Jérémy Baudry committed
34

Jérémy Baudry's avatar
Jérémy Baudry committed
35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67

        !stat=nf90_put_att(ncid,ID,"units", values)
        stat=nf90_enddef(ncid)
        call handle_err(stat)
        stat=nf90_put_var(ncid,om_varID,omega)
        call handle_err(stat)
        stat=nf90_put_var(ncid,spectreID,E)
        call handle_err(stat)
        stat=nf90_put_var(ncid,Dmax_varID,Dmax)
        call handle_err(stat)
        stat=nf90_put_var(ncid,Dave_varID,Dave)
        call handle_err(stat)
        stat=nf90_put_var(ncid,t_varID,time)
        call handle_err(stat)
        stat=nf90_put_var(ncid,x_varID,x_axis)
        call handle_err(stat)




stat=nf90_close(ncid)
call handle_err(stat)
        
contains

subroutine handle_err(stat)
        integer, intent(in) :: stat

        if(stat /=nf90_noerr) then
        write(*,*)trim(nf90_strerror(stat))
        STOP 'unable de write the netcdf output' 
        end if
end subroutine handle_err
Jérémy Baudry's avatar
Jérémy Baudry committed
68 69 70 71



end subroutine write_output