Changeset 435 for branches

Timestamp:

06/22/23 17:12:00 (12 months ago)

Author:

dumas

Message:

new output in output_antarcti_mod to track changes in Antarctica drilling sites | input file : Ice core list in icecore-Ant.dat and variables descriptions in icecore-var.dat

Location:

branches/GRISLIv3/SOURCES

Files:

• Ant40_files/output_anta40_mod-0.4.f90 (modified) (6 diffs)
• Fichiers-parametres/icecore-Ant.dat (added)
• Fichiers-parametres/icecore-var.dat (added)

branches/GRISLIv3/SOURCES/Ant40_files/output_anta40_mod-0.4.f90

@@ -62 +62 @@
 integer,dimension(nvar) :: varID
 integer :: nbtimeout  ! index time output
-
 real,dimension(nvar) :: var_shortoutput
 
+! variables netcdf icecore (ic)
+integer,parameter :: ncicecoreout=1 ! 1 sorties netcdf icecore / 0 no output
+integer :: nbic ! number of ice core site in output
+integer,parameter :: nvaric=10 ! nombre de variables dans le fichier ice core
+real,parameter :: dticecoreout=50. ! sorties ice core tous les dticecoreout
+integer :: ncidic
+integer :: timeDimIDic
+integer :: timeVarIDic
+integer :: stationDimIDic
+integer :: stationVarIDic
+integer :: station_strlenDimIDic
+integer :: lonVarIDic
+integer :: latVarIDic
+integer,dimension(nvaric) :: varIDic
+integer :: nbtimeoutic  ! index time output
+real,allocatable,dimension(:,:) :: var_icoutput
+character(len=100), allocatable, dimension(:) :: site_nameic
+integer, allocatable, dimension(:) :: iic,jic ! coordonnees i,j des sites ice core
 
 CONTAINS
@@ -75 +92 @@
 character(len=100),dimension(nvar) :: long_name
 character(len=100),dimension(nvar) :: units
+
+! ice core variables
+character(len=100),allocatable,dimension(:) :: ic_name ! ice core name
+real,allocatable,dimension(:) :: ic_lon ! ice core lon
+real,allocatable,dimension(:) :: ic_lat ! ice core lat
+character(len=100),dimension(nvaric) :: namevaric ! name, standard_name, long_name, units
+character(len=100),dimension(nvaric) :: standard_nameic
+character(len=100),dimension(nvaric) :: long_nameic
+character(len=100),dimension(nvaric) :: unitsic
+integer :: err ! recuperation erreur
+integer :: ii,jj
+integer :: nloop ! nbr de tour dans la boucle while (evite boucle folle)
+real :: dist_lonlat, dist_lonlat_new, distance ! distance to lon lat ice core site
+
+
 
 !ndisp sorite courte tous les ndisp
@@ -106 +138 @@
   enddo
   close(568)
+  
 ! Fichier Netcdf initMIP
 ! creation du fichier Netcdf :
@@ -136 +169 @@
   enddo
 
-! fin de la definition du fchier :
+! fin de la definition du fichier :
   status=nf90_enddef(ncid)
   if (status /= nf90_noerr) call handle_err(status)
@@ -143 +176 @@
   corrsurf(:,:)=1.
 endif
