Changeset 93 for trunk/SOURCES/climat-perturb_mod-0.4.f90

Timestamp:

11/28/16 16:03:50 (8 years ago)

Author:

dumas

Message:

First version with Schoof flux parameterisation at the grounding line. | New module furst_schoof_mod.f90 | New flag Schoof in grdline namelist (see in SOURCES/Fichiers-parametres/A-LBq15_param_list_Schoof.dat)

File:

• trunk/SOURCES/climat-perturb_mod-0.4.f90 (modified) (3 diffs)

trunk/SOURCES/climat-perturb_mod-0.4.f90

@@ -14 +14 @@
 
 use module3d_phy,only:nx,ny,S,S0,Tann,Tjuly,precip,acc,Ylat,num_forc,num_param,num_rep_42,dirforcage,dirnameinp,tafor,time,sealevel,coefbmshelf
+use netcdf
+use io_netcdf_grisli
 
 implicit none
@@ -41 +43 @@
 !! Routine qui permet d'initialiser les variations temporelles des variables climatiques
 !>
-  subroutine input_clim !routine qui permet d'initialiser les variations temporelles des variables climatiques
+
+  subroutine input_clim
 
     implicit none
+                character(len=100) :: precip_file        ! fichier precipitations
+    character(len=100) :: temp_annual_file   ! fichier temperature annuelle
+    real               :: coef_dens          ! pour corriger si donnees en eq. eau
+    logical            :: temp_param         ! si utilisation de temperature parametree
+    real*8, dimension(:,:), pointer :: data_2D => null() ! donnees lues dans le netcdf
+    
+    
     character(len=8) :: control      !label to check clim. forc. file (filin) is usable
     character(len=80):: filin
     integer ::  err                       !< pour l'allocation des tableaux
