Changeset 1356 for trunk/LMDZ.TITAN/libf/phytitan/newstart.F

Timestamp:

Oct 8, 2014, 9:26:28 AM (10 years ago)

Author:

slebonnois

Message:

SL: update to newstart/start2archive tools in Venus+Titan / additional diagnostics in radiative fluxes for Titan

File:

• trunk/LMDZ.TITAN/libf/phytitan/newstart.F (modified) (8 diffs)

trunk/LMDZ.TITAN/libf/phytitan/newstart.F

@@ -22 +22 @@
       USE infotrac
       use cpdet_mod, only: ini_cpdet,t2tpot
+      use exner_hyb_m, only: exner_hyb
+      use exner_milieu_m, only: exner_milieu
 
       implicit none
@@ -135 +137 @@
       integer, dimension(4) :: start,counter
       REAL phisinverse(iip1,jjp1)  ! geopotentiel au sol avant inversion
-      logical topoflag,albedoflag
+      logical topoflag,albedoflag,razvitu,razvitv 
       real    albedo
       
@@ -969 +971 @@
 c--------------------------------------------------------------
 
+!!! ATTENTION TEMPORAIRE
+c     ps(:,:)=146700.
+
       ptotal =  0.
       DO j=1,jjp1
@@ -1004 +1009 @@
       CALL pression(ip1jmp1, ap, bp, ps, p3d)
          if (disvert_type==1) then
-           CALL exner_hyb(  ip1jmp1, ps, p3d,alpha,beta,pks, pk, pkf )
+           CALL exner_hyb(  ip1jmp1, ps, p3d, pks, pk, pkf )
          else ! we assume that we are in the disvert_type==2 case
-           CALL exner_milieu( ip1jmp1, ps, p3d, beta, pks, pk, pkf )
+           CALL exner_milieu( ip1jmp1, ps, p3d, pks, pk, pkf )
          endif
       
@@ -1012 +1017 @@
 
 c ATTENTION: peut servir, mais bon...
-c modif: profil uniforme
 c     do l=1,lmold
 c      do j=1,jmold+1
 c       do i=1,imold+1
+c modif: profil uniforme
 c          told(i,j,l)=told(1,jmold/2,l)
+c mean T profile:
+c       told(i,j,l) = 142.1*exp(-((p3d(i,j,l)/100.+21.45)/40.11)**2.) 
+c    .              + 106.3*exp(-((p3d(i,j,l)/100.-3183.)/4737.)**2.)
 c       enddo
 c      enddo
@@ -1040 +1048 @@
 c on assure la periodicite
       teta(iip1,:,:) =  teta(1,:,:)
+
+! RESETING U TO 0: may be done through run.def
+       razvitu = . FALSE .
+       CALL getin('razvitu',razvitu)
+       razvitv = . FALSE .
+       CALL getin('razvitv',razvitv)
 
 c calcul des champ de vent; passage en vent covariant
@@ -1057 +1071 @@
      &                   rlonuold,rlatvold,rlonu,rlatv)
       call scal_wind(us,vs,unat,vnat)
+! Reseting u=0
+      if (razvitu) then
+           unat(:,:,:) = 0.
+      endif
       write (*,*) 'unat ', unat (1,2,1)    ! INFO
       do l=1,llm
@@ -1070 +1088 @@
       write (*,*) 'ucov ', ucov (1,2,1)  ! INFO
 c     write(48,*) 'ucov',ucov
+! Reseting v=0
+      if (razvitv) then
+           vnat(:,:,:) = 0.
+      endif
+      write (*,*) 'vnat ', vnat (1,2,1)    ! INFO
       do l=1,llm
         do j = 1, jjm