source: branches/iLoveclim/SOURCES/no_spinup_mod.f90 @ 91

Last change on this file since 91 was 77, checked in by dumas, 8 years ago

Merge branche iLOVECLIM sur rev 76

File size: 2.9 KB
Line 
1!> \file no_lakes.f90
2!! This module is to choosen when no spinup is needed
3!<
4
5!> \namespace no_lakes
6!! This module is to choosen when no spinup is needed
7!! \author ...
8!! \date ...
9!! @note  used modules
10!! @note  - use module3D_phy
11!<
12
13
14module no_spinup
15
16use module3d_phy
17
18contains
19!-------------------------------------------------------------------------------------
20!> SUBROUTINE: init_spinup
21!! Routine permettant d'initialiser les spinup
22!>
23subroutine init_spinup
24
25implicit none
26
27namelist/spinup/ispinup
28
29    ! put here which type of spinup
30    ! ispinup = 0                          ! run standard
31    ! ispinup = 1                          ! ispinup = 1 -> saute icethick3, diagnoshelf,
32    ! diffusiv, isostasie   pour equilibre temperature
33    ! en prenant les vitesses calculees aux premiers
34    ! pas de temps 
35    ! ispinup = 2                          ! ispinup = 2 -> fait la conservation
36    !                                      ! de la masse avec les vitesses de bilan
37    ! ispinup = 3                          ! fait le calcul des temperatures
38    !                                      ! avec les vitesses de bilan
39 
40
41
42! lecture des parametres
43
44rewind(num_param)        ! pour revenir au debut du fichier param_list.dat
45read(num_param,spinup)
46
47write(num_rep_42,*)'!___________________________________________________________'
48write(num_rep_42,*)'!        spinup          module no_spinup               '
49write(num_rep_42,spinup)                   ! pour ecrire les parametres lus
50write(num_rep_42,*)
51write(num_rep_42,*)
52write(num_rep_42,*)'! ispinup = 0     run standard'
53write(num_rep_42,*)'! ispinup = 1     equilibre temperature avec vitesses grisli '
54write(num_rep_42,*)'! ispinup = 2     conservation de la masse avec vitesses bilan '
55write(num_rep_42,*)'! ispinup = 3     equilibre temperature avec vitesses bilan'
56write(num_rep_42,*)
57
58 
59if (ispinup.ge.2) then
60   write(6,*)' ispinup > 1 must be called with module spinu_vitbil '
61   write(6,*) 'and dragging modules must be removed (see module_choix)'
62   stop
63endif
64
65
66if (itracebug.eq.1)  call tracebug(' fin routine init_spinup de no_spinup')
67
68end subroutine init_spinup
69!-------------------------------------------------------------------------------------
70subroutine force_balance_vel
71
72if (itracebug.eq.1)  call tracebug(' Subroutine force_balance_vel dans no_spinup')
73
74end subroutine force_balance_vel
75!-------------------------------------------------------------------------------------
76subroutine limit_coef_vitbil
77
78if (itracebug.eq.1)  call tracebug(' Subroutine limit_coef_vitbil dans no_spinup')
79
80end subroutine limit_coef_vitbil
81!-------------------------------------------------------------------------------------
82subroutine calc_coef_vitbil
83
84if (itracebug.eq.1)  call tracebug(' Subroutine limit_coef_vitbil dans no_spinup')
85
86end subroutine calc_coef_vitbil
87!-------------------------------------------------------------------------------------
88
89
90end module no_spinup
Note: See TracBrowser for help on using the repository browser.