Changeset 1130 for trunk/LMDZ.MARS/libf/phymars/newstart.F

Timestamp:

Dec 20, 2013, 4:04:56 PM (11 years ago)

Author:

emillour

Message:

Mars GCM:
Series of changes to enable running in parallel (using LMDZ.COMMON dynamics);
Current LMDZ.MARS can still notheless be compiled and run in serial mode
"as previously".
Summary of main changes:

• Main programs (newstart, start2archive, xvik) that used to be in dyn3d have been moved to phymars.
• dyn3d/control.h is now module control_mod.F90
• rearanged input/outputs routines everywhere to handle serial/MPI cases. physdem.F => phyredem.F90 , phyetat0.F => phyetat0.F90 ; all read/write routines for startfi files are gathered in module iostart.F90
• added parallelism related routines init_phys_lmdz.F90, comgeomphy.F90, dimphy.F90, iniphysiq.F90, mod_grid_phy_lmdz.F90, mod_phys_lmdz_mpi_data.F90, mod_phys_lmdz_mpi_transfert.F90, mod_phys_lmdz_omp_data.F90, mod_phys_lmdz_omp_transfert.F90, mod_phys_lmdz_para.F90, mod_phys_lmdz_transfert_para.F90 in phymars and mod_const_mpi.F90 in dyn3d (for compliance with parallel case)
• created generic routines 'planetwide_maxval' and 'planetwide_minval', in module "planetwide_mod", that enable obtaining the min and max of a field over the whole planet.

EM

File:

• trunk/LMDZ.MARS/libf/phymars/newstart.F (moved) (moved from trunk/LMDZ.MARS/libf/dyn3d/newstart.F) (16 diffs)

trunk/LMDZ.MARS/libf/phymars/newstart.F

@@ -17 +17 @@
 ! to use  'getin'
       use ioipsl_getincom, only: getin 
-      use infotrac, only: iniadvtrac, nqtot, tnom
-      use tracer_mod, only: noms, igcm_h2o_vap, igcm_h2o_ice
+      use infotrac, only: iniadvtrac, nqtot, tname
+      use tracer_mod, only: noms, igcm_dust_number, igcm_dust_mass,
+     &                      igcm_ccn_number, igcm_ccn_mass,
+     &                      igcm_h2o_vap, igcm_h2o_ice
       use surfdat_h, only: phisfi, z0, zmea, zstd, zsig, zgam, zthe,
      &                     albedodat, z0_default
       use comsoil_h, only: inertiedat, layer, mlayer, nsoilmx
+      use control_mod, only: day_step, iphysiq, anneeref
+      use phyredem, only: physdem0, physdem1
+      use iostart, only: open_startphy
+      use comgeomphy, only: initcomgeomphy
       implicit none
 
@@ -36 +42 @@
 #include "comvert.h"
 #include "comgeom2.h"
-#include "control.h"
+!#include "control.h"
 #include "logic.h"
 #include "description.h"
 #include "ener.h"
 #include "temps.h"
-#include "lmdstd.h"
+!#include "lmdstd.h"
 #include "comdissnew.h"
 #include "clesph0.h"
@@ -142 +148 @@
 c variables diverses
 c-------------------
-      real choix_1
+      real choix_1 ! ==0 : read start_archive file ; ==1: read start files
       character*80      fichnom
       integer Lmodif,iq
@@ -183 +189 @@
       preff  = 610.    ! for Mars, instead of 101325. (Earth)
       pa= 20           ! for Mars, instead of 500 (Earth)
+
+! initialize "serial/parallel" related stuff
+      CALL init_phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
+      call initcomgeomphy
 
 ! Load tracer number and names:
@@ -316 +326 @@
 c-----------------------------------------------------------------------
       if (choix_1.eq.0) then
+         ! tabfi requires that input file be first opened by open_startphy(fichnom)
+         fichnom = 'start_archive.nc'
+         call open_startphy(fichnom)
          call tabfi (nid,Lmodif,tab0,day_ini,lllm,p_rad,
      .            p_omeg,p_g,p_mugaz,p_daysec,time)
       else if (choix_1.eq.1) then
+         fichnom = 'startfi.nc'
+         call open_startphy(fichnom)
          call tabfi (nid_fi,Lmodif,tab0,day_ini,lllm,p_rad,
      .            p_omeg,p_g,p_mugaz,p_daysec,time)
@@ -461 +476 @@
         txt=" "
         write(txt,'(a1,i2.2)') 'q',iq
-        if (txt.eq.tnom(iq)) then
+        if (txt.eq.tname(iq)) then
           count=count+1
         endif
@@ -471 +486 @@
       if (count.eq.nqtot) then
         write(*,*) 'Newstart: updating tracer names'
-        ! copy noms(:) to tnom(:) to have matching tracer names in physics
+        ! copy noms(:) to tname(:) to have matching tracer names in physics
         ! and dynamics