+
+
+if (ncicecoreout.eq.1) then  ! ecriture netcdf Ice Core
+! Fichier Netcdf ice core Antarctique
+! Lecture du fichier avec la liste des sites de forages
+! infos sur les variables :
+  open(568,file=trim(dirsource)//'/Fichiers-parametres/icecore-Ant.dat',status='old')
+  ! lecture en-tete
+  read(568,*)
+  read(568,*)
+  read(568,*)
+  read(568,*)
+  read(568,*)
+  ! lectuire du nombre de sites de forages:
+  read(568,*) nbic
+  read(568,*)
+  ! allocation des variables dependant du nombre de site de forage (nbic)
+  if (.not.allocated(iic)) then
+     allocate(iic(nbic),stat=err)
+     if (err/=0) then
+        print *,"erreur a l'allocation du tableau iic dans init_outshort",err
+        stop 4
+     end if
+  end if
+  if (.not.allocated(jic)) then
+     allocate(jic(nbic),stat=err)
+     if (err/=0) then
+        print *,"erreur a l'allocation du tableau jic dans init_outshort",err
+        stop 4
+     end if
+  end if
+  if (.not.allocated(ic_name)) then
+     allocate(ic_name(nbic),stat=err)
+     if (err/=0) then
+        print *,"erreur a l'allocation du tableau ic_name dans init_outshort",err
+        stop 4
+     end if
+  end if
+  if (.not.allocated(ic_lon)) then
+     allocate(ic_lon(nbic),stat=err)
+     if (err/=0) then
+        print *,"erreur a l'allocation du tableau ic_lon dans init_outshort",err
+        stop 4
+     end if
+  end if
+  if (.not.allocated(ic_lat)) then
+     allocate(ic_lat(nbic),stat=err)
+     if (err/=0) then
+        print *,"erreur a l'allocation du tableau ic_lat dans init_outshort",err
+        stop 4
+     end if
+  end if
+  if (.not.allocated(var_icoutput)) then
+     allocate(var_icoutput(nbic,nvaric),stat=err)
+     if (err/=0) then
+        print *,"erreur a l'allocation du tableau var_icoutput dans init_outshort",err
+        stop 4
+     end if
+  end if
+! lecture des infos sur les variables :
+  do k=1,nbic
+     read(568,'(a100)') ic_name(k)
+     !ic_name(k)=trim(ic_name(k))
+     read(568,*) ic_lon(k)
+     read(568,*) ic_lat(k)
+     read(568,*)
+  enddo
+ close(568)
+ ! recherche du point i,j correspondant aux coordonnées lon,lat de chaque site de forage
+ do k=1,nbic
+    dist_lonlat = 1.e7 ! initialisation
+    dist_lonlat_new = 1.e6 ! initialisation
+    ! on part du centre de la grille
+!    jj=int(ny/2)
+!    ii=int(nx/2)
+    jj=6
+    ii=6
+    nloop=0
+    do while (dist_lonlat_new .LT. dist_lonlat .AND. nloop .LT. 10000)
+       nloop=nloop+1
+       dist_lonlat = dist_lonlat_new
+       ii = iic(k)
+       jj = jic(k)
+       do j=max(jj-5,1),min(jj+5,ny)
+          do i=max(ii-5,1),min(ii+5,nx)
+             distance = abs(xlong(i,j) - ic_lon(k)) + abs(ylat(i,j) - ic_lat(k))
+             if (distance .LT. dist_lonlat) then
+                dist_lonlat_new = distance
+                iic(k) = i
+                jic(k) = j
+             endif
+          enddo
+       enddo
+    enddo
+    print*,'init_outshort : ice core ', trim(ic_name(k)), nloop
+    print*,'init_outshort lon,lat', ic_lon(k),ic_lat(k)
+    print*,'init_outshort i,j',iic(k),jic(k)
+    print*,'init_outshort lon lat GRISLI',xlong(iic(k),jic(k)),ylat(iic(k),jic(k))
+ enddo           
+          
+
+
+! creation du fichier Netcdf :
+  status=nf90_create(path = trim(dirnameout)//'icecore'//runname//'.nc', cmode = nf90_clobber, ncid = ncidic)
+  if (status /= nf90_noerr) call handle_err(status)
+
+  ! definition des dimension :
+  status=nf90_def_dim(ncidic, name="time", len=NF90_UNLIMITED, dimid=timeDimIDic)
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_def_var(ncidic, name="time", xtype=nf90_float, dimids=(/ timeDimIDic/), varid=timeVarIDic) 
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_put_att(ncidic, timeVarIDic, "standard_name", "time")
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_put_att(ncidic, timeVarIDic,"units", "years since 0000-01-01 00:00:00")
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_put_att(ncidic, timeVarIDic,"calendar", "360_day")
+  if (status /= nf90_noerr) call handle_err(status)
+  
+  ! dimension "station" Sites de forages
+  status=nf90_def_dim(ncidic, name="station", len=nbic, dimid=stationDimIDic)
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_def_dim(ncidic, name="name_strlen", len=len(ic_name(1)), dimid=station_strlenDimIDic)
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_def_var(ncidic, name="station_name", xtype=nf90_char, dimids=(/ station_strlenDimIDic, stationDimIDic /), varid=stationVarIDic) 
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_put_att(ncidic, stationVarIDic, "standard_name", "station_name")
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_put_att(ncidic, stationVarIDic, "long_name", "station name")
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_def_var(ncidic, name="lon", xtype=nf90_float, dimids=(/ stationDimIDic/), varid=lonVarIDic) 
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_put_att(ncidic, lonVarIDic, "standard_name", "longitude")
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_put_att(ncidic, lonVarIDic, "long_name", "station longitude")
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_def_var(ncidic, name="lat", xtype=nf90_float, dimids=(/ stationDimIDic/), varid=latVarIDic) 
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_put_att(ncidic, latVarIDic, "standard_name", "latitude")
+  if (status /= nf90_noerr) call handle_err(status)
+  status=nf90_put_att(ncidic, latVarIDic, "long_name", "station latitude")
+  if (status /= nf90_noerr) call handle_err(status)
+
+  
+! infos sur les variables :
+  open(568,file=trim(dirsource)//'/Fichiers-parametres/icecore-var.dat',status='old')
+! lecture en-tete
+  read(568,*)
+  read(568,*)
+  read(568,*)
+  read(568,*)
+  read(568,*)
+! lecture des infos sur les variables :
+  do l=1,nvaric
+    read(568,'(a100)') namevaric(l)
+    read(568,'(a100)') standard_nameic(l)
+    read(568,'(a100)') long_nameic(l)
+    read(568,'(a100)') unitsic(l)
+    read(568,*)
+ enddo
+ close(568)
+  
+! definition des variables de sortie :
+  do l=1,nvaric  ! boucle sur le nbr de variable a definir
+    status=nf90_def_var(ncidic, name=trim(namevaric(l)), xtype=nf90_float, dimids= &
+          (/ stationDimIDic, timeDimIDic /), varid=varIDic(l))
+    if (status /= nf90_noerr) call handle_err(status)
+    status=nf90_put_att(ncidic, varIDic(l), "standard_name", trim(standard_nameic(l)))
+    if (status /= nf90_noerr) call handle_err(status)
+    status=nf90_put_att(ncidic, varIDic(l), "long_name", trim(long_nameic(l)))
+    if (status /= nf90_noerr) call handle_err(status)
+    status=nf90_put_att(ncidic, varIDic(l), "units", trim(unitsic(l)))
+    if (status /= nf90_noerr) call handle_err(status)
+    status=nf90_put_att(ncidic, varIDic(l), "coordinates", "lon lat")
+    if (status /= nf90_noerr) call handle_err(status)
+  enddo
+
+! fin de la definition du fichier :
+  status=nf90_enddef(ncidic)
+  if (status /= nf90_noerr) call handle_err(status)
+  
+  nbtimeoutic = 0 ! initialisation compteur sorties axe time
+  ! ecriture des noms des sites de forage dans la variable site :
+!  do k=1,nbic
+!     status=nf90_put_var(ncidic, stationVarIDic, trim(ic_name(k)),start=(/k/),count=(/1/))
+!     if (status /= nf90_noerr) call handle_err(status)
+!  enddo
+  print*,'ic_name',ic_name   
+  status=nf90_put_var(ncidic, stationVarIDic, ic_name) !,start=(/1/),count=(/nbic/))
+  if (status /= nf90_noerr) call handle_err(status)
+
+  status=nf90_put_var(ncidic, lonVarIDic, ic_lon,start=(/1/),count=(/nbic/))
+  if (status /= nf90_noerr) call handle_err(status)
+     
+  status=nf90_put_var(ncidic, latVarIDic, ic_lat,start=(/1/),count=(/nbic/))
+  if (status /= nf90_noerr) call handle_err(status)
+!     print*, 'debug 2 init_outshort'
+endif
+
+
 
 end subroutine init_outshort
@@ -319 +551 @@
 
 
+! ecriture netcdf ice core output (tous les 50 ans)
+if (ncicecoreout.eq.1 .and. (mod(abs(time),dticecoreout).lt.dtmin)) then  ! ecriture netcdf
+   nbtimeoutic=nbtimeoutic+1
+   do k=1,nbic ! boucle sur les stations
+      ! Surface elevation
+      var_icoutput(k,1)=S(iic(k),jic(k))
+      ! Ice thickness
+      var_icoutput(k,2)=H(iic(k),jic(k))
+      ! Bedrock elevation
+      var_icoutput(k,3)=Bsoc(iic(k),jic(k))
+      ! Surface Mass Balance
+      var_icoutput(k,4)=Bm(iic(k),jic(k))
+      ! Basal melt
+      var_icoutput(k,5)=Bmelt(iic(k),jic(k))
+      ! Surface velocity in x
+      var_icoutput(k,6)=Uxbar(iic(k),jic(k))
+      ! Surface velocity in y
+      var_icoutput(k,7)=Uybar(iic(k),jic(k))
+      ! Velocity magnitude
+      var_icoutput(k,8)=sqrt(Uxbar(iic(k),jic(k))**2 + Uybar(iic(k),jic(k))**2)
+      ! Surface Temperature
+      var_icoutput(k,9)=Ts(iic(k),jic(k))
+      ! Accumulation
+      var_icoutput(k,10)=Acc(iic(k),jic(k))
+   enddo
+  ! ecriture de time
+   status=nf90_put_var(ncidic, timeVarIDic, time, start=(/nbtimeoutic/))
+   if (status /= nf90_noerr) call handle_err(status)
+   
+   do l=1,nvaric  ! boucle sur le nbr de variables a ecrire
+      status=nf90_put_var(ncidic, varIDic(l), var_icoutput(:,l),start=(/1,nbtimeoutic/))
+      if (status /= nf90_noerr) call handle_err(status)
+   enddo
+   status=nf90_sync(ncidic)
+   if (status /= nf90_noerr) call handle_err(status)
+endif
+
+
+
 end subroutine shortoutput