-    integer :: i
+    integer :: i,j
+    
+    namelist/climat_acc_T_gen/precip_file,coef_dens,temp_annual_file
+
+428 format(A)
+    rewind(num_param)                     ! pour revenir au debut du fichier param_list.dat
+    read(num_param,climat_acc_T_gen)
+
+    write(num_rep_42,428)'!___________________________________________________________' 
+    write(num_rep_42,428)'!  module  lect_clim_acc_T_ant_gen                          '
+    write(num_rep_42,climat_acc_T_gen)
+    write(num_rep_42,428)'!___________________________________________________________' 
+
+
+    ! precipitation
+    precip_file  = trim(dirnameinp)//trim(precip_file)
+    
+    call Read_ncdf_var('precip',trim(precip_file),data_2D)
+    precip(:,:)=data_2D(:,:)
+    !call lect_datfile(nx,ny,precip,1,precip_file)                  
+
+    precip(:,:)=precip(:,:)*coef_dens
+    acc(:,:)=precip(:,:)
+
+    if ((trim(temp_annual_file).eq.'no').or.(trim(temp_annual_file).eq.'NO')) then 
+       temp_param=.true.
+    else
+       temp_param=.false.
+    end if
+
+    !    temperature en surface 
+
+    test_param: if (.not.temp_param) then
+       temp_annual_file = trim(dirnameinp)//trim(temp_annual_file)
+
+
+                         call Read_ncdf_var('Tann',trim(temp_annual_file),data_2D)
+                         Tann(:,:)=data_2D(:,:)
+!       call lect_input(3,'Tann',1,Tann,temp_annual_file,trim(dirnameinp)//trim(runname)//'.nc')
+       !call lect_datfile(nx,ny,Tann,1,temp_annual_file)               ! temperature annuelle
+
+    else                        !    parametrisation de Fortuin pour la temperature annuelle.
+
+       do j=1,ny
+          do i=1,nx
+
+7            if (s0(i,j).le.200.) then                                    ! shelfs
+                tann(i,j)=49.642-0.943*abs(ylat(i,j))
+             else if ((s0(i,j).gt.200.).and.(s0(i,j).lt.1500.)) then      ! pente
+                tann(i,j)=36.689-0.005102*s0(i,j)-0.725*abs(ylat(i,j))
+             else if (s0(i,j).ge.1500.) then                              ! plateau
+                tann(i,j)=7.405-0.014285*s0(i,j)-0.180*abs(ylat(i,j))
+             endif
+          end do
+       end do
+    end if test_param
+
+    ta0(:,:)=tann(:,:)
+
+
+    !           pour la temperature d'ete, idem parametrisation huybrechts
+    do j=1,ny
+       do i=1,nx
+
+          tjuly(i,j)=tann(i,j)-17.65+0.00222*s0(i,j)&
+               +0.40802*abs(ylat(i,j))
+       end do
+    end do
+
+
+
+
 
     ! Lecture du forcage
@@ -138 +219 @@
 
 
-!!!!!!!! ATTENTION AJOUTE POUR TEST MAIS A REMETTRE AU PROPRE PLUS TARD C. DUMAS !!!!!!!!!!
-!!!!!!!! ancien input_climat_ref de lect_clim_act_anteis
-!     accumulation de Philippe
-      filin='accumHUY40km.dat'
-      call lect_eis(nx,ny,precip,filin,DIRNAMEINP)
-!====================================== La reponse est 42 ===========
-      write(num_rep_42,*) 'fichier accum : ', filin 
-
-!     cas particulier de Vostok
-      ivo=101
-      jvo=62
-      do j=jvo-1,jvo+1
-        do i=ivo-1,ivo+1
-            precip(i,j)=0.02         ! valeur plus faible a Vostok.
-        end do
-      end do
-      acc(:,:)=precip(:,:)
-
-!    temperature en surface :
-!    parametrisation de Fortuin pour la temperature annuelle.
-      do j=1,ny
-        do i=1,nx
-
-            if (s0(i,j).le.200.) then   ! shelfs
-               tann(i,j)=49.642-0.943*abs(ylat(i,j))
-            else if ((s0(i,j).gt.200.).and.(s0(i,j).lt.1500.)) then ! pente
-               tann(i,j)=36.689-0.005102*s0(i,j)-0.725*abs(ylat(i,j))
-            else if (s0(i,j).ge.1500.) then        ! plateau
-               tann(i,j)=7.405-0.014285*s0(i,j)-0.180*abs(ylat(i,j))
-            endif
-
-            ta0(i,j)=tann(i,j)
-!           pour la temperature d'ete, idem parametrisation huybrechts
-            tjuly(i,j)=tann(i,j)-17.65+0.00222*s0(i,j)&
-                         +0.40802*abs(ylat(i,j))
-        end do
-      end do
-!!!!!!!! FIN MODIF  TEMPORAIRE !!!!!!!!!!
-
+!cdc Commente pour être compatible avec lecture fichiers Cat Schoofing
+
+!~ !!!!!!!! ATTENTION AJOUTE POUR TEST MAIS A REMETTRE AU PROPRE PLUS TARD C. DUMAS !!!!!!!!!!
+!~ !!!!!!!! ancien input_climat_ref de lect_clim_act_anteis
+!~ !     accumulation de Philippe
+!~       filin='accumHUY40km.dat'
+!~       call lect_eis(nx,ny,precip,filin,DIRNAMEINP)
+!~ !====================================== La reponse est 42 ===========
+!~       write(num_rep_42,*) 'fichier accum : ', filin 
+
+!~ !     cas particulier de Vostok
+!~       ivo=101
+!~       jvo=62
+!~       do j=jvo-1,jvo+1
+!~         do i=ivo-1,ivo+1
+!~             precip(i,j)=0.02         ! valeur plus faible a Vostok.
+!~         end do
+!~       end do
+!~       acc(:,:)=precip(:,:)
+
+!~ !    temperature en surface :
+!~ !    parametrisation de Fortuin pour la temperature annuelle.
+!~       do j=1,ny
+!~         do i=1,nx
+
+!~             if (s0(i,j).le.200.) then   ! shelfs
+!~                tann(i,j)=49.642-0.943*abs(ylat(i,j))
+!~             else if ((s0(i,j).gt.200.).and.(s0(i,j).lt.1500.)) then ! pente
+!~                tann(i,j)=36.689-0.005102*s0(i,j)-0.725*abs(ylat(i,j))
+!~             else if (s0(i,j).ge.1500.) then        ! plateau
+!~                tann(i,j)=7.405-0.014285*s0(i,j)-0.180*abs(ylat(i,j))
+!~             endif
+
+!~             ta0(i,j)=tann(i,j)
+!~ !           pour la temperature d'ete, idem parametrisation huybrechts
+!~             tjuly(i,j)=tann(i,j)-17.65+0.00222*s0(i,j)&
+!~                          +0.40802*abs(ylat(i,j))
+!~         end do
+!~       end do
+!~ !!!!!!!! FIN MODIF  TEMPORAIRE !!!!!!!!!!
+!cdc fin Commente pour être compatible avec lecture fichiers Cat Schoofing
 
 return