-        tnom(1:nqtot)=noms(1:nqtot)
+        tname(1:nqtot)=noms(1:nqtot)
       endif
 
@@ -497 +512 @@
      $ tracer'
       write(*,*) 'q=profile    : specify a profile for a tracer'
+      write(*,*) 'freedust     : rescale dust to a true value'
       write(*,*) 'ini_q        : tracers initialization for chemistry
      $ and water vapour'
@@ -723 +739 @@
           write(*,*) 'Which tracer name do you want to change ?'
           do iq=1,nqtot
-            write(*,'(i3,a3,a20)')iq,' : ',trim(tnom(iq))
+            write(*,'(i3,a3,a20)')iq,' : ',trim(tname(iq))
           enddo
           write(*,'(a35,i3)')
@@ -730 +746 @@
           read(*,*) iq
           if ((iq.ge.1).and.(iq.le.nqtot)) then
-            write(*,*)'Change tracer name ',trim(tnom(iq)),' to ?'
+            write(*,*)'Change tracer name ',trim(tname(iq)),' to ?'
             read(*,*) txt
-            tnom(iq)=txt
+            tname(iq)=txt
             write(*,*)'Do you want to change another tracer name (y/n)?'
             read(*,'(a)') yes 
@@ -768 +784 @@
              write(*,*) 'Which tracer do you want to modify ?'
              do iq=1,nqtot
-               write(*,*)iq,' : ',trim(tnom(iq))
+               write(*,*)iq,' : ',trim(tname(iq))
              enddo
              write(*,*) '(choose between 1 and ',nqtot,')'
@@ -777 +793 @@
                cycle
              endif
-             write(*,*)'mixing ratio of tracer ',trim(tnom(iq)),
+             write(*,*)'mixing ratio of tracer ',trim(tname(iq)),
      &                 ' ? (kg/kg)'
              read(*,*) val
@@ -787 +803 @@
                ENDDO
              ENDDO
-             write(*,*) 'SURFACE value of tracer ',trim(tnom(iq)),
+             write(*,*) 'SURFACE value of tracer ',trim(tname(iq)),
      &                   ' ? (kg/m2)'
              read(*,*) val
@@ -804 +820 @@
              write(*,*) 'Which tracer do you want to set?'
              do iq=1,nqtot
-               write(*,*)iq,' : ',trim(tnom(iq))
+               write(*,*)iq,' : ',trim(tname(iq))
              enddo
              write(*,*) '(choose between 1 and ',nqtot,')'
@@ -814 +830 @@
              endif
              ! look for input file 'profile_tracer'
-             txt="profile_"//trim(tnom(iq))
+             txt="profile_"//trim(tname(iq))
              open(41,file=trim(txt),status='old',form='formatted',
      &            iostat=ierr)
@@ -831 +847 @@
                    q(:,:,l,iq)=profile(l+1)
                  enddo
-                 write(*,*)'OK, tracer ',trim(tnom(iq)),' initialized ',
-     &                     'using values from file ',trim(txt)
+                 write(*,*)'OK, tracer ',trim(tname(iq)),
+     &               ' initialized ','using values from file ',trim(txt)
                else
                  write(*,*)'problem reading file ',trim(txt),' !'
-                 write(*,*)'No modifications to tracer ',trim(tnom(iq))
+                 write(*,*)'No modifications to tracer ',trim(tname(iq))
                endif
              else
                write(*,*)'Could not find file ',trim(txt),' !'
-               write(*,*)'No modifications to tracer ',trim(tnom(iq))
+               write(*,*)'No modifications to tracer ',trim(tname(iq))
              endif
              
+c       q=profile : initialize tracer with a given profile
+c       --------------------------------------------------
+        else if (trim(modif) .eq. 'freedust') then
+          do l=1,llm
+            do j=1,jjp1
+              do i=1,iip1
+                if (igcm_dust_number .ne. 0) 
+     &          q(i,j,l,igcm_dust_number)=
+     &                q(i,j,l,igcm_dust_number) * 1e-3 ! grosso modo
+                if (igcm_dust_mass .ne. 0) 
+     &          q(i,j,l,igcm_dust_mass)=
+     &                q(i,j,l,igcm_dust_mass)   * 1e-3 ! grosso modo
+                if (igcm_ccn_number .ne. 0) 
+     &          q(i,j,l,igcm_ccn_number)=
+     &                q(i,j,l,igcm_ccn_number) * 1e-3 ! grosso modo
+                if (igcm_ccn_mass .ne. 0) 
+     &          q(i,j,l,igcm_ccn_mass)=
+     &                q(i,j,l,igcm_ccn_mass)   * 1e-3 ! grosso modo
+              end do
+            end do
+          end do
+
+         ! We want to have the very same value at lon -180 and lon 180
+          do l = 1,llm
+             do j = 1,jjp1
+                do iq = 1,nqtot
+                   q(iip1,j,l,iq) = q(1,j,l,iq)
+                end do
+             end do
+          end do
 
 c       ini_q : Initialize tracers for chemistry