Changeset 248 for trunk/SOURCES/climat_forcage_mod.f90

Timestamp:

01/28/19 14:06:08 (5 years ago)

Author:

dumas

Message:

climat_forcage_mod reads the topo corresponding to the reference climate in clim_ref_file

File:

• trunk/SOURCES/climat_forcage_mod.f90 (modified) (7 diffs)

trunk/SOURCES/climat_forcage_mod.f90

@@ -6 +6 @@
 ! possibilite d'utiliser autant de snapshots que l'on veut (ntr)
 ! fichiers de forcage, Snapshot et valeurs de lapse rate definis dans fichier param
-  use module3d_phy,only:nx,ny,sealevel,sealevel_2d,S,S0,Tmois,Tann,Tjuly,acc,num_forc,num_param,num_rep_42,tafor,coefbmshelf,PI,ro,time,dirforcage,dirnameinp
+  use module3d_phy,only:nx,ny,sealevel,sealevel_2d,S,Tmois,Tann,Tjuly,acc,num_forc,num_param,num_rep_42,tafor,coefbmshelf,PI,ro,time,dirforcage,dirnameinp
   use netcdf
   use io_netcdf_grisli
@@ -41 +41 @@
                      ! 2 for steady state using snap 2 ...
                      ! 0 for transient
-
-! S0 topo de ref correspondant a la climato est lu dans lect_topo 
   
 ! variables pour climato mensuelle :  
   real,dimension(nx,ny,12) :: t2m0   !< initial monthly air temperature (climato)
   real,dimension(nx,ny,12) :: pr0    !< initial monthly precip (climato)
+  real,dimension(nx,ny) :: orog0     !< topo de la climato
   
 ! climat des snapshots :
@@ -53 +52 @@
   real,dimension(nx,ny,ntr) :: orog_ss    !< topo for every GCM snapshot 
   
-  real,dimension(nx,ny,12,ntr) :: t2m_sstopo0 !< t2m GCM snapshot on S0
-  real,dimension(nx,ny,12,ntr) :: pr_sstopo0  !< pr GCM snapshot on S0
+  real,dimension(nx,ny,12,ntr) :: t2m_sstopo0 !< t2m GCM snapshot on orog0
+  real,dimension(nx,ny,12,ntr) :: pr_sstopo0  !< pr GCM snapshot on orog0
   
   real :: PYY                        !< constante pour calcul de la fraction solide des precips
@@ -87 +86 @@
   call Read_Ncdf_var('t2m',trim(DIRNAMEINP)//TRIM(clim_ref_file),tab3d)
   t2m0(:,:,:)=tab3d(:,:,:)
+  call Read_Ncdf_var('orog',trim(DIRNAMEINP)//TRIM(clim_ref_file),tab2d)
+  orog0(:,:)=tab2d(:,:)
   
   filtr(1)=TRIM(DIRNAMEINP)//'Snapshots-GCM/'//TRIM(forcage_file1)
@@ -125 +126 @@
   do m=1,12
     do k=1,ntr
-      t2m_sstopo0(:,:,m,k)=t2m_ss(:,:,m,k)-lapserate*(S0(:,:)-orog_ss(:,:,k))
+      t2m_sstopo0(:,:,m,k)=t2m_ss(:,:,m,k)-lapserate*(orog0(:,:)-orog_ss(:,:,k))
       pr_sstopo0(:,:,m,k)=pr_ss(:,:,m,k)*exp(rappact*(t2m_sstopo0(:,:,m,k)-t2m_ss(:,:,m,k)))
     enddo
@@ -246 +247 @@
   integer :: ift                     !< indice correspondant au pas de temps du fichier de forcage
   real :: TEMP                       !< calcul du nbr de jour < psolid
-  real,dimension(nx,ny) :: deltaZs   ! diff temp par rapport a topo S0 
+  real,dimension(nx,ny) :: deltaZs   ! diff temp par rapport a topo orog0 
   real,dimension(nx,ny,12) :: prmois ! precip sur topo Zs
   
@@ -350 +351 @@
   !$OMP WORKSHARE
   Zs(:,:)=max(sealevel_2d(:,:),S(:,:))
-  deltaZs(:,:)= -lapserate*(Zs(:,:)-S0(:,:)) ! difference de temperature entre Zs et S0
+  deltaZs(:,:)= -lapserate*(Zs(:,:)-orog0(:,:)) ! difference de temperature entre Zs et orog0
   !$OMP END WORKSHARE
 ! correction d'altitude