NOTAS DE AULA DO CURSO
1. INSTALAÇÃO DO MODELO WW3
1.1 PROGRAMAS NECESSÁRIOS
–
–
–
–
–
Compilador necessário do fortran: IFORT, PGF90, GFORTRAN, etc.
Sistema operacional : UNIX/Linux
Software para visualização: GRADS ou Matlab ou Phyton, etc
Cadastrar-se no site do NCEP
https://lstsrv.ncep.noaa.gov/mailman/listinfo/ncep.list.wwatch3.users
Baixar o modelo wwatch3.v4.18.tar.gz da página
http://polar.ncep.noaa.gov/waves/wavewatch/distribution com usuário e senha cadastrado no
site:
Index of /waves/wavewatch/distribution
Name
Last modified
Size
Parent Directory
WW3_GMD_genetics.v1.4.tar.gz
gridgen.v2.1.0.tar.gz
gridgen1.0/
gridgen1.1.tar.gz
gridgen3.0.tar.gz
v3.14-14/
v3.14/
v4.18/
19-Mar-2014
29-Apr-2013
09-Jun-2009
27-Jan-2011
01-Aug-2014
02-Oct-2013
07-Aug-2009
21-Mar-2014
7.0M
555M
555M
556M
-
16:19
18:56
17:51
15:42
20:23
20:47
14:57
13:59
Na pasta v.4.18 encontra-se wwatch3.v4.18.tar.gz
1.2 INSTALAÇÃO
Criar uma pasta para instalar o modelo:
> mkdir ww3v418
Descompactar o programa do modelo na pasta ww3v418:
> tar -zxvf wwatch3.v4.18.tar.gz
(Não jogar o original fora, por se dar errado começar do zero)
–
–
> ls -l ww3v418
total 464380
-rw------- 1 rosio users 217262 Mar 19 21:17 guide.v4.18.pdf
-rwx------ 1 rosio users 39546 Mar 19 21:18 install_ww3_tar
-rw------- 1 rosio users 3542701 Mar 19 21:18 manual.v4.18.pdf
-rw------- 1 rosio users 69949440 Mar 19 21:18 wwatch3.v4.18.cases.tar
-rw------- 1 rosio users 135895040 Mar 19 21:18 wwatch3.v4.18.model.tar
-rw------- 1 rosio users 123146240 Mar 19 21:18 wwatch3.v4.18.regtests.tar
-rw-r--r-- 1 rosio users 142723764 Ago 19 16:43 wwatch3.v4.18.tar.gz
–
Verificar as permissões dos arquivos, especialmente o install_ww3_tar ou digitar o comando
:
> chmod -R +x * (libera a permissão de execução)
> ls -l
total 464380
-rwx--x--x 1 rosio users 217262 Mar 19 21:17 guide.v4.18.pdf
-rwx--x--x 1 rosio users 39546 Mar 19 21:18 install_ww3_tar
-rwx--x--x 1 rosio users 3542701 Mar 19 21:18 manual.v4.18.pdf
-rwx--x--x 1 rosio users 69949440 Mar 19 21:18 wwatch3.v4.18.cases.tar
-rwx--x--x 1 rosio users 135895040 Mar 19 21:18 wwatch3.v4.18.model.tar
-rwx--x--x 1 rosio users 123146240 Mar 19 21:18 wwatch3.v4.18.regtests.tar
-rwxr-xr-x 1 rosio users 142723764 Ago 19 16:43 wwatch3.v4.18.tar.gz
–
Executar :
> ./install_ww3_tar
Aparece na tela:
===================================
------ Installing WAVEWATCH III v.4 -----===================================
Script for installing package from tar files.
Requires files in same directory as script.
Continue? [y|n]y
===================================
------ Installing WAVEWATCH III v.4 -----===================================
from tar source
This installation requires a configuration file (wwatch3.env).
The current version allows two types of env files:
- A local [L] wwatch3.env (Allowing multiple independent installations).
- A generic [G] dot-file .wwatch3.env (Old-fashioned option).
[L] Installs new, uses existing or updates env file in current directory.
[G] Installs new, uses existing or updates env file in home directory,
(home is presumably /home2/rosio}).
Type your choice now: L
Installing in
/home2/rosio/ww3v418
OK ? [y/n] y
Setting up environment variables.
Previous setup file not found. Variables will be set to defaults.
(User must check to see if these setting are appropriate.)
Creating wwatch3.env locally (also in home if G option chosen).
Printer (listings)
: printer
FORTRAN comp. (aux only) : f77
C Compiler (aux only) : cc
Scratch directory
: /home2/rosio/ww3v418/tmp
Save source code
: yes
Save listings
: yes
Update settings ? [y/n]y
Modifying set-up
Type n new settings, or press ENTER to keep [current ones]:
Printer for listings [printer] :
Compiler for aux. [f77] : ifort
Compiler for aux. [cc] : cc
Scratch space [/home2/rosio/ww3v418/tmp] :
Save source code files (*.f) [yes] :
Save listing files [yes] :
Modified settings:
Printer (listings)
: printer
FORTRAN comp. (aux only) : ifort
C Compiler (aux only) : cc
Scratch directory
: /home2/rosio/ww3v418/tmp
Save sources
: yes
Save listings
: yes
New settings OK ? [y/n] y
Continue with actual implementation ? [y/n]y
- No final da instalação aparecem as seguintes pastas:
> ls -l
total 3728
drwxr-xr-x 2 rosio users 4096 Ago 19 16:57 arc
drwxr-xr-x 6 rosio users 4096 Ago 19 16:56 aux
drwxr-xr-x 2 rosio users 4096 Ago 19 16:56 bin
drwxr-xr-x 6 rosio users 4096 Ago 19 16:57 cases
drwxr-xr-x 2 rosio users 4096 Ago 19 16:56 exe
drwxr-xr-x 3 rosio users 4096 Ago 19 16:56 ftn
-rwx--x--x 1 rosio users 217262 Mar 19 21:17 guide.v4.18.pdf
drwxr-xr-x 2 rosio users 4096 Ago 19 16:56 inp
lrwxrwxrwx 1 rosio users
21 Ago 19 16:57 install_ww3_tar -> ./bin/install_ww3_tar
-rwx--x--x 1 rosio users 3542701 Mar 19 21:18 manual.v4.18.pdf
drwxr-xr-x 2 rosio users 4096 Ago 19 16:56 mod
drwxr-xr-x 2 rosio users 4096 Ago 19 16:56 obj
drwxr-xr-x 40 rosio users 4096 Ago 19 16:57 regtests
drwxr-xr-x 2 rosio users 4096 Ago 19 16:57 tmp
drwxr-xr-x 2 rosio users 4096 Ago 19 16:57 work
-rw-r--r-- 1 rosio users 282 Ago 19 16:56 wwatch3.env
2. COMPILAÇÃO
2.1 Entrar na pasta /home2/rosio/ww3v418/bin:
> cd bin
2.2 Antes de compilar: copiar os arquivos comp e link segundo os compiladores do fortran
disponível (Intel, Portland, gfortran, etc); modificar o switch.
comp
Para escolher o compilador da Intel (ifort):
> cp comp.Intel comp
link
Para escolher o link da Intel:
> cp link.Intel link
* Abrir os arquivos comp e link e editar as seguintes linhas:
– comp
trocar mpiifort para mpif90 nas linhas abaixo:
linha 98 comp=mpiifort
linha 99 which mpiifort 1> /dev/null 2> /dev/null
–
link
trocar mpiifort para mpif90 nas linhas abaixo:
linha 105 comp=mpiifort
linha 106 which mpiifort 1> /dev/null 2> /dev/null
switch
Para compilar o modelo em série, o arquivo switch deve conter o flag SHRD:
> vi switch
F90 NOGRB NOPA LRB4 SHRD PR3 UQ FLX2 LN1 ST2 STAB2 NL1 BT1 DB1 MLIM TR0
BS0 IC0 REF0 XX0 WNT1 WNX1 CRT1 CRX1 O0 O1 O2 O3 O4 O5 O6 O7 O11 O14
Para compilar o modelo em paralelo, o arquivo switch deve conter o flag MPI DIST :
> vi switch
F90 NOGRB NOPA LRB4 MPI DIST PR3 UQ FLX2 LN1 ST2 STAB2 NL1 BT1 DB1 MLIM
TR0 BS0 IC0 REF0 XX0 WNT1 WNX1 CRT1 CRX1 O0 O1 O2 O3 O4 O5 O6 O7 O11 O14
2.3 Exportar os caminhos dos diretórios
gedit runpath.sh
#!/bin/bashrc
export PATH=$PATH:/(local onde está a pasta ww3v418)/ww3v418
export PATH=$PATH:/(local onde está a pasta ww3v418)/ww3v418/bin
export PATH=$PATH:/(local onde está a pasta ww3v418)/ww3v418/exe
export WWATCH3_NETCDF=NC3
export NETCDF_CONFIG=/(local onde está o netcdf)/netcdf/bin/
export NETCDF_LIBDIR=/(local onde está o netcdf)/netcdf/lib/
export NETCDF_INCDIR=/(local onde está o netcdf)/netcdf/include/
chmod +x runpath.sh
source runpath.sh
Executar w3_make
–
para ver se os arquivos foram criados na compilação:
ls exe
total 57328
-rwxrwxr-x. 1 luis luis 2520943 Sep 4 10:15 ww3_grid
-rwxrwxr-x. 1 luis luis 2347979 Sep 4 10:15 ww3_strt
-rwxrwxr-x. 1 luis luis 2226715 Sep 4 10:15 ww3_bound
-rwxrwxr-x. 1 luis luis 2384330 Sep 4 10:15 ww3_prep
-rwxrwxr-x. 1 luis luis 3357224 Sep 4 10:16 ww3_shel
-rwxrwxr-x. 1 luis luis 2796824 Sep 4 10:16 ww3_outf
-rwxrwxr-x. 1 luis luis 2425957 Sep 4 10:16 ww3_outp
-rwxrwxr-x. 1 luis luis 1565437 Sep 4 10:16 ww3_trck
-rwxrwxr-x. 1 luis luis 2684179 Sep 4 10:16 ww3_grib
-rwxrwxr-x. 1 luis luis 2362509 Sep 4 10:16 ww3_gspl
-rwxrwxr-x. 1 luis luis 2709094 Sep 4 10:17 ww3_gint
-rwxrwxr-x. 1 luis luis 2718820 Sep 4 10:17 gx_outf
-rwxrwxr-x. 1 luis luis 2302981 Sep 4 10:17 gx_outp
-rwxrwxr-x. 1 luis luis 3863479 Sep 4 10:18 ww3_multi
-rwxrwxr-x. 1 luis luis 3873243 Sep 4 10:18 ww3_sbs1
-rwxrwxr-x. 1 luis luis 4442281 Sep 4 10:18 ww3_prnc
-rwxrwxr-x. 1 luis luis 5133278 Sep 4 10:18 ww3_ounf
-rwxrwxr-x. 1 luis luis 4635242 Sep 4 10:18 ww3_ounp
-rwxrwxr-x. 1 luis luis 4307415 Sep 4 10:18 ww3_bounc
-rw-rw-r--. 1 luis luis 283 Sep 4 10:18 exec_type
3. TESTANDO O EXEMPLO PADRÃO (MODELO WW3)
SERIAL
•
para rodar o exemplo padrão:
cd ww3_serial/work
../exe/ww3_grid
- irá gerar o mod_def.ww3
../exe/ww3_strt
- irá gerar o restart
../exe/ww3_prep
- irá criar o ice.ww3 e o wind.ww3
../exe/ww3_shel
•
para gerar as saídas das variáveis:
../exe/ww3_outf
../exe/ww3_outp
../exe/ww3_ounf
4. TESTANDO O EXEMPLO DO CURSO (MODELO WW3) – SERIAL
•
•
Mapear o python e o matlab no bashrc:
gedit ~/.bashrc
export PATH=$PATH:/(local onde está o python)/anaconda
export PATH=$PATH:/(local onde está o python)/anaconda/bin
export PATH=$PATH:/(local onde está o python)/anaconda/pkgs
export PATH=$PATH:/(local onde está o python)/anaconda/include
export PATH=$PATH:/(local onde está o python)/anaconda/lib
#
alias py='/(local onde está o python)/anaconda/bin/ipython'
alias python='/(local onde está o python)/anaconda/bin/ipython'
alias mat='matlab -nosplash -nodesktop'
•
source ~/.bashrc
4.1 Gerar as grades:
cd ww3v418
mkdir gridgen
Para baixar o gridgen é necessário ter usuário e senha para acessar o NCEP/NOAA, para isto basta
acessar o site http://polar.ncep.noaa.gov/waves/wavewatch/license.shtml e clicar em: “I agree to
comply with the license terms.”, ao final da página, e solicitar usuário e senha.
Para baixar o gridgen3.0.tar.gz, acessar:
http://polar.ncep.noaa.gov/waves/wavewatch/distribution/
usuario:
senha:
copiar o gridgen3.0.tar.gz para esta pasta e descompactá-lo:
tar -xpvf gridgen3.0.tar.gz
Fazer o download do arquivo do curso: arquivos.tar.gz
- Extrair o arquivos.tar.gz
cd examples
copiar os scripts que conterão as informações sobre a grade a ser formada para
a pasta examples:
cp arquivos/dia2/create_grid.* .
(pasta arquivos disponível para Dowload no site)
•
• para criar a grade global:
cp create_grid.m cria_grade_global.m
cd gridgen
mkdir output
cd examples
• modificar as seguintes informações no cria_grade_global.m:
gedit cria_grade_global.m
% Path to directories
bin_dir='/(local onde está o gridgen)/gridgen/bin'
ref_dir='/(local onde está o gridgen)/gridgen/reference_data'
out_dir='/(local onde está o gridgen)/gridgen/output'
% Grid definition
dx=2.0
dy=2.0
lon1d=[0:dx:360]
lat1d=[-80:dy:80]
fname='global'
boundary='low'
*ver melhor tipo para fronteira
ref_grid='etopo1' ou 'etopo2'
*ver melhor dado de batimetria
LAKE_TOL= 1
IS_GLOBAL=1
*quando o domínio é global =1, com domínio regional = 0
• usar o matlab para gerar a malha global:
mat < cria_grade_global.m
•
os dados gerados estão na pasta output
• para criar a grade regional:
cp create_grid.m cria_grade_southatl.m
• modificar as seguintes informações no cria_grade_southatl.m:
gedit cria_grade_southatl.m
% Path to directories
bin_dir='/(local onde está o gridgen)/gridgen/bin'
ref_dir='/(local onde está o gridgen)/gridgen/reference_data'
out_dir='/(local onde está o gridgen)/gridgen/output'
% Grid definition
dx=0.5
dy=0.5
lon1d=[-69:dx:-6]
lat1d=[291:dy:358]
fname='southatl'
boundary='low'
ref_grid='etopo1' ou 'etopo2'
LAKE_TOL= 1
IS_GLOBAL=0
% ver melhor tipo para fronteira
% ver melhor dado de batimetria
% quando o domínio é global =1, com domínio regional = 0
•
usar o matlab para gerar a malha regional:
mat < cria_grade_southatl.m
•
•
os dados gerados estão na pasta output
criar a pasta ww3_curso dentro da pasta ww3v418:
cd ww3v418
mkdir ww3_curso
cd ww3_curso
mkdir fix
cd fix
•
copiar os resultados da pasta output para a pasta fix:
cp ../../gridgen/output/global.* .
cp ../../gridgen/output/southatl.* .
•
para a grade global:
mv global.meta ww3_grid.inp.global
•
modificar as seguintes informações no ww3_grid.inp.global:
gedit ww3_grid.inp.global
* há um exemplo desse arquivo .inp na pasta ww3_serial/inp (/ww3_grid.inp)
contendo as informações de cada flag
$ Grid name
'GRADE CURSO INPE GLOBAL 2X2'
$ Frequency increment factor and first frequency (Hz)
1.1 0.05 25 24 0.5
$ Set model flags
FTTTFT
$ Set time steps
3600. 1800. 1800. 30.
$ Start of namelist input section
&MISC CICE0 = 0.25, CICEN = 0.75, FLAGTR = 4 /
END OF NAMELISTS
$ Define grid
'RECT' T 'SMPL'
181 81
120.00 120.00 60.00
0.0000 -80.0000 1.00
$ Bottom Bathymetry
-0.10 2.50 40 0.001000 1 1 '(....)' 'NAME' 'global.depth_ascii'
$ Sub-grid information
50 0.010000 1 1 '(....)' 'NAME' 'global.obstr_lev1'
$ Mask Information
60 1 1 '(....)' 'NAME' 'global.maskorig_ascii'
•
criar o link:
ln -s ww3_grid.inp.global ww3_grid.inp
•
executar o ww3_grid:
(local onde está o ww3_serial)/ww3_serial/exe/ww3_grid
mv mod_def.ww3 mod_def.global
•
para a grade regional:
mv southatl.meta ww3_grid.inp.southatl
cp arquivos/dia2/maskFinalAtlSul.txt .
modificar as seguintes informações no ww3_grid.inp.southatl :
gedit ww3_grid.inp. southatl
* há um exemplo desse arquivo .inp na pasta ww3_serial/inp (/ww3_grid.inp) contendo
as informações de cada flag
$ Grid name
'GRADE CURSO INPE SOUTH ATLANTIC'
$ Frequency increment factor and first frequency (Hz)
1.1 0.05 25 24 0.5
$ Set model flags
FTTTTT
$ Set time steps
900. 450. 450. 30.
$ Start of namelist input section
&MISC CICE0 = 0.25, CICEN = 0.75, FLAGTR = 4 /
END OF NAMELISTS
$ Define grid
'RECT' T 'NONE'
135 127
30.00 30.00 60.0
291.0000
-69.0000
1.00
$ Bottom Bathymetry
-0.10 10.0 20 0.001000 1 1 '(....)' 'NAME' 'southatl.depth_ascii'
$ Sub-grid information
30 0.010000 1 1 '(....)' 'NAME' 'southatl.obstr_lev1'
$ Mask Information
40 1 1 '(....)' 'NAME' 'maskFinalAtlSul.txt'
•
•
remover e criar o link:
rm ww3_grid.inp
ln -s ww3_grid.inp.southatl ww3_grid.inp
•
executar o ww3_grid:
(local onde está o ww3_serial)/ww3_serial/exe/ww3_grid
mv mod_def.ww3 mod_def.southatl
4.2 Gerar as forçantes:
•
para gerar a grade para as forçantes:
cp arquivos/dia2/ww3_grid.inpe.cfsr .
rm ww3_grid.inp
ln -s ww3_grid.inp.cfsr ww3_grid.inp
•
executar o ww3_grid:
(local onde está o ww3_serial)/ww3_serial/exe/ww3_grid
mv mod_def.ww3 mod_def.cfsr
•
para gerar os arquivos com as forçantes vento e gelo:
cp /home/curso01/arquivos/dia3/ww3_prnc.inpe.icecfsr .
cp /home/curso01/arquivos/dia3/ww3_prnc.inpe.wndcfsr .
•
editar o ww3_prnc.inpe.wndcfsr:
gedit ww3_prnc.inpe.wndcfsr
- alterar o nome do arquivo para: wnd10m.gdas.200304.nc
•
editar o ww3_prnc.inpe.icecfsr:
gedit ww3_prnc.inpe.icecfsr
- alterar o nome do arquivo para: icecon.gdas.200304.nc
•
cria a pasta cfsr no ww3_curso
cd ww3_curso
mkdir cfsr
cd cfsr
•
criar os seguinte links para o vento:
ln -s /home/curso01/arquivos/dia3/vento/wnd10m.gdas.200304.nc
ln -s ../fix/ww3_prnc.inpe.wndcfsr ww3_prnc.inp
ln -s ../fix/mod_def.cfsr mod_def.ww3
•
executar o ww3_prnc para o vento:
(local onde está o ww3_serial)/ww3_serial/exe/ww3_prnc
mv wnd.ww3 wnd.cfsr
rm ww3_prnc.inp
•
criar os seguinte links para o gelo:
ln -s /home/curso01/arquivos/dia3/vento/icecon.gdas.200304.nc
ln -s ../fix/ww3_prnc.inpe.icecfsr ww3_prnc.inp
•
executar o ww3_prnc para o gelo:
(local onde está o ww3_serial)/ww3_serial/exe/ww3_prnc
mv ice.ww3 ice.cfsr
4.3 Passos finais para executar o modelo:
•
criar a pasta work no ww3_curso
cd ww3_curso
mkdir work
cd work
•
trazer para a pasta work os arquivos necessário para a rodada:
cp arquivos/dia3/ww3_multi.inp.umdia .
•
editar o ww3_multi.inp.umdia:
gedit ww3_multi.inp.umdia (* duvidas, acessar ww3_multi.inp na pasta INP do modelo)
$
21F1TT
$
* Caso queira visualizar os resultados em pontos específicos, inserir no arquivo
ww3_multi.inp as linhas sombreadas em cinza, na seguinte ordem:
$ WAVEWATCH III multi-scale input file
$ -----------------------------------21F1TT
'cfsr' F F T T F F F
$
'global' 'no' 'no' 'cfsr' 'cfsr' 'no' 'no' 'no' 1 1 0.00 1.00 F
'southatl' 'no' 'no' 'cfsr' 'cfsr' 'no' 'no' 'no' 2 1 0.00 1.00 F
$
20030401 000000 20030402 000000
$
F F
$
20030401 000000 3600 20030402 000000
N
WND HS FP DP
$
20030401 000000 3600 20030402 000000
$
-49.86667 -31.56667
'p101_boiRS'
-43.46000 -23.16000
'p504_RioJa'
$
0.0
0.0
'STOPSTRING'
$
20030401 000000
0 20030402 000000
20030402 000000
1 20030402 000000
20030401 000000
0 20030402 000000
20030401 000000
0 20030402 000000
$
'the_end' 0
'STP'
$
ln -s ww3_multi.inp.umdia ww3_multi.inp
ln -s ../fix/mod_def.cfsr mod_def.cfsr
ln -s ../fix/mod_def.global mod_def.global
ln -s ../fix/mod_def.southatl mod_def.southatl
ln -s ../cfsr/wnd.cfsr wnd.cfsr
ln -s ../cfsr/ice.cfsr ice.cfsr
•
executar o ww3_multi:
(local onde está o ww3_serial)/ww3_serial/exe/ww3_multi
•
criar a pasta teste1 no ww3_curso para colocar os resultados
cd ww3_curso
mkdir teste1
cd teste1
mv out_* .
mv restart* .
mv log.* .
cp /home/curso01/arquivos/dia3/ww3_ounf.inp .
4.4 Visualizar os resultados:
•
editar o ww3_ounf.inp:
gedit ww3_ounf.inp (* duvidas, acessar ww3_ounf.inp na pasta INP do modelo)
$
20030401 000000 3600 9999
$
•
para o global, criar os links:
ln -s ../fix/mod_def.global mod_def.ww3
ln -s out_grd.global out_grd.ww3
•
executar o ww3_ounf
(local onde está o ww3_serial)/ww3_serial/exe/ww3_ounf
rm mod_def.ww3 out_grd.ww3
mkdir global
mv *.nc global/
•
para o regional, criar os links:
ln -s ../fix/mod_def.southatl mod_def.ww3
ln -s out_grd.southatl out_grd.ww3
•
executar o ww3_ounf:
(local onde está o ww3_serial)/ww3_serial/exe/ww3_ounf
mkdir regional
mv *.nc regional/
5. TESTANDO O EXEMPLO DO CURSO (MODELO WW3) – PARALELO
export PATH=$PATH:/(pasta ww3_mpi)/d-WW3_mpi
export PATH=$PATH:/(pasta ww3_mpi)/d-WW3_mpi/bin
export PATH=$PATH:/(pasta ww3_mpi)/d-WW3_mpi/exe
export WWATCH3_NETCDF=NC3
export NETCDF_CONFIG=/opt/netcdf/3.6.3/intel/bin/
export NETCDF_LIBDIR=/opt/netcdf/3.6.3/intel/lib/
export NETCDF_INCDIR=/opt/netcdf/3.6.3/intel/include/
•
Para rodar o modelo em paralelo executar os passos dos itens 4.1, 4.2 e 4.3 (exceto o
passo da execução).
• Criar o arquivo machines para o seu cluster (exemplo):
robura:24
•
Executando o ww3_multi (Paralelo):
Sintaxe:
(pasta do mpich2)/mpirun -machinefile (file_machine) -np (num_processadores) ww3_multi
Exemplo:
/opt/mpich2/mpirun -machinefile machines -np 12 ww3_multi
•
Para visualizar os resultados executar o passo 4.4.
6. EXTRAINDO PONTOS DOS RESULTADOS DO MODELO E COMPARANDO
COM DUAS BOIAS (EXEMPLO SERIAL)
• Acessar a pasta teste1, copiar e linkar:
cd teste1
cp /arquivos/teste_ST4/ww3_ountp.inp .
ln -s out_pnt.southatl out_pnt.ww3
• Executar:
ww3_outp
• Criar pasta teste2
mkdir (dir_local)/ww3_curso/teste2
cp (dir_local)/ww3_curso/teste2/work/out* .
mv teste2 teste_ST4
mkdir teste_ST2
cp ../teste1/ww3_altnf.inp
cp ../fix/mod_def.southaltl
mv mod_def.southaltl mode_def.ww3
ln -s out_grd.doutjatl out_grd.ww3
•
Executar ww3_ounf.
ww3_ountf
•
Copiar e linkar
cp (pasta raiz)/ww3_curso/teste_ST4/ww3_outp.inp
ln -s /arquivos/ww3_curso/teste_ST4/tspoint_wndwave.py
mv tab50.ww3 southatl.tab
•
Executar python
python tspoint_wndwall.py
•
Copiar mod_def.global
cp ../fix/mod_def.global .
mv mod_def.global mod_def.ww3
ln -s out_grd.global out_grd.ww3
•
Executar ww3_ounf
ww3_ountf
•
Executar ww3_outp
ln -s out_pnt.global out_pnt.ww3
ww3_outp
•
Mover também os pointers criados – pointparameters_001_southatl.txt
cp /arquivos/dia4/wave_AxysRS.txt .
py
from pylab import *
import numpy as np
ls
hsb=np.loadtxt('wave_AxysRS.txt')
hsb
hs=hsb[:,5]
ls
hsm_satl=np.loadtxt('satl_pntpar_001.txt')[:,5]
hsm_glo=np.loadtxt('glo_pntpar_001.txt')[:,5]
plot(hsb,'x',hsm_glo,'r',hsm_satl, 'b'; show ())
cd ../teste_ST2
ln -s '/arquivos/teste_andre/out_*
* repetir os passoas acima
•
Teste de Regressão
cd ..
mkdir regressao
cp /arquivos/dia4/resultComp_AxysMinuanoRS.txt .
mat (este comando chama o matlab)
file1 = load('resultComp_AxysMinuanoRS.txt')
vcfsr = file1(:,8);
x = [1:1453];
vboia = file1(:,15);
plot(x,vcfsr,'r',vboia,'b')
cp /arquivos/dia4/regressao.m .
cp /arquivos/dia4/aplicaRegLinCFSR.py .
cp /arquivos/dia4/wnd10mRegLin.200304_teste.nc .
cd cfsr
rm wnd10.2003.nc
ln -s /arquivos/dia4/wnd10mRegLin.200304_teste.nc wnd10m.gdas.200304.nc
mv wind.cfsr wind_cfsr.old
cp /arquivos/dia2/ww3_prnc.inp.cfsr
more ww3_prnc.inp (conferir se as duas linhas abaixo estao corretas)
uwnd vwnd
'wnd10m.gdas.200304.nc'
•
Executar
ww3_prnc
ln -s wind.ww3 wind.cfsr
cd ../work
ln -s ./cfsr/wind.ww3 wind.cfsr
mpd &
(pasta_mpich2)/mpirun -machinefile machine -np 2 ww3_